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16754 POST- BUCK NG. ANA ALY OF GUYED TOWERS 


——-— ar WORDS: Buckling; Cables; Guy wires; Horizontal loads; Poles; 


= "ABSTRACT: The p post- buckling behavior of short towers, supported by a sinaie.s set ra 
| guy wires is investigated. Using the principle of stationary potential energy, a set f 


' 
| ‘non- linear algebraic equations are obtained for both fixed and hinged-base towers. The 
| equations are solved numerically for several actual systems. The results of the 
iD investigation include the shape of the post-buckling curves for towers subjected to 


_ vertical loads, the effect of cable pretension on the post-buckling behavior and - 
effect of horizontal on n the ne capacity of the towers | to ‘supprot vertical loads. 
Massachusetts, Amherst, Mass. 01003), and Ling, David, “Post- -Buckling Analysis of — 
| Bose Towers,” Journal of the Structural Division, ASCE, Vol. 107, No. ST12, Proc. | 


16768 NONLINEAR DYNAMIC ANALYSIS BY MODAL gi 0 wewssioy or 


KEY WORDS: Analysis; Approximation method; Cables; ‘Dynamic 
‘yespense; Dynamics; Structural engineering; Superposition (mathematics); 


ABSTRACT: modified superposition is developed and used to 
evaluate the nonlinear dynamic response of a cable network composed of a linkage of 
 euuighe members with lumped masses. The method is compared with Newmark’s B 


method for accuracy and the time interval requirements for stability. To show th 
usefulness of the method a number of parameter studies are presented using the cable 


REFERENCE: Geschwindner, Louis F., Jr. (Asst. Ail Architectural Engrg., 
Pennsylvania State Univ., University Park, Pa.), “Nonlinear Dynamic Analysis by 


Modal Superposition,” Journal of the Structural Division, ASCE, Vol. 107, No. ST12, 


roc. - Paper 16768, December, 1981, pp. 2325-2336 


16771 OPTIMIZATION WITH FREQUENCY CONSTRAINT 


KEY WORDS: "Allowable stress; Eigenvalues; Frequency _ 


ABSTRACT: The us use of ‘allowable stress” algorithms for structural 
optimization problem with a single frequency constraint and a constant mass matrix is 

| discussed. The mathematical properties of the eigenvalue problem are reviewed. An 

_ Allowable stress algorithm is stated and its properties discussed. Two simple examples — 
are presented. The primary concern of this paper lies in placing allowable stress- — 


REFERENCE: ‘Spillers, V WwW William R. (Prof., Dept. of Cv. -Engrg., Rensselaer 
Polytechnic Inst., Troy, N.Y. 12181), Singh, Sohan, and Levy, Robert, “Optimization 
with Frequency Constraint,” Journal of the Structural Division, ASCE, Vol. 107, No. 


ST12, Proc. Paper 16771, December, 
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OLD WOOD JOISTS 


Ww ORDS: Buildings; Joists; ; Nondestructi tests; Reinforceme 
Statistical analysis; Strength of materials; Stress waves; tee 


xt 
Timber construction; Visual inspection __ 


_ ABSTRACT: Two nondestructive techniques are used to evaluate the strength of wood | 
joists in an existing structure. The traditional visual grading and longitudinal stress 

} & approaches are employed, with the latter technique providing better insight into — 


the strength of the floor system. Further refinements of the stress wave method are © 


REFERENCE: Lanius, Ross M., Prof., Dept. of Civ. and Environmental 
H Bg.” Univ. of New Haven, Ww. Haven, Conn.), Tichy, Robert, and Bulleit, William © 
7, No. ST12, Proc. ‘Pa er 16767, December, 1981, P- 2349- 2363 
a 10 per | Pp 


16770 FRAMES UNDER WIND DISTURBANCES 


KEY W ORDS: Collapse; Dynamics; Frame design; Framed structures; 
Frames; Load factors; Numerical Parametric ‘equations; 
ABSTRACT: AY parametric numerical analysis is on strain- multi- 
story frames subjected to wind disturbances generated by amplifying a recorded 
fluctuating wind pressure. The effect of spreading and diminishing of the strain- 
hardening regions has been taken into account in the one dimensional finite element 
_ formulation employed here. The dynamic behavior including in particular the dynamic | 
to collapse behavior is extensively examined and the details are presented for the first 
i time. The effect of a fluctuating wind disturbance on a frame may be represented by 
the required overall plastic strength, the therefore measured by a multiplier (load — 
factor) of the basic design wind loads. There exists a critical frame which just exhibits 
the dynamic collaps behavior under a wind disturbance. The corresponding critical 
oles of the multiplier has been called the dynamic collapse gust load factor. The 
factor provides a rational basis for a practical safety factor and is therefore 
recommended for a more rational method of wind resistant design of ductile —— 


REFERENCE: Morisako, ‘Kiyotaka (Asst., Dept. of Architecture, Kyoto Tech. Univ., — 
Kyoto, Japan), Ishida, Shuzo, Nakamura, Tsuneyoshi, and Ishizaki, Hatsuo, “Dynamic 
Collapse of Frames Under Wind Disturbances,” Journal of the Structural Division, — 
ASCE, Vol. 107, No. ST12, Proc. Paper PP. 
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“16753 BUCKI CKLING OF REINFORCED | CONCRETE SHELLS wird 
ORDS: Concrete (reinforced); Domes (structural); 
Engineering mechanics; Loads (forces); Plastic properties; Shells (structura 
BSTRACT: The special phenomena of reinforced concrete shell buckling is described 
on the basis of the latest results of research. The method presented Fey oll 
a consideration the special properties of the reinforced concrete including, plasticity, 
creep, and quantity of reinforcement. Snap-through critical loads calculated by the 
described method are compared to the results of model tests and to the loads and — 

_ ¢ritical loads of erected domes. The shell buckling stiffness was seen to be composed of 
‘a two parts; the bending and the axial stiffness of the reinforced concrete cross section. Fe - 
The critical load is reduced by the plasticity and the imperfections. The critical load of 
4 he reinforced concrete shells decreases by the craks of concrete and increases by the 


REFERENCE: Dulacska, Lead. Engr., Planning Office, BUVATI of 
Budapest, HUNGARY), “ Buckling of Reinforced Concrete Shells,” Journal of the 
Structural Division, ASCE, Vol. 107, No. ST12, Proc. Paper 16753, December, 1981, 
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16752 SHEAR REINFORCEMENT FOR CONCRETE SLABS 
' 
Deformation; Ductility; Economics; Flat concrete plates; Flexural strength; th; 
"Reinforced concretes; Shear properties; Shear strength (concrete) 
| ABSTRACT: A review of the different types of shear reinforcement for concrete re flate 
{ plates is presented and the difficulties encountered when these are used in practice are 
eee Practical shear reinforcement should efficienty increase the shear strength 
and ductility of the slab-column connection, and should be easy to fabricate and install. 
= New Type of shear reinforcement which satisfies these requirements is presented. An 
individual reinforcing element is in the form of vertical bars with mechanical 
\ a _ anchorages at the two ends. The size and shape of the anchors are examined as well as 
_! the arrangement of the shear elements around the column. From the results of 40 full- 
“sie tests, rules for the design of the new type of shear reinforcement are given in a 
form of suggested revisions to the curren 
; x REFERENCE: Dilger, Walter H. (Prof. of Civ. Engrg., The Univ. of Calgary, Calgary, . 
Alberta, Canada, T2N 1N4), and Ghali, Amin, “Shear Reinforcement for Concrete 
Slabs,” Journal of the Structural Division, ASCE, Vol. 107, No. ST12, Proc. Paper 
16752, December, 1981, pp. 2403-24200 


16769 COMPOSITE DECK- REINFORCED SLABS et a wore 
KEY WORDS: Composite structures; Concrete (reinforced); Floors; Load 
limits; Shear strength (concrete); Slabs; Structural engineering BA 

ABSTRACT: The behavior of composite deck- reinforced floor systems under 
concentrated loads is described. Laboratory experiments show that these floor systems _ 4 


may have a considerably larger capacity to support concentrated loads than indicated 
by existing design techniques. There are two sources of this additional strength. First, 
continuity of the slab and slip restraint of stud shear connectors increase the shear 
bond capacity. Further benefit is obtained through the distribution of the load effects 
over the weak direction of the one way slab. The utilization of these benefits may 
_ produce significant economics in the design of floors to resist concentrated loads. bess a 
REFERENCE: Roeder, Charles W. (Assoc. Prof., Dept. of Civ. Engrg., Univ. | - = 
Washington, Seattle, Wash. 98195), “Point Loads on Composite Deck-Reinforced _ 
Slabs,” Journal of the Structural Division, Vol. 107, No. ST12, 
December, 1981, PP. 2421-2429 
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Customary-SI Conversion Factors 
In accordance with the October, 1970 outen ¢ of the ASCE Board of Direction, which ‘stated 
that all publications of the Society should list all measurements in both U.S. Customary and 
SI (International System) units, the following list contains conversion factors to enable readers 
i to compute the SI unit values of measurements. A complete guide to the SI system and its 
use has been published by the American Society for Testing and Materials. Copies of this 


; publication (ASTM E-380) can be purchased from ASCE at a price of $3.00 each; orders must a 


_ All authors of Journal papers are being asked to prepare their papers in this dual-unit format. 
r ; ‘To provide preliminary assistance to authors, the following list of conversion factors and guides — 


are recommended by the ASCE Committee on Metrication. 


yards (yd) F meters (m) 0.914 


square inches (: (sq in.) square millimeters (mm (m 
«Square feet (sq ft) square meters (m » 0.093 
square yards (sq yd) square meters (m 0.836 
square miles miles) square kilometers ters (km_* 
cubic inches (cu cubic millimeters (mm 
cubic feet (cu cubic meters (m 
cubic yards (cu yd) cubic meters 
(ib) mass kilograms (kg) 
pound force (lbf) newtons (N) 
pounds per square foot pascals (Pa) 
pounds per square inch iach kilopascals 
gallons (gal) fitters 
acre-feet (acre-ft) cubic meters 
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by a ‘single set of cables. The tall towers, which | are » used primarily < as radio a 
a television antennas, and for which wind and gravity are the primary loads, 7 
_ have been studied extensively and significant information about ee behavior = 
i is available in the literature (3,4,6). By comparison, short towers, which are 7 ; 
, used frequently today to support wind turbines and offshore platforms and 
are subjected to large concentrated and horizontal loads" 
ar 


- aim was to obtain the critical loads of axially compressed short towers supported : 
bya single set of guy wires. It is the purpose of this paper to extend of 
results of that investigation by describing the - elastic | post- -buckling. behavior 
of short towers. It is by now a well-known fact that the elastic collapse load — 
of an actual imperfect structure can be less than, equal to, or greater than — 
- the critical load, depending on the post-buckling characteristics of the system. _ 
4 As a consequence, a knowledge of of those characteristics is essential for ame 


‘The n nonlinearities. considered in the finite deflection analysis « described eit 
_ include axial shortening of the mast due to bending and axial stress, and changes 
in slope and length of the cables. The cables are, however, assumed to remain 
i Straight at all times, i.e., cable sag is neglected. Based on a more precise — 
investigation, which accounts for the actual shape of the cables, it has been 
concluded that cable sag is relatively unimportant for short towers (5). = 


ive. Using the theorem of stationary potential energy, the governing equations 


obtained in closed form and then solved by numerical analysis. 
Since the behavior of a guyed tower depends on the support conditions at 
q the base of the mast, the relative stiffness of mast and cables, the initial slope 
€ the cables, and the pretension in the cables, the effect of each of these 


‘Assoc. Prof. of Civ. Engrg. , Univ. of Massachusetts, Amherst, Mass. 01003 2 a 
Grad. Student in Civ. Engrg., Univ. of Massachusetts, Amherst, Mass. 01003 0 9 vy 
Note.—Discussion open until May 1, 1982. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, + 
ASCE. Manuscript was submitted for review for possible publication on December 22, J 
_ 1980. This paper is part of the Journal of the Structural Division, Proceedings of the J 
_ American Society of Civil Engineers, ©ASCE, Vol. 107, ters ST12, December, 1981. — 
ISSN 0044-8001 /81 /0012- 2313 / $01.00." Y 


ANALYSIS OF GUYED TOWERS 
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7 ~ 
factors is s investigated. In addition, with four cables, , 90° 


Derivation of Equations 
Consider the guyed tower depicted in Fig. (a). ‘The mast has length, iL 
_ of inertia, J; cross-sectional area, A; and modulus of elasticity, E. ms 
Each cable has an initial length, _S, and makes an initial angle, a, with the 


a _ horizontal. They have a cross-sectional area, A., and modulus of elasticity, — 
The ¢ consists of a vertical P, a horizontal load, H, applied 


4 
at the top of the tower. In addition, the cables may be steak to an nen 
- tension, T, which induces an axial compression, C, in the mast. Se eee, 
To obtain the differential equations and boundary conditions of the system, 


the principle of stationary sain — is utilized. The total potential energy 


in | wales U,, = the strain energy in the mast; U. = the strain energy in n the — 


= the potential of the external loads. Only the increase in energy 
jaan that existing due to cable pretension is considered. . 


BESS 

| 


ite 
in which u and w w = the displacements in aN z and x pe eer as indicated 
> in Fig. 1(5); and ’ y = differentiation with respect to z. In view of Eq. 2, the 
of gponing eqpations in 8. ant te 


ia which K = A.E,/S = stiffness of each ca cable; 4, = 
any cable; and m = the number of cables. 
ony, 
in which COS a; the plane containing any 
and the x-z plane as shown in Fig. l(c); and 8, and 5, = the horizontal and — 
vertical displacements of the top of the pole as indicated i in Fig. 1(6). ane, OF 
ge = PS, whee stffaes: tot © 
tt should be noted that Pi is positive when acting down. 
Combining the engucations in Eqs. 3 4, and $ oe the | total I potential ¢ energy 


II = P8, ul + (w" Mistanco die 


The chat this expression has: a stationary leads to the dif ff erential 


=0. 


and to the boundary fre Br dering acy 


] 
| 
| 
| 


7 


+ 


+w’ (w’) 

which ‘the notation [ Ve _ denotes that the is atz= 
‘a Eqs. 10 and 11 are conditions of equilibrium for forces in the vertical a 
lateral direction, seapectively, at the top of the mast. In addition to these boundary 
conditions, w, = u, = 0, and either w” = 0 for a hinged | base, or w,, = 0 
In the form presented above, Eqs. 8 and 9 are nonlinear and coupled in 
_ u and w. However, following a procedure outlined in Ref. 2, it is possible 
to the equations considerably. Integr: 8 gives 


= 


In view of Eq. 12, the axial is the entire length of the 
mast. Thus, Eq. 13 is valid not only at z = = L but for all values of z. ‘dente ni 
‘Elw" +|P+ +T)|w" =0; or w' = (4) 


L+ 
in which = — . 


Eq. 9 has thus been police from a nonlinear equation in u and w, toa a 


16 for u integrating th the expression. Thus 


ee 


| 
_ _ Making use of Eq. 13, it is also possible to remove u from Eq. 11. Combining 
bined by solving Eq. 


ety 

_ The final step consists of solving Eq. 14 and substituting the result into ool 
17 and 18. The solution of Eq. 14 fy) 


g w=A,sinpz+A, COs pz + A,z+. 

Making 1 use of the appropriate boundary conditions for either a hinged ora 
_ fixed base tower, the constants in Eq. 19 can be evaluated, and the resulting 2 


_ expression and its derivatives substituted into Eqs. 17 and 18. This leads to — =4 


Ely?) - — | 


two” sets of nonlinear algebraic equations in 5, and 8, ' whose solutions give 
the load-deflection characteristics of the towers. eninge © 
‘Fixed the base of the mast is fixed, the final equations 


+ tan 
tanpl—p 


_ it can bend like a hinged-hinged column, in - addition to rotating. The former 
of these two modes occurs only when the stiffness of the cable system is , 


_ small compared to that of the mast. Since this will rarely be the case for an — 
actual tower, only the bending mode will be ‘considered. ‘The equations governing 


4 Using the iteration numerical solutions of the 
nonlinear equations in 8, and 8,, derived in the previous section, have =e 
obtained. Although the natural procedure in dealing with equations of this type 

7 is to vary the applied load and solve for 8, and 8,, it is necessary, in order — 
to obtain convergence in the vicinity of the maximum load, to vary 6, and — 
solve for the corresponding load and 8,. If, during any iteration cycle, the - 

7 force in a cable becomes negative, indicating the cable to be slack, the term = | 

2 corresponding to that cable is dropped from the solution. All solutions have 

been out for at tower 100 ft 5 m) hi high, a 6. in. 3 


| 
nT sin a « > 
: 


diam pipe mast. ‘Both the he cables are steel with E = ; 
ksi (200 x 10° kN/m’), and E -= =o x 10° ksi (165 x 10° kN/m a 
been assumed that the plane containing H, if a horizontal load is present, ond - # 
the plane in which the mast bends with or without H, coincides with the x-z _ 
i, as defined in Fig. - This assumption results in the most critical atin 


é Fixed Base.—Fig. 2 depicts typical post-buckling curves for fixed-base towers. 


The ratio of the applied load, P, to the critical load of a fixed- -hinged column, 
Be is plotted against the ratio of the lateral displacement, 5,, to the tower 
- height, L. The ‘cables that support the mast are assumed to | have a a small amount 


1 


2 2 -—Post- Post-Buckling Curv Curves: for Fixed- -Base 


: However, compared to the applied load, = the tension in the cables at the ra 


onset of buckling is assumed to be negligible. 2: 
> It is evident from the curves that buckling commences at a load very dese 
_ tothe critical load of a fixed-hinged column. This result agrees with the conclusions 
obtained using a linear analysis (1). In that analysis it was shown that only 
a minimum cable stiffness is required to produce a lateral restraint at the top 
- of the mast that ‘is nearly equivalent to a hinged support. The curves also 
indicate that the vertical load, P, which the tower can support, decreases gradually 
with increasing lateral deflection subsequent to the onset of buckling, and that 
“this decrease in load becomes more pronounced as the initial slope of the cables, 


q@, increases. M beck ods ovine 


a As the mast | bends, an increasing tension force develops in the cables ony 


a side opposite to the deflection. As a consequence, the mast is subjected a 4 
toa aecmenies force equal to the vertical components of these cable tension | 


Bretension, su eit tO Prevent te Mast OUCKIINP as a 
8 
_ 


— 


_ forces, and the applied load, P, that the tower can n support is accordingly reduced. _ 

- ‘The larger the initial slope of the cables, the larger the vertical components | 
of the cable tension forces, and the larger the reduction in P with increasing — 
lateral deflection. Except for towers with extremely flexible cables, which are 

1 rarely used, changes in the diameter of the cables or in the number of cables — s 4 

a do not t significantly a affect | ‘the shape of the post- buckling curve of the oe J 


= characteristics of the tower. It indicates that the presence of cable pestedaion : 
_ reduces the critical load by an amount very nearly equal to the vertical components _ 
of the cable tension forces. In other words, for the system described in Fig. — 

_ 3, the vertical components of the pretension forces add up to approximately 

0. 3P, and the critical load is thus equal to 0. 7P, 
Comperioen of curve in Fig. 3, with hess in Fig. 2, also shows that 
the presence of cable pretension markedly affects the shape of the post- buckling 


3.—Effect of Cable Pretension on Post- -Buckling Curves of Fixed- Base Towers 
_ curve. Instead of a gradual decrease in fi with increasing deflection, as was _ 
the case when no pretension existed, the curves now exhibit a sharp increase bs 
: ‘in P as the mast begins to buckle, followed by a gradual decrease in P as 
7 - deflection continues to increase. The system thus exhibits considerable 
_ post-buckling strength initially before it eventually becomes unstable. =| ‘a 
The difference in | behavior between towers that possess initial cable tension 
= those that do not can be explained by considering the relative values of : 
8, and 8, for the two cases in Table 1. For small deformations, the tower 
with initial cable tension, 7, undergoes an axial shortening, 5,, that is three | 
_ times as large as the lateral displacement, 8,. As a consequence, the cables 


cable tension, the lateral displacement, 8,, exceeds the axial shortening, 5 


of the mast, causing a buildup of tension in the cables sociating aad deformation. 


z? 


a increase in P. By comparison, in the tower no 


= 
| 


This —_— the mast to an increasing compression, G, ond results in a reduction . 
_of the applied load, P, that the tower can support. > > 
_ Fig. 4 depicts the effect of a horizontal load on the ability of a tower to 
a vertical load. The horizontal load is assumed to be applied first, and y 
-.. vertical load thereafter. The solid curves in the figure represent the behavior 
of three- and four-cable systems with H = 0.2P,. For comparison, systems ~ 
_ with zero horizontal load are also inten. ‘They ‘ are indicated by the dashed 
_ For each of the two cases sssauiinid di presence of horizontal load 
; significantly reduces the ability of the tower to resist vertical load. The four-cable 


ABLE 1.—Displacement at Topof Tower 


d ne. 4. ‘.—Effect of Horizontal | Load on | Carrying | Capacity of Fixed- -Base Towers 
system with a horizontal load is able to support a maximum load of 0.67P,, 
_ and the three-cable system has its carrying capacity reduced to 0.51 P,. A horizontal 7 


load, like a lateral deflection, causes tension forces to be developed in the = 


aft cables. These forces, in turn, induce a compression force in the mast whose 
he reduces the ability of the mast to support an externalload, P. 

It is evident that a greater reduction in the capacity to resist vertical load ~ 
Ss for the three-cable system than for the four- cable system. This is rom 
to the fact that a given horizontal load induces a . vertical compression in the 


- equal to 2H tan a for three cables, and V 2 H tan a for four cables. ‘aa 
An increase ee othe slope always increases the compression induced in the 


| 
| 
| 
=) 
| 
q 


Thus, the in the carrying capacity of 
_ the tower due to the presence of a horizontal force will increase as a increases. 

_ As indicated in the figure, the difference between the critical load of a tower 
without horizontal load, and the maximum carrying capacity with horizontal 
load, can be subdivided into two parts labeled r, and r,. The first of these 
is the reduction in carrying capacity that results from the unstable post-buckling 

_ characteristics of the system. In other words, at the lateral deflection caused 
_ by H, the tower can resist less vertical load than it can at zero lateral displacement. — - 
_ The second reduction, r,, is due to the compression induced in the mast by | 
the cable tension forces required to resist H. While it is possible to obtain “ 
a very good approximation of r, by assuming the mast and the aft cables to 

4 behave like a truss, the value of r, can be obtained only from a nonlinear > 

yy In the foregoing analysis of the effects of horizontal forces, it was assumed _ 

‘that the cables had zero pretension. If large cable pretensions are present they 
do not ot materially alter the effect of H. A horizontal force will still cause a 


FIXED BASE HINGED ‘BASE 
6. 5.— —Deformation ¢ of Guyed sd jon! 
‘reduction i in the maximum carrying capacity « of the tower becanes still i induces 
a sizable compression os ed | 
Hinged Base.— -Except for the fact that the critical load of hinged-base towers 7 
is approximately half as large as that of fixed-base towers, the behavior of . 
_ the two is very similar. The most significant difference between the two is 
that the decrease in P, with increasing lateral deflection, is much smaller for — 
hinged-base towers than for fixed-base towers. In other words, the post-buckling i. 
curve of a hinged- tower is almost horizontal. This difference in 
is largely due to the different ways in which the two systems deform. As shown Jf 
in Fig. 5, the fixed-base tower has a larger lateral deflection at its upper end 
than the hinged-base tower. As a consequence, buckling induces larger cable 
tensions and more mast compression in the fixed-base tower than in the 7 
hinged-base tower, and the former has a more unstable post- eating curve > 
a Numerical Example. .—Consider a a fixed base tower consisting of a 6.625 in. 
aed (168.3 mm) diam pipe, and three 0.25 in. (6.4 mm) diam cables. The mast and 
, cables have the following properties: L = 100 ft (30.5 m); J = 28.1 in. *(, 690 


mm‘); A = 5.58 in.’ mm x ksi (200 


| 
i 
4 


= “049 | in. ? 3.2 mm’); E. =e 10° ksi (165 x 10° KN/m’); a = 45 
Ref. 1 indicates that only a minimum cable stiffness is necessary to produce 
a lateral restraint at the top of the mast equivalent to a hinged support. Thus, 
3 the critical load of the given tower is close to that of a hinged-fixed column, 
ie., P, = 0. 99P, = 11.7 kips (52 KN). Since the cables are already producing» 
— the maximum lateral restraint nt ‘possible, the critical load of the tower can | be 


diameter from 6.625 in. (16.83 cm) to 8.625 in. (21. 91 cm) increases the oritical 

load from 11.7 kips (52 KN) to 30.2 kips(134.4KN), 
- According to Fig. 3, an initial cable tension of T = 1.2 K (5.3 kN) reduces 
r - the critical load from 11.7 kips (52 kN) to 8.2 kips (36.4 kN). However, as _ 
_ soon as the mast begins to bend, the carrying capacity of the tower increases 
| to 11.2 kips (50 kN). T hus, even a relatively large cable pretension does not 
significantly decrease the carrying capacity of the tower, tt” tedt 
- Finally, Fig. 4 shows that a horizontal load of H = 2.4 kips (10.7 KN) reduces 
= the vertical load that the tower can support by 50%. By comparison, if a is 
i increased from 45°-75°, it takes only a horizontal load of H = 0.6 6 kips (2. 7. 


: kN) to reduce the ability of the tower to resist vertical loads by 50%. 


g. Except for systems in which the ratio of cable stiffness to mast stiffness 

. is extremely small, the critical loads of hinged and fixed-base towers with small - 

_ cable pretension forces are very - close to the critical loads of hinged-hinged __ 
and hinged-fixed columns, r respectively. 
2. Towers with small cable pretension forces have an unstable post- -buckling © 

_ curve. However, except for very steep guy wires, the decrease in vertical load ‘ 
with increasing deflection is relatively small. As a consequence, these systems — 

_ will not be “‘imperfection-sensitive’’ and the actual buckling load of a real tower 
can be expected to be fairly close to the theoretically obtained critical load. 
_ 3. The presence of sizable cable pretension forces decreases the load at a 

; __ buckling begins by an amount equal to the sum of the vertical components 
of the cable forces. However, as soon as the tower bends, the cable tension a 
forces decrease rapidly with a corresponding increase in the capacity of the 
tower to support vertical load. The net effect is that a tower with large cable 
= forces can support a - vertical load that is only slightly smaller than — 

7 the load a tower without sizable cable pretension forces can support. EES 


: : 4. The presence of a horizontal load results in a significant decrease in the _ 
: capacity of the tower to resist vertical load. Since this decrease is primarily - 
caused by the compression that the horizontal load induces i in the mast, increases — 


7 
of the tower to resist vertical load. 
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A = cross-section area of meat 
A. = cross-section area of cable; this pws 


CC = axial compression induced in mast by cable NE 
» ‘= horizontal projection of cable length; 
= modulus of elasticity of mast; 
} = lateral load applied at top mast; 


vertical load applied at top of mast; aw 
tag P, critical load of fixed- hinged c column; a 
= reduction in load carrying capacity due to load; 
U = strain energy in cables; _ the 
= strain energy in mast; 
displacement in z-direction; 
displacement in x-direction; 
coordinate axis normal to mast; _ tes 
= coordinate axis along mast; 
= initial angle between cable and horizontal plane; 
B, = angle between ith cable and x-z plane; 
; 


cable; and 
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NONLINEAR DYNAMIC ANALYSIS ‘BY 


“ structural response to dynamic loads. In some instances this allows og 
_Tefinement of our understanding of the behavior of a particular system, while — 
in others, it is is almost imperative if if we are to have even a general appreciation — 
of the response of a particular system. An example of the first case might 
involve the inclusion of partial rigidity of joints in a steel frame, while the | 
“second case might relate to the geometric nonlinearities associated with cable 


_ Numerous methods for dealing with the nonlinear dynamic analysis problem 


proposed by Newmark | (7), Houbolt 6), and Park (8). Explicit methods such 


& been presented in the literature. Some implicit approaches have been pay 


as central difference predictors and the Runge-Kutta approach are also used. 
Each has its advantages and disadvantages. A good review of some of these — 
_ The method presented here uses the modal superposition approach, corrected - 
to account for the nonlinear "aspects of the system response. The method is 
_ presented by examining the results for cable network problems, where the Pare _ 
in geometry under load provide the rm of the system. 
Motion 
‘The: equations of motion for a nonlinear ne be written, in matrix 


+ + + +(R@)) = 
which = diagonal mass matrix; [K] = linear stiffness matrix; 
= nonlinear terms, as a function | of displacements; (P(t)} = dynamic loads, a 

as a function of time; and 8, $= = displacements and accelerations, ‘Tespectively. 
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-" t obtained as a combination of the independent patterns described a ‘the 


mode shapes and corresponding modal amplitudes. They may | be written as 


= 
which = the matrix of mode | shapes; and Yo and = 
_ emplitedes. ‘Substitution of Eqs. 2 and 3 into Eq. 1, and eo to 
a 


if the right hand side of Eq. 4 is defined as — 


4 


ast 
+ [K] [>] (Y@) = = {Q(4,5)} magne 40 (6) 


Since mode shapes from the linear analysis will be orthogonal, these equations 


may be uncoupled if the system is premultiplied by the transpose of the nth FS 
IM] + [0] (YO) = COGS) 
Due to the the orthogonality of of re shapes, all 1 but the nth nth terms | vanish, which — 


¥,(t) + (X)7[K1(X), Y,@) = (X}7(00, 


Ea. 8 may y be written as 


= 


and (18) = (29 


The: terms. M, and K, = the gen generalized mass and ‘stiffness. They a are re the 
- same for this nonlinear formulation as they would be for the linear a an 
The generalized load, Q,(t,8), however, is different and is the term through 
the nonlinearity is is modeled. Substitution of Eq. 5 inte into Eq. 12 yields 
4 ‘If the nonlinear contribution were known ahead of a the generalized load, 
S 6s €@, (t, 8), could be determined, and Eq. 9 could be solved as a typical linear a 
gingle-degree-of-freedom problem using any acceptable method. The nonlinear 
contribution, {R(8)}, is generally not known a priori and may be different at “ 
any time, ¢. Thus, it can be said that {R(8)} is, actually, indirectly a function 
of time. In order to factor the time function out of Eq. 13, the nonlinear 
. 2 contribution is treated as a known load. Eq. 13 is then multiplied and ue 


by {X}7(P,}, in which | {P, 7° the vector of load a. Thus, Eq. 13 
becomes 


4 

| = 
+ K,, Y(t) = (0,8)... 
ayn 
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The ‘brackete 


d term in Eq. 14 represents the variation i in a tiene, fe, e (t)s of re 


effective load, with the other term representing the generalized load, 


. The s solution « of Eq. 9, as ‘as given by the Duhamel nel integral, i is given | by 


. 15 into Eq. 16, and expansion of the yields 


2d) * 


\ fe (7) sin w, tdt 


in which (PF), = the participation factor for mode m and remains constant 
for any time, ¢; and (IAF), = the Instantaneous Amplification Factor. This 

IAF * will differ from the linear IAF due to the inclusion of the {R(5)} terms 
through the time variation of the effective load. 


_ After the natural frequencies and mode shapes are determined for a linear s 
formulation, the nonlinear modal s ‘superposition scheme may be employed a 
determine the nonlinear forced vibration response. As can be seen from Eq. y 
17, the participation factors remain constant throughout the entire analysis and 
the same as for a linear analysis. 
Since the nonlinear terms are a function of displacement, and the displacements — 
are small at the beginning of the respons., an assumption that {R(8)} = 0. 0.0 
for the first time station is generally a good one i} ® a 
sg The evaluation of the IAFs, the bracketed terms in Eq. 17, for a general 
time function, requires a numerical integration scheme if a closed solution for 
the integrals is not available. This holds for either the linear or nonlinear 
formulation. With the nonlinear problem, the effective time function, 


(RO). 2 ‘Assuming a parabolic variation of the. effective time function, the integrals 
in Eq. 17 are evaluated and the IAFs for this trial are obtained. The modal 7 
amplitudes, Y,(t) and Y,(t), are then obtained in the same manner as for 


a linear analysis. The corresponding displacements and accelerations and 
7 5, are calculated ‘Eqs. and 


1 
— 
| | 
| 


with the assumed values. For the first time station, the ‘nonlinear terms were 
taken initially as zero. The calculated terms are small and the difference between 
the assumed and calculated values is usually less than the set tolerance. If — 
that is the case, the calculated displacements are accepted as correct and the 
calculated nonlinear terms are used t to establish the continuous effective time 
4 For the next time station, an on Seienigiiais is needed for the nonlinear terms. 
— It is possible to assume these terms to be zero; however, assuming them to 
be the same as those at the end of the calculations for the preceding time 
_ Station will somewhat speed convergence. With these newly assumed erm, 
> “the process is repeated with trial IAFs and the resulting displacements and 
accelerations. When the difference between the assumed and calculated nonlinear — 
terms is less than the permitted tolerance, the displacements for this —_ ne are s 4 
Essentially, this is the solution of a linear problem for each ‘iteration, — 
the nonlinear terms applied as part of the time dependent loading. The modal 2 
method fora a linear problem allows for the evaluation of system 


nonlinear terms, are then ovelonted from Eq. 


in cases ses where the integrals of Eq. 16 « can n be evaluated. For the general loading 
§ function, the integrals must be numerically evaluated. Thus, there will be some 
5 increment of time used to evaluate the IAFs. For the nonlinear problem, the > 
- effective loading function is not known prior to the evaluation of the nonlinear os 


_ terms. This means there is a greater concern over the the time increments for 


= 


_ +The structure studied is the Aden Airways Building, Fig. 1, presented by a 
- _ Knudson (6), and previously used by the writer for parameter studies and forced 
vibration analysis (2,3). The static equilibrium configuration is given t 
in which a = b = 57.5 ft (17.5 m); C, = ong in x direction aii 5 ft 3.81 — ; 
- 4 ro = rise in z direction = 17.5 ft (5.33 m); - and f = elevation at «= 
z = 0) = 12.5 ft (3.81 m). 
The additional structural data needed for an are as follows: 
Size of network—115 ft x 115 ft (35.1 m x 35.1 m); 4 
Number of cables each way—7; 
Nominal diameter—I- 13/16 i in. n. (46. 0 
Mass ‘peru unit length—0. 209 Ib-sec’ / ft? (9.98 kg/m); 


= 
| 


wor un 


Breaking strength—404 kips 800 KN); [at 


Modulus of elasticity — —24,000 ksi (16,600 MN /m 


Horizontal component of initial cable tension: 


(longitudinal)—45.0 kips (200 kN); 
(transverse)—32.14 kips (143 kN); 


iven in Table” A complete of the modes and frequencies for this 


g 
network has been presented by the writer and West (2). They have also presented — 


forced vibration for this ‘Newmark’s B method (3). 1 The 


TABLE 1.—Natural Periods | for Aden ‘Airways Network amy 


(2) 3) (4) (5) 


0.022 #| Si 


oi 
0. 7 54 


_ formulation of the various components of Eq. | follows the development presented ~ 
_ Time Increment.—Of particular interest in the development of any dynamic : 
analysis is the required time increment for system stability and convergence. 
_ Knudson (6) and Goudreau (4) have presented studies relating to the selection 
Of the time ¢ increment and the other factors associated w with Newmark’s B method. = 


| 
4 Number 
0.006 
8 | o19 =| 33 | oon | ss | 0.005 
0.087 42 | (0.008 = 0.005 
0.062 47 | 0.006 72 0.004 
OF a Value 5 = the Well KNOWN linear aCceicration Method, 
method, is unstable for a time increment, At, greater than 0.5517, in which 


pares 


= the pe period with | the largest of 
at the Aden Airways network is subjected to a 2 kip transverse load (+z 
- direction i in Fig. 1), with a duration of 0.08 sec as shown in Fig. 2, the response - 

- has been shown to be very nonlinear (3). In order to investigate the influence _ 
“of the time interval, this system was studied with the Newmark linear and — 
nonlinear 1 methods, and the proposed method. The vertical displacements at 
joints | and 25, at a time ¢ = 0.06 sec, are the same e for both nonlinear approaches. _ 
The difference is seen in the number of iterations required to complete the — 
analysis. “tah shows the relation between the number of iterations required a 


z(trang) 


4 = 


= 


FIG. 1.—Aden Airways Network Geometry (1.0 ft = 0.305 m) 
requires the same number of iterations fc or either the linear c or nonlinear analysis. 

It is also seen that for a time interval greater than the previously discussed 
stability limit, this approach is unstable. The proposed method requires only 

two cycles, even for a time interval of approximately three times the Newmark © 

_ Time Dependent Displacements. .—The dynamic response of the Aden Airways — 
network to loadings with the time variation of Fig. 2 has been presented in » 


the literature (3). For cases where the loading is antisymmetric, a significant 
nonlinear response is seen. Figs. 4-6 give the response of selected joints of 


= the Aden Airways network subjected to three different patiorns. 


- 
| 
“| 
Ny 
z=. 


a is ee: for a vertical load of 2 kips per per joint applied in - 
fond oF? ke pattern about the main longitudinal cable, case I; a transverse ; 
load of 2 kips per joint in the z direction, case II; and, finally, a combination — 
of the 2 kip transverse loading, with the 2 kip vertical load, case III. _ The 
function is an impulse with the duration of 0.08 sec. __ 


Ample Versus Time fe for Load 


& NONLINEAR 


@ proposen NONLINEAR 


FIG. 3.—Iterations Versus Time Increment for Proposed Method and Newmark’s 


7 For the loading of case I, , joints 10-16 are not loaded, while an upward force , 
of 2 kips is applied to joints 1-9, and a downward force of 2 kips is applied © 
- to i 17- 25. The dominant mode is mode 14 with a period T = 0.100 sec. — 
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the corresponding mode | numbers, a are e given in in Table 
_ The time dependent vertical displacements associated with joints | and 25, 
the joints which exhibit the greatest degree of nonlinearity, are presented in 


Fig. 4. In all cases, the analysis has been carried out using all modes. ‘The . 
figure shows the displacements for the nonlinear ‘Tesponse, and the displacements _ 


JOINT 25LINEAR 


JOINT 25 LINEAR 


“eid 


JOINT 25 NONLINEAR 


ERTICAL DISPLACEMENT (ft.) 
co} 


from the linear modal superposition said The difference in stiffness between 
; q an upward displaced joint and a Cowawaed baal joint results i ina significant 4 


_ The upward force applied at joint | is resisted by a longitudinal cable which 
_ it tends to put into reduced tension. Thus, the settee are greater for — 


| 
| 
| 
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Ge joint 1 in the nonlinear analysis than they are in the linear analysis. The onpenie 
is true for the downward load on n joint 


= hg the nonlinear seepense, is also evident in Fig. 5. As for « case I, the linear ‘and 
y nonlinear response is given. For the transverse loading, mode bs is again the 


_ The loadings from the previous two cases are combined so that their meen 


be —Participation Factors for Loading Case 


Mode number Participation factor 


0.0559 


JOINT 25 NONLINEAR 


ODS 


ISPLACEMENT (ft.) 


VERTICAL 


FIG. 6. —Vertical Displacement Versus Time, Case Ill (1.0 ft = = 0. 305 m) Prd: D 
will reinforce each other. The results are presented i in Fig. a Again, the linear 


e - nonlinear vertical displacements are ‘presented. The amount by ' which the 
_ nonlinear solution varies f rom the linear solution i is even en greater is in ie combined : 


response of cable networks is adequate; however, the results presented oceal 
and previously (3), indicate that there are loading cases where consideration | 
of the nonlinear response is necessary for an accurate analysis. 


A direct cmgeeaen of the displacements obtained fr rom the linear r acceleration 


| | 
4: 
| 
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_ Time Variation of Forcing Functions. —The variation of load magnitude, with __ 
time for four simplified time functions, yields the dynamic response indicated : 
ie: Fig. 7. The loading pattern of case I, with three triangular pulses and the 


og 


of Fig. 7 which are the vertical displacements for joint 1. A comparison of i 

“the four curves indicates a response to the time variation as might be expected he 

fora linear system, the are modified to reflect t the nonlinearity 


w al 
m 
) 
; ; FIG. 7.—Vertical Displacement Versus Time for Four Variations of Forcing Function __ ; 
(1.0 ft = 0.305 j 
8.—Vertical Displacement Versus Time for Four Load Magnigudegett 
| 
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a few of the cables have a “very low | tension. 
a the load magnitude is Pama to 4, 6 and 8 kips, and the elements are 
‘felis to take compression, the resulting response is as indicated in Fig. 
_ 8. The time increment for these results was taken as 0.00333 ‘sec. The analysis 


two iterations ateach time station. || 
The lit linear ar response, as predicted by tl the modal superposition approach, ; simply 


can n be seen that the nonlinear response reflects the changes over time of the y 
= system stiffness. Thus, the period of the response is not the same for all load — : 


-a new method for nonlinear dynamic analysis. Using this modified modal 
superposition approach, a number of response studies are carried out pie 
_ simplified forcing functions and various loading patterns. be 
‘The proposed method has been shown to successfully determine the saaiiaiis 
dynamic response for systems which exhibit significant geometric nonlinearities. 
_It accomplishes the analysis with much larger time increments and significantly 
fewer iterations than Newmark’ s B method, using the linear acceleration formula- 


This study has presented a method for handling the; response 
of systems with geometric nonlinearities. Its applicability to problems with 
material nonlinearities needs to be treated. The method converged well — 4 
- remained stable for the examples presented. Further study i is required to ms geo - 
theoretical limits for stability and i 


The required these studies were out using the 
370-3303 system of The Pennsylvania State University Computation Center, 
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following symbols are used i in this paper: 


one-half of longitudinal span; 


C= = sag in x direction; watt wie 
ff = elevation at x = z= =0; ony 
= time variation of effective load; 


= instantaneous amplification factor; 


« stiffness; 
= mass matrix; 
M, = generalized mass; 


P. F. = participation factor; 


= vector of load amplitudes; 


vector of terms nonlinear in 


= 
ewe 

= mode shape vector; 


longitudinal coordinate; 
modal displacement amplitude; 
“modal acceleration amplitude; 
vertical coordinate; 


transverse coordinate; 
acceleration; and ne a 


ody 
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OPTIMIZATION W ITH FREQUENCY | ConstRAINT 
William R. R. Spitters," M. ASCE, Sohan Singh,” and Levy” 


- This paper is concerned with the most simple structural optimization problem 
with frequency constraints: Find a minimum weight truss for which the funda- 

_ mental frequency is greater than some given value, given the structural configura- _ 
tion and the masses. This is, of course, not a new problem and the ee 
here is not with the development of an algorithm (although certainly there are 

_ implications in this area). The primary concern lies with attempting t to explain w 


why ‘ ‘allowable s stress type”’ algorithms. are successful when applied to problems _ 
_ There will be no attempt made here to review comprehensively the literature — 
on optimization with frequency constraints. The interested reader could well 
begin with the recent paper of Cassis and Schmit (3), or the recent review _ 
of f Venkayya (8), but will | not, in any case, find it difficult to get into the | 
literature. Neither is there concern here with finding a workable algorithm, 
since several different algorithms have been used to deal with this type of t 
pore tars Fox and Kapoor (4), for example, point out that the derivatives of 
eigenvectors and eigenvalues with respect to member stiffness were available 
in 1962, making any type of Newtonian algorithm feasible at that time, and 
Young and Christiansen (9) discuss an early aerospace application. In fact, 
since the problem under discussion is a basic nonlinear programming problem—a _ 
linear objective function with a single equality constraint—it is a candidate 
for solution by many of the classical procedures. 


__ The concern here lies with ‘ “allowable stress” algorithms in in rn —_ 


applied to problems of frequency constraints. “Allowable stress” algorithms 
have proven to be generally effective in dealing with a wide class of structural — 
4 optimization problems. Their principal feature is that they tend to reduce the . 
al computations required for optimization to the repeated applications of methods 

of structural analysis, with repeated | simple intermediate redesign steps. The ; 
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_. Venkayya and his colleagues at Wright-Patterson Air Force Base have been 
~ some of the most active proponents of “‘allowable stress’’ methods for optimiza- 
en, ond the scope of their applications (8) is quite remarkable. While this 
_ entire approach must on the whole be regarded as effective, ‘‘allowable — 
> algorithms can have some theoretical difficulties when applied to situations of 
multiple loading conditions and, naturally enough, don’t seem applicable to ne 
_ There has been some concern over “‘allowable stress’’ algorithms when it } 
comes to motivation. They do not seem to relate directly to any of the standard 


_ mathematical programming algorithms and, for the most part, very little has : 


been shown formally concerning their convergence. This, of course, takes nothing 
: away from them since, in any actual application, the engineer can see whether 


or not the weight is being reduced. 


_ This paper then attempts to determine why ‘“‘allowable stress’’ algorithms 
work when they are applied to problems with frequency constraints. It extends 
earlier efforts (5-7) concerned with placing “allowable stress’ algorithms within 


‘the realm of classical mathematical programming theory, and is built around 
‘the problem of minimizing the structural weight, given the laa, when the 
fundamental frequency of a structure is specified. 
_ Following this introduction a formal problem statement is developed, the 


cigenvalue problem i is discussed, Properties of the algorithm (such as monotone 


“two examples ar are e presented. It is s hoped that | through this ‘simple optimization oo ; 
problem it will be possible to see more clearly the mechanisms at work and, 
thus, develop some intuition concerning both ‘allowable stress’” algc algorithms and 

the structures to which they are applied. 


Prosiem STATEMENT 


matrix, given the matrix, so that the fundamental (lowest) frequency of the 
wr 


Structure remains above some specified value. This will b be written as 


q. to ‘the truss problem in which vector of the 


Ss of the member allowable length changes, i.e., (A2), = (o,L,/E)* in 
which o,, L,, and E = the allowable stress, member length, and Young’ s modulus, 
‘respectively; K = vector of member stiffnesses, i.e., (K), = / A, E/L, in Seu 

= the area of member i; = ,(A2), (K), = (02 /E)%,A,L,is proportional 
to the sum of the volumes. of the individual members; and a” = square _ 
the fundamental frequency which is to be larger than the given value ¢. —_ 


in which K, = swrecteral stiffness matrix; M = mass matrix (which j is ora 
7 and 6 = eigenvector. Eq. 2 is called a generalized eigenvalue ees to distin- 


| 
i 
| 
q 
; fa iS paper 1s concerned with the problem of selecting an optimal stiffness s - 
i 


guish it the standard problem Ax = 


FREQUENCY RESTRANET OPTIMIZATION 
_ In any case, to solve Eq. | it is necessary to find the member stiffucsses 


& which minimize the structural volume and satisfy a asingle frequency constraint. 
a is assumed that the « changes in structural volume do not significantly affect 
the mass matrix, M, which is taken to be fixed during the optimization. This - 
is certainly the case for structures, such as buildings or bridges, when the — 


With regard to other types of structures, wr is clearly an identical plane — 
Other formulations tend to be 


(stress, strain) finite element analog to Eq. 
For a general discussion of the structural eigenvalue problem, ec reader 


mass is primarily due to nonstructural elements such as floor ee ee 


can consult the recent paper of Bathe and Wilson (1). The ones of this 

section is simply to collect some results which will be required later, 


if only the fundamental frequency is of interest, t, Eq. 2 is S equivalent to o the 


but the inequality i is valid for any 8. 
Several useful features of the eigenvalue alli now can 1 be listed: 
> Scaling—If the structural stiffness matrix is ‘multiplied by a scalar a = 
0, so are the frequencies, K, a K, > w’ a This follows directly 
 ° Concavity—Beckenbach and Bellman (2) give an elegant proof of the old — 
result Courant that a of the matrix K,, orthat 


(aK! 
n which K’ and K” = two different structural stiffness matrices. (This result 
, the be conven- 


3. Inverse Iteration—Given « w? 


wen 
ed 


This is analyzed i in in Ref. 1 
Functions—Let a be a scalar. An Any scalar hous 


‘te 


f(ax) = 
is one to be homogeneous: of degree n. It follows that 


1 
7 
| 
| 
| 
&g 
_ 4 in which H = Hessian matrix of the function f. Since the structural stiffness — 


‘ vanadate the vector Vw ? then must be homogeneous of degree zero. ek 

§. Derivatives of w?—The derivative of w?, with respect to any member 

__ stiffness, (K),, can be obtained in a relatively simple manner and is an identifiable 

term. Let the 0/0(K), be indicated a comma. 


= 02 (K,, ,M)3 + (K, -o™M *M)8, 


asa 


q. 10 can be simplified further. First: ane a, it meen convenient | to normalize 


= 


FIG. 1.—Concave Funct 


it is convenient to matrix K, as 


in which each term in the sum represents the contribution of | 


ceding property 
a 


‘the matrix K (see Ref. 6 for the de detailed | of this 


follows directly that" Aside & of ‘ 


N,& = the member displacement (length change for the truss). , od 
We Eq. 14 is a simple, but important, result. It identifies the derivative of wo” 
a with respect to the member stiffness, (K),, as equal to the square of the member 
_ displacement associated with the normalized eigenvector,8. 8” 
6, Concavity and Homogeneity—Let f(x) bea concave, homogeneous function | 
of degree 1 so that f = . FF -m. Furthermore, let x, and x, be two points on 
the surface f(x,) = = f(x.) = ¢ (Fig. 1), and 71 to Vf, tt follows ‘that 


One; possible approach to the solution of Eq. | is through the use — Lagrange : 
multipliers which reduce constrained optimization problems to dati 
ber this method, the Lagrangian is firstformed as 


4 


‘Note: It is only necessary to consider the case of equality constraints in Eq. ; 

1, since any solution with w’? > c can be reduced in weight t by ‘scaling to 
the surface w = c. Let or tefer to the objective fu function 

and the asterisk denote an optimal pc point 


The optimal Lagrange multiplier, 

b 
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4 
| 
| 
| n be computed by multiplying Eq. 19 
(23) 


v 


Eqs. 19 and 20, the on problem of of Eq. at has been reduced 

to a nonlinear system which requires that Vo* be parallel to A2 - This fact 
g the basic motivation for treating Eq. | as an ‘‘allowable stress’’ problem. _ 
Since Eq. 14 shows that ‘Vor ‘is proportional to the member displacement, 

and since the eigenvector, 8, can be assigned any length, 5 could be sealed 

to give \ = 1, which would imply all members to be fully stressed for an 7 4 
- optimal dete. That scaling will not be used here since it tends to mask the ©} 
basic features of the eigenvalue problem. In any case, the concept of a “fully coe 

stressed’’ algorithm remains useful and will be pursued in the next section. ot 

_ From another point of view there is little reason to think of Eq. | as an 

allowable stress problem. Beyond Eq. 19, allowable stress problems are consider- 

ably different from Eq. 1 since they concerned with nonhomogeneous 
equilibrium equations rather than homogeneous eigenvalue p problems. This simply 

, makes a study of “‘allowable stress’’ algorithms more interesting since they 

ovest then have a broader scope to ney. problems which are quite different. 


now can be shown that Eq. | is closely related to other problems with 


Problem |.—Given a solution to Eq. 1, the solution for —_— other frequency, — 


_¢’, can be obtained simply by scaling K — K’ = Kc’ /¢. (This f follows ‘directly: 
from the fact that Vw? is homogeneous of degree 0. Don: 
_ Problem 2.—Givena solution of Eq. |, the solution of the problem of maximizing 

the frequency, given the volume, follows tty, by ane scaling K, to the 


Let K be arbitrary, but such that it satisfies w?(K) = c, and K, solve Eq. a 
1. Then it is possible to select an a < 1 such 
in which w? wf = ¢ is concave. Now scale point K down to aK which has the 
same ¢ weight as K*. Since is homogeneous of degree 2 follows that 


25) 


‘)= = (K) = w*(aK) . . 
it, _As indicated in Fig. 2, to solve Eq. l, it is necessary to find a a point ‘K, 7 
on the surface w* = ¢ such that A?. The “allowable stress’” ‘scheme 
to be discussed here iterates in the following manner, Starting with an arbitrary i 


‘ 
| 
Step 1. Given K,,,, solve the for « ow and « w? , and scale 


Step 2. 2. Compute | the new So that 


| 
| 
| 


is approximately u until the | is solved, it is 
compute. replacing Vw, in its approximate value: 


= 3 
—= are several properties of this algorithm which can now be demonstrated: 
_ 1. The optimal point is a fixed point of the algorithm, i.e., K,, so 


~ 


K.,)- This follows directly from Eqs. 27-28 in a somewhat trivial manner, 
_ Property = 6,,,/,. This follows from the fact that 


_ Property 3. In the vicinity of ga ae > c. This property | uses ses the linear 

‘terms ina a Taylor series expansion and is, thus, a local theorem. . Since w? 


toe 


vanes es of Constant 
, Weight 
a - 


A 


FIG. 2 —Surtace of Constant Frequency 
isa eee function of degree 1 in K, it follows, from a Taylor series 


| 
a 
| 
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jr 


(using property 4) 


vector equation by iterating, following Eq. 28 which 
square roots. It is thus Schwarzian i in the sense that the square roots are compatible ‘ 
with Schwarz’s inequality. Furthermore, it possesses the peculiar feature of — 
equating the scalar product of two vectors term This is is unusual 
but can be motivated somewhat by the 


— 


(40) 
in which K = an arbitrary point on the surface w? = c; and K is: again in optimal. 7 


In some sense, Eq. 31 involves K , as the unknown, given K and an approximation 7 7 

of w* . The square root terms of the algorithm can be! forced out of the ateasimaaad -hand 

_ side oft this by writing 


_ which might be regarded a a useful device when K > »K. . Finally, the agree 


behaves like a Newtonian algorithm but does not require as many derivatives. a 
Figs. 3 and 4 show two truss problems to which this ‘‘allowable stress”’ 
algorithm has been applied. In both cases wo? = 1,757 (rad/sec)’, and in both ee 
the convergence is quite good. aa the case 3, convergene 


e is 


| = 
aut | 
(9) 
4 
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71, while Fig. 4 converges to 101, 310. ) What is not oo in these figures 
is that fact that they were computed twice: Once as described above, and once 
inverse iteration with almost identical results. 


to the method of inverse iteration rather than some conned eigenvalue routine, oa 
since inverse iteration can be carried out using a structural analysis computer ; 
- after only minor modifications. In terms of the algorithm stated above, _ 
the only change required, in addition | to performing an inve rse iteration (a structural 
analysis) rather than an exact eigenvalue analysis, is to ‘replace at each step a 


the term Min by its Rayleigh quotient approximation: _ 


> Win) _M 


ny M 
q common 


paper has examined an “allowable. stress’ algorithm for structural 

P optimization with frequency constraints. Several properties of this algorithm 

_ have been developed in an attempt to place it within the framework of classical _ 
mathematical programming. For the two examples considered, inverse iteration, 
which has obvious appeal to to the structural engineer, is seen to perform as 
well as the ‘‘exact’’ eigenvalue routine for the purposes of optimization. —__ 


It is hoped that, — analysis such as this paper, an improved understanding | 
of ‘‘allowable stress’’ 


can be improved and their use extended. — - eee. 


Appenoix |.—REererence 


od Bathe, K.-J. ons Wilson, E. L., ‘‘Eigensolution of Large Structural Systems with 
_ Small Bandwidths,”” Journal of the Structural Division, ASCE, Vol. 99, No. ST6, 


2. Beckenbach, E. F., and Bellman, R., Springer-Verlag, 1971. 
3. Cassis, J. H, and Schmit, L. A., Jr., “Optimum Structural Design with Dynamic 
7 ea Journal of the Structural Division, ASCE, Vol. 102, No. ST10, Oct., 


. L., and Kapoor, M. P., 


AIAA 6, 12, Dec., 


“Rates of Change of Eigenvalues an and Eigenvectors,”’ 


1968, pp. 2426-2429. 
Spillers, W. R., , and Farrell, J., “An Absolute Value Linear Programming Problem,”’ 


- Journal of Mathematical Analysis and Applications, 28, 1, Oct., 1969, pp. 153-158. _ 
6. Spillers, W. R., Automated Structural Analysis: An Introduction, Pergamon Press, 
7. Spillers, W. R., Iterative Structural Design, North Holland, Amsterdam, 1975. ary a 
¥. B. “Survey of Optimization Techniques in Structural | Design,”” 

Memorandum TM-FBR-28-43, Wright-Patterson AFB, Ohio, 1978. 
9. Young, J. W., Jr., and Christiansen, H. M., ‘Synthesis of a Space Truss Based on on 7 
_ Dynamic Criteria,’’ Journal of the Structural Division, ASCE, Vol. 92, No. ST6, Dec., 

1966, PP. 425-442. 

| 


Appenoix Il. —Notanon : 


symbols 


per 


) | 

| 
| | 

| 
| 
q 


vector of member 
= structural stiffness matrix; 
mass matrix; ‘ox Wi 


vector of square of ‘member ‘ ‘allowable”’ 


Square of fendementel frequency. 
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evaluate the strength of old buildings, the number of structural investigations — 

_ of this type is increasing. This paper provides architects and engineers with — be 
a procedure for examining wood joists in an existing structure, and for a Poa : 
reasonable judgements concerning the structural adequacy of those joists. 
' os In this study, a a | wood framed barn constructed in Lewd was examined. Tw wo 


method and a longitudinal st stress “wave . method. "Utilizing ign techniques | to 
develop a procedure for examining the joists, the basic questions are: (1) How 
is joist strength determined from either grading technique; and (2) how do the 


Little has been published in the past several decades addressing the problem — 
of evaluating old wood structures. Lee (14, 15), and Wen, et al. (20) report — ‘ 
using longitudinal stress waves to evaluate the strength of i in situ wood. However, a 


“J 
_ much of the work in this area has relied on a visual assessment of sof Caling 


More recently, dynamic measurement techniques for the evaluation of the 


dynamic approach can take several forms, 
~ longitudinal stress waves. Following is an overview of the elementary aspects 
_ of the longitudinal stress wave propagation theory as it applies to the ie 
= of the modulus of elasticity in wood. ‘ 
A simplification of the governing differential equation for wave propagation 
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STRENGTH OF OLD Woop Joists 
wy 
abrupt 
Recyciing, conservation—the now Deing heard in the building 
9 - field. Old buildings, historic buildings, and some not-too-old buildings are being % 
be an, altered and revamped for new uses. In all cases, the structural adequacy must — 
he evaluate While chitects and ensineers have wave heen called 
— 


between the stress wave propagation rate, the material’ s dynamic modulus one 7 
elasticity, and the mass density: 
which c = stress wave propagation rate (a. jeec); E a = dynamic — 


(psi); and p= mass density (Ib sec ?/in.* = 


| 


LWAVE FRONT PLOTTE 
AT TIME INTERVALS > 


(OR) 


wath te | TRIGGER 


ow we 4 
‘FIG. 2—Schematic of Stress Wave Timing Sytem 
dynamic modulus of elasticity follows from Eq. 1 


E,=c’p .. 
- To determine E,, a stress wave must be imparted to the . member. This can 
7 be done by directly impacting the board, i.e., with a hammer, or by an ultrasonic — 
7 pulse. In either case, the wave travels through the material in much the same 
Be as a wave travels through a pool of water. When an impactor such 
as a hammer is used, a compression wave is induced into the material. This’ 
= results in a a differential pattern of motion in- the wood particles which 


is the propagating dynamic stress wave. 


DECEMBER 1981 
through an isotropic, elastic, homogeneous materi 
| 
| 
| 
q 


The cc wave traveling through the is in nature, 
g i.e., a wave front does not separate or disintegrate at discontinuities. The = 
_ propagation rate is greatest along the major axis of a piece of lumber (i.e., “d 

along the grain). In the transverse direction, propagation rate is slower due 

to the changes i in density ! from earlywood to latewood wood to su ‘summer 


seduces | depending or on 1 the severity of the "deviation. If 0: one were to look at 
the components of the stress wave in the direction of travel and along the | 
direction of the cross grain, a relationship might be established as to the degree 
of resistance a grain deviation would produce. As a wave encounters a knot, 
there. are local grain deviations around the knot the abrupt ‘discontinuity 
at the knot. Fig. i ‘shows a plot « of a wave front passing through a piece of _ a 
_ Structural lumber containing knots. The portion of the wave that contacts the 
discontinuity lags behind the rest of the wave, and, in the case of ak knot, 7 
it mgs itself around the knot and eventually becomes continuous again. _ 
-A common method for obtaining the propagation time is through the = 
of accelerometers ¢ clamped against the member. An accelerometer is a device ~ 
which contains piezoelectric material (usually quartz) that is activated as the 
wave passes the point of contact. Fig. 2 is a schematic representation of a 
typical stress-wave timing system. 
_ The stress wave causes the first accelerometer to produce on electrical signal © 
which turns on the timer. . The electrical signal produced by the second accelerom- 
ec r and the time interval i is stored. ‘With this time interval — 
and the board ae the board’ s apparent E, can be calculated. It should - 
be noted that the value obtained is an average E, over the distance between 
accelerometers; furthermore, it is = an apparent E a because Eq. 2 need 


barn was examined we in 1925, for the College of 

_ Agriculture, Washington State University, Pullman, Washington. It was primarily 
used as an animal shelter on the ground floor, and for hay storage on the 
second floor. The inspection was confined to the nominal 2 x 12 in. (50 x 
300 mm) floor joists in the south bay of the barn where hay storage was believed 

to be the primary use. ones ‘he feild which wert ceed 
_ Two methods of inspection were performed to determine the allowable working — 
stresses for the joists: (1) The conventional visual grading technique; and Q) 
the stress-wave technique. The implementation of the stress-wave technique 

to evaluate in-place strength requires a means of relating the strength properties 
to the stiffness pronerties of the individual members. Through the use of a 

_ linear regression, the dynamic modulus, E,, was adjusted to the static modulus, — 

q E,,, which was ultimately related to an allowshie'e extreme fiber stress in bending, 
ae corroborate the findings of others concerning the correlation between E, 
and E,, a random 30 piece sample of Douglas-fir (Pseudotsuga menziesii), — 
- nominal 2 x 6 in. (50 x 150 150 mm) wood joists were examined in the the laboratory. 


| 
{ 
| | 
| 
| 
; | = governing wave equation, Eq. 2 has been shown to be useful in determining 
E,, for structurallumber (5,89). 
— 
£ 
— 


id — 


2 
ale iit tested by both the stress-wave method, and by static third 


point loads to determine E, and E,, respectively. The regression run on this 
data was the vehicle used to determine the static modulus of elasticity, E,, 


- sample). Independent subsamples of 5, 10, 20, and 30% were > chosen from the 
7 parent sample, with replacement, to allow for comparisons in later analyses. ae 
_ The intent here was to determine an appropriate sample size for this type of *§ 


FIG. 3. Wave System: (a). Timing Unit with 


a of this type, it would been far better to sample 
— 100% of the joists rather than only 50%, but due to time and financial constraints ? ; 


a With the parent sample distinctly numbered, visual assessment of the joists" 


- 


value. Next, stress wave time intervals for ‘es joists were recorded. Two 
; accelerometers were clamped to the joist with a specific distance between them. | 


_ The signal was introduced by impacting a specially designed anvil (Fig. 3(a)) 
a hammer. The start accelerometer was located least 6 in. (150 


— 
23 
— 
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of the member was determined. Small blocks were cut from overhangs | on n about 

_ two-thirds of the samples for determination of density and moisture content. _ 
The moisture content of all joists in the sample was measured by using a moisture : 

: "meter. The densities were compared to values listed in Ref. 21. Evaluation 

of this data results in reasonable allowable design values and some indication 4 
of an an acceptable sample size. The: re 

Data Anatysis Weners 


The data are discussed in two distinct sections. The first section conten 
of data obtained from laboratory tests of 30 joists. 


og 
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of 


1G. 4. —Regression of E, Versus E, (Laboratory Tests) 


f elasticity were found using the stress- -wave equipment, and by a third point — 
loading static bending test. The > second section consists of "the stress wave = 
visual grading data obtained in the field which were used to determine the 
moduli of elasticity and allowable bending stress of the barn joists. matin ofl 
_ Section One.—The wood tested in the laboratory was nominal 2 x 6 in. (50 
al x 150 mm) Douglas-fir (Pseudotsuga menziesii) joists consisting of visually graded — 
No. | and No. 2 lumber. The longitudinal stress wave propagation time was 
taken over a 90 in. (2.29 m) s span for all boards. The density was determined — 
by weighing and measuring overall dimensions of each sample, and moisture 
content was measured using an electrical resistance meter. Since all boards 
had moisture contents of approx 7%, it was not necessary to adjust the mechanical 
properties to a standard moisture content. The dynamic moduli were calculated 


Each board was was statically tested as a joist in third point loading. The cen center-point © 


| 
ISTSTRENGTH 
in hows the timer, accelerometers, and clamps = J 
_ were used. Ihe time for the leading edge of the first pulse to travel between _ Z ; 
i 
| 


‘deflection this used to E,. The E, 
-. were modified to a span to depth ratio of 21, and a uniformly distributed load 

_ A linear regression was performed to relate E, and E,. The plot of the 
data and the the regression | line a are shown in Fig. 4. ‘The Tegression equation is is. Jat 


= 0.101 x 10° + 0.688E, or 
_ in which E, = static modulus of elasticity in bending « on edge; and E, = dynamic 


_ The coefficient of determination is, r? = 0.680. This coefficient of determination oe 


is comparable to to previously reported re results (9,12,19). 
TABLE 1 1.—MACHINE STRESS-RATED LUMBER 2” and Lessi in Tiknoae 2 ag and wiser 
a Vetuse in Pounds per Square inch” ave 


Bending Modulus of “Tension parallel parallel 
Single Elasticity “E” 


same as assigned by ASTM methods to visually graded lumber. — | a 
pil «Fb design values are applicable to lumber loaded on edge. When loaded flatwise, 
the values may be increased by multiplying by the following factors: 


Reprinted under Permission of the Western Wood Products Association. = | ——~ 


Section Two.—The longitudinal stress-wave data, taken in the field, consisted 

of a random sample of approx 50% of the total number joists. All but one 

of the data points were considered valid. This one specimen had already failed 
and in any structural refurbishing would need to be replaced. ‘The stress wave al a 
propagation time was taken over spans of 78, 84, or 144 in. (1.98, 2.13, or 
+t 66 m), depending on the length of the joists or other obstacles. Some of a 
the joists in the : sample had short overhangs from which a specimen could 

2. taken for density determinations without reducing the strength of the floor. 
these specimens, the average density was = 0.586 g/cm’. Using the 


! 4 
1650 | 1900 1,800,000 
2250 2600 | 1,900,000 = | 
2550 2950 2,100,000 205000 80 
Compression perpendicular to grain (Fc 1) and horizontal shear (Fv) values are the a - 
7 


"stress wi wave data, and the average density for the sample, the dynamic modulus - 


elasticity, E,, was calculated from Eq.2, 4 
AD approximate average moisture content for the joists of 15% was measured, y . 
using an electric resistance moisture meter, with the appropriate corrections 
_ for temperature and species as per the manufacturer’s specifications. This value s 
_ high for a covered, unheated structure in the month of vou 
ip =o a second moisture content measurement was made by oven drying pal 
- 3 “a cut ‘from the joist overhangs. This average moisture content | t of | 13. nied 


can 1 be used for practical purposes. The joists in a the ‘barn w were a . mixture of - 
4 Douglas-fir and Western larch (Pseudotsuga menziesii and Larix 


ay 


FIG. 5. —Stress Wave E, upon Visual Grade E, 
- With E,, for the sample determined, E, values were obtained by using th he 
- segression Eq. 3. The E, obtained was for a moisture content of 13.6%, and 
was adjusted to a moisture content of 19%, which is the maximum allowable 
in use moisture content for tabulated design values (7). The average E, was 
1.639 x 10° psi (11.3 GPa), and s, the standard deviation n of th the mirror image 
the lower half of the distribution (16), was 0.186. 
All the joists in the sample were visually graded in order to compare E, 
from the stress-wave method, to that derived from the visual grading method. ; 
Fig. 5 shows that the average stress wave, E,, is less than the E, assigned 
using visual grading criteria, i.e., E, stress wave is conservative, in this case, 
to the visual assessment of ay Itis is possible to utilize these test 


| 

> 


“results t to construct a machine-stress-rated (MSR) grade of —_ ie a floor 
system. MSR lumber is classified into F, — E groupings (Table 1). The F, 
= allowable fiber stress in bending) that is associated with each E class 
is derived from a regression function of bending strength versus modulus of 


There are some physical requirements which must also be satisfied, for a 4 


‘Particular se set of | lumber | ton meet an n MSR grade, €.8-, the average edge 1 modulus 


assigned to that grade, E, (16). A second requirement is that the lower <i 
percentile of the sample modulus of elasticity distribution must equal or exceed | 
82% of the modulus of elasticity assigned to that grade (i.e., E.,, = 0.82 E “3 
(16). The factor 0. 82 assures that the grade being assigned has no more than Pi 
an 11% coefficient of variation. The final physical ‘Tequirement refers to the 
allowable visual characteristics for the grade. This category includes s strength 
reducing defects and appearance characteristics. The strength reducing charac- 

h teristics are of primary importance. Appendix I describes all the visual charac- 

_ teristics which can override the mechanical assessment. Because we are concerned 

_ with existing structures, the researcher may choose to ignore : some of the visual 
restrictions not relating to his particular case. 


Application of the aforementioned criteria to the floor joist data suggests 
that a 2050F,, 1.6E£ (14.1 MPa, 11.0 GPa) MSR grade for repetitively apened 
joists might be assigned for 


AND RECOMMENDATIONS 


_ ‘The intent of this study was to develop a methodology for determining strength 
of wood members within an existing structure. Three alternative methods are 
: proposed. The first, and probably the most common in practice today, is the 
process of visually grading the members and using the grading rules to determine — 
allowable design values. The second and third methods utilize the stress wave 
timing equipment to obtain E, which can be related to E,. These last two | 
methods differ in their resultant allowable stress determinations, i.e., one uses 7 
the MSR grade rules, and the other goes directly to curves relating E, to mode 
of rupture (MOR), (ultimate bending stress). a 
Method 1—Visual Strength Assessment. —This technique requires strict ad- 
herence to a grading rule manual. Size and location of defects must be accuratel 
characterized to ensure against high grading or low grading the material. -_ 
common problem is the subjectiveness of this approach. Experienced perenne 
must do the grading. Frequently the material in older structures can be weathered, 
‘Tough sawn, or both, which makes precise identification of the strength reducing 
_ characteristics difficult. Another drawback of this system is the inability to 
look at all faces and edges of the member. Often the top edge of a joist is 
inaccessible; consequently, one is not able to do a complete job of visual © 
Once the visual grade has been assigned it is a simple matter of reading 
_ the allowable stresses from the grading rules. Using current grading rules, the 
joists i in this : study were a grade mixture, with 95% No. | or better. weg No. 
2 joist and one No. 3 joist were foundinthe sample. 
_ The discovery of the No. 3 joist in the sample penalized all of nie ian P 
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joists, w: one wished to assign dictate stresses to the joists. If this joist 
_ is strengthened, then the joists may be graded as No. 2 and better, with F, 
and E values of 1,450 psi (10.0 MPa) (repetitively spaced), and 1.7 x 10° psi 
(11.7 GPa), respectively. Alternatively, if the two joists below the No. | grade 
(were strengthened, the floor system could be assigned No. 1 stresses for 
_ Tepetitively spaced members, i.e., 1,750 psi (12.1 MPa) for F,, and 1.8 x 10° 
Method 2—Stress Wave—MSR Assessment.— Utilizing the stress-wave tech-— 
7 _ nique removes the personal subjectiveness from the process. It gives a mechanical _ f 
/ means of characterizing defects. A visual override must still be employed along 
_ with other limitations that are listed in the grading rules (11). Members with 
- obvious failures must be noted and replaced. The stress-wave technique is also — 
_ valuable for indicating defects not readily visible, such as incipient decay, sloping 
grain, and insect infestations. The extent of decay in large wood arches by 
 stress-wave techniques has been reported (13). 
When using stress-wave equipment, one must he aware of the fact that a_ 
hole or a partial failure (horizontal split) in the member will have little effect 
on the wave propagation rate. The wave will tend to travel in the sound es iy 
essentially ignoring this. type of defect. It should also | be noted that the E, 7 


derived from these types of tests is an average E,, over the distance between - 

accelerometers. If necessary, the accelerometers can be placed at short intervals 
_which would allow the determination of a low point E, value. The stress-wave _ 
assessment gives the investigator a feel for the distribution of material properties. 7 

oA procedure that might be followed would begin with calibration of the | 

- stress- wave equipment to confirm the E, versus E, relationships. Inspect all 

as members for decay, failure, and other limitations. The writers recommend that — 

all joists that are decayed, have men or are below the desired grade be 


“ Upon elimination of the edu failed specimens and those with visual 


be determined. Species can ‘be confirmed by a a local Wood T echnologist, — 
_ US. Forest Products Laboratory, or through an identification key in a basic 
wood technology book (6,18). The most accurate methods of density and moisture 
me determination would be with actual wood blocks removed from some “ 
safe location. If this is impossible, the writers found that Ref. 22, and the i 
electrical resistance meter both gave conservative measures of these properties. + 4 
A note about densities taken from Ref. 21 should be made here. Handbook _ 
- values which most appropriately characterize the material being investigated _ 
_ should be used, e.g., a structure constructed in 1925 would be built with material 


nearest that date should be used. This is done in consideration of our changing | 

resource base. The densities from Ref. 21 are conservative because the bulk _ 
om (wood fiber and moisture) (10) is used in Eq. 2, while densities in 

the Wood Handbook are of the wood only. I, vaereweVormasy corrections can = 


from a pre-1925 resource consequently, the edition of the Wood Handbook 


“This i is obtained as described in the procedure section of this report. Using 
-. = this information, E, is computed and modified to an E, value as shown in 
Fe 4. = sample can then be analyzed according to the procedure previously 


| 
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outlined, an and a a MSR grade c can be assigned. 
ih recognition of the fact that MSR machines monitor the low point, E, 7 7 
of each piece, as well as its average, it may be necessary to provide some — 
"assurance that the sample does not have an unrealistic E value based only 
on the average E obtained from the stress-wave technique. If the low point 
cannot be determined, a conservative approach to’ “assigning F,, values would 
be to > assign a MSR based on a minimum or near minimum E, value 
& from the sample. For this study one might address this by applying a nonparametric 


to the the second order statistic for 


. Tables ‘exist such that one can select the enter statistic, 


ow 
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i yee FIG. 6.—Confidence Limit for Mean Value of MOR given an Release 


- with its associated confidence level commensurate with the objectives of the 


a There are, however, two inherent problems with utilizing the MSR approach. 
?— The first is that the MSR allowable stress values include a size correction which 
_ is monitored at the sawmill through daily in-plant quality control procedures. " 
For visually graded lumber, the size effect is incorporated into the published 
= values (2), but for MSR lumber the size effect has not been explicitly ~e 
identified. A second problem relates to the visual grading requirements which © 
may override the grade as determined by E alone om A suitable correction — 


| 
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Nonetheless, the MSR table for are based. on an adjusted lower 5% 
level, using the reduction factors specified in Ref. 2. 
Method 3—Stress ‘Wave-MOR Assessment.—This method is a more empirical 
_ approach not relying upon the grading rules manual. This technique relies solely m4 
~ upon the grading rules manual. This technique relies solely on the relationships _ 
between E, and E,, and E, versus MOR. For method 3 to be used by the ied 
or engineer, further refinement is mecessary. Allowebie stresses must 


be derived. This ci can n be done i in n the following manner. aa 


regression line fitting the data by ‘minimizing ‘the sum of squared residuals 
Fe ste about that line. Often it is useful to compute confidence limits 

; for the mean response (in our case MOR), given a selected independent variable, 
_E,. Pictorially this can be venunears by Fig. . 6. Using ro 6, it is — 


DISTRIBUTION OF MOR 


ABOUT 


et 


AB TE 7. —E, versus MOR for Common Species of Structural Le 


ok A nice feature about this method is that the — can ogee a = 


_of confidence suitable to his particular needs, i.e ., it may be necessary to choose 


8 very high degree of confidence in his analysis; therefore, he may choose _ 


a Another interesting feature about this mathad is seen upon noting the a 


‘master « curves,”’ i.e., , one or two 
curves for a all structural species, , might be | developed; thus, greatly simplifying 
_ the design process. The allowable bending stress a rom m method 3 can be compared 


| 
| 
a | are only slight differences in the slopes and locations of the lines for the six 
i — 
a 
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which X, = = 6,300 psi at =" 4 x 10° (Fig. 7 
7 -~ipsi (17); and ¢, ,. = 1.645 at the lower 5% exclusion level (normal distribution). 
_ The 1.15 factor is the allowable stress increase for repetitive member use. F, 
is found to be 2,000 psi (13.8 
The writers feel method 3 is of special interest because of the direct nature 
‘ in which allowable stresses can be derived. We feel further development of | 
ts technique is warranted and encourage additional investigation toward this _ 
g Seecntiiete for estimating the strength of in situ structural members are i in 
demand because of increased restorations, refitting, and repair of existing 7 
structures. This initial study presents procedures for estimating the strengt 


h 

of repetitively spaced wood joists, based on on a ty measures and the 

modulus of elasticity. 4 


E of Design Values 


F, for Repetitively 
in pounds 


| per = | in mega- n 
square inch in Remarks 


(6) 


71% confidence 
of E values 
equal or the 


second order 


Statistic 


050 «(14.1 | 1.6 x 11.0 by machine-stress- 
ol ¥ on rated procedure 


Table 2 summarizes the allowable values by the different 
methods. Method 2, 
to or larger than method 1, the eal strength as assessment. . Method 3 yielded 
‘comparable results but needs additional 
_ One important advantage of stress wave measurement is the ability 1 to detect 
the low strength members. These members can be reinforced or replaced and 
the overall structural capacity can be increased. Phas 
_ After eliminating failed joists by visual inspection it was was found that a random : 
= sample of 10% of the joists was statistically adequate for this investigation. a 
a However, it may be impossible to randomly select the specimens for the sample. 
For practical purposes all joists need to be examined, whenever possible, in 


‘The nature of the Nondestructive e Testing (NDT) stress-wave inspec- 


— 
| J 


| 


tion the subjective: necessary procedures res are re- 
_ quired. Future work is needed to refine the instrumentation into a convenient _ a 
package that can be used easily in the field by investigators. The method of me 
a arriving at allowable stress for in situ joists needs further study. Future test 
_ programs should include field testing and then the removal of joists toa laboratory — 
for st static > testing. T The results will provide direct correlation between and 
_ A question of greater concern in a timber structure is the ‘strength of large 
: "primary members—solid beams and glulam beams. Solid beams are found in 
_ ~pre- -World War II structures, while glulam beams have had increased use in s 
- post-World War II structures. Is the response of large members to stress waves 
similar to the response of nominal 2 in. (50 mm) joists? 
> It is obvious that no one method is a panacea, but it is felt that the e stress-wave 


4 techniques show considerable merit and deserve further 


x The writers wish to thank R. Pellerin and R. J. Hoyle, Jr., Washington State a 
University, Puliman, Wash. for their comments and suggestions during this study a 
The writers also wish to acknowledge the constructive comments wl the ASCE A 
_ reviewers, including Donald M. _ Onysk of Forintek Canada Corp., , who were 
neat first writer wishes to » acknowledge the receipt of a Professional Develop- _ 4 
ment Grant from the National Science Foundation which enabled him to take i : 


_ Mechanically stress-rated lumber must be well manufactured and visually 

graded to limit certain characteristics even though the actual strength is not. 


affected. All pieces | shall be Visually graded to assure that the characteristics 


as Checks—Sessoning checks not limited. Through checks at ends limited as 


a hake—On ends limited to 1/2 the thickness. Away from ends ak 
shakes up to 2 feet long, well separated. If not through, single 

_ Shakes may be 3 feet or up to ‘1/4 the length, whichever is greater = 
"Skips—Hit and miss, and in addition 5% of the pieces may be hit or miss _ 


or heavy y skip n not longer than 2 feet. Section 720. 0.000(e), 


- Splits— —Equal in length to 1 1/2 times the width of the piece. its 


- Wane—1/3 thickness, 1/3 width, 5% of pieces may have wane up | to 2/3 
thickness and 1/2 the width for 1/4 the length. 


In addition to the visual limitations listed, knots, knot holes, burls, distorted a : 
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grain or decay partially or wholly at edges of wide faces, must not occupy © 
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velocity; 
a E = allowable design modulus of elasticity; 
modulus of elasticity assigned grade; 
"static modulus of elasticity in bending on edge; 
average edge modulus of elasticity; pet 
= allowable design extreme fiber stress in bending; niet “ 
standard deviation of mitror image of lower half distribution; 
material mass density. ah 
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UNDER WIND DISTURBANCES pia 


Yo By Kiyotaka Morisako,' Shuzo Ishida,’ Tsuneyoshi 


— the evaluation of dude wind loads, a considerable wanes of attempts | 
hs been made to take into account the dynamic behavior of structures when | 


subjected to fluctuating wind pressure. In the latter half of the 1960’s, procedures 
for evaluating design wind loads for structures as a function of their mass 
distribution , flexibility, damping, and size were proposed by, e.g., Davenport 
spectral approach for the estimation 1 of the response of a linear system. 
_ They differ in detail, e.g., the estimation of spatial characteristics of the turbulence 
by Vellozzi et al. is different from those in the aerodynamic admittance by — 
_ Davenport and by Vickery. Each of them has proposed a factor defined as 
_ the ratio of the maximum expected value of a structural displacement response 
to its mean value. The factor has been called “gust loading factor (2),”’ or 
4 “gust response factor (14).’” Many studies (e.g., 6, 7, 13) have since been carried 
out with the intent to place the evaluation of design wind loads on a more 
_ Only two papers have dealt with the plastic behavior of structures under 
Fe "wind action. Vickery’s investigation on the behavior of simple elastic-perfect 
‘plastic: Structures is an approximate linear spectral approach and includes 7 
numerical analysis (16). He has proposed ‘ “the damage rate defined as the rate ma, 
of accumulation of permanent set.’’ The * ‘rate of accumulation of plastic damage,” ” 


¥ similar to the damage rate, has been examined by Wyatt and May ona bilinear 


— 


oe In order to establish a rational plastic design method of wind- resistant buildings, , = 
itis s indispensible | to accumulate a large number of data on dynamic « haracteristics — 
> 0 of wind, and to accurately investigate dynamic collapse behavior of structures. ‘i 


= 'Asst., Dept. of Architecture, Kyoto Tech. Univ., Kyoto, Japan. 
*Prof., Dept. of Architecture, Kyoto Tech. Univ., Kyoto, Japan. eth bosaces how 
hal *Prof., Dept. of Architecture, Kyoto Univ., Kyoto, Japan. a at . 
*Prof. * Disaster Prevention Research Inst., Kyoto Univ., Kyoto, Japan. 9 e.. 
-Note.—Discussion open until May 1, 1982. To extend the closing date one acu, _ 
a written request must be filed with the Manager of Technical and Professional Publications, . 
_ ASCE. Manuscript was submitted for review for possible publication on September 3, 
1980. This paper is part of the Journal of the Structural Division, Proceedings of the | 
American Society of Civil Engineers, ©ASCE, Vol. 107, No. ST12, December, 1981. 
0044-8001 /81/0012-2365/$01.00. 
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, the of d on a strong wind, and the development 
of the methods. of dynamic large- -deflection analysis, enable one to carry” out 
numerical analysis of dynamic collapse behavior of multistory frames. subjected | 
‘to strong winds. Some characteristics of the dynamic collapse behavior of | ad 
_ strain-hardening frames have been enslyzed by the senior writers (3,5). The 
_ numerical method used by them has been developed on the basis of the combined 
“nonlinear finite element method, and of the combined use of the appa matrix 
This paper clarifies what | strengths are to be assigned to the iene of a 
a frame so that it can withstand a fluctuating wind disturbance. _ The effect | 
of a wind disturbance on the collapse of the frame is estimated with respect - Bate 
to the following two scales on the overall strength of the frame (Fig. 1): (1) J 
a by means of simple plastic theory; and (2) 8 by means of the refined numerical 
ervey The concept of the “dynamic collapse gust load factor’’ 


by the senior writers (5) is examined more closely. A plane six- story frame 


__ STATIC 


by a refined 
4 numerical analysis 


etude FIG. 1 —Numerical Approach to Dynamic Collapse Gust Load Factor ‘diveda 


of one bay is adopted here as a model frame. Model wind disturbances are 


ag 
‘the numerical approach, as shown in n Fig. 1, several frames with different ‘overall ‘ 
strengths are tested consecutively as to whether each will survive or collapse 
= a model disturbance. From this approach, a critical frame with the boundary 
strength is found. This boundary strength provides a rational basis for a safety — 
- factor for the frame to o be designed with respect to the wind disturbance. . The 
a more sound foundation for determining the safety factor of structures subjected 
Model Disturbances Wind Fluctuation. — Fora simulation 
= of dynamic collapse processes, it is desirable to apply realistic wind — 
disturbances to frames. The fluctuating wind pressure record of Typhoon No. 


) 
| 
| 
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on 22, 1972 (4), been utilized here. This record was taken at 
7 a height of 10 m above ground in the Shionomisaki Wind Effect Laboratory _ 

of the Disaster Prevention Research Institute of Kyoto University. The maximum 
z: 10 min mean wind pressure was 126 N/m? (12.8 kgf/m_’). In the present study, 
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the or original record was so amplified that the resulting 10 min mean wind velocity 
at a height of 10 m above ground would amount to 30 m/sec. It is assumed, _ 
for the sake of this numerical analysis, that the windward wall is subjected 


to the wind pressure in proportion to (h/h,)°° in which A, and h represent 


ec: 
Bea 


‘the total wind pressure on the windward wall is uniformly distributed over 


= the numerical analysis, three segments of the amplified wind pressure 
record have been utilized as model disturbances. Each segment represents a 
_ typical fluctuating pattern which can be frequently found in this wind pressure 
record, and which would give the highest mean wind pressure among similar 
segments. The three model disturbances are shown by solid lines in Fig. 2. 

‘The data have been digitized with the sampling interval of 0.04 sec. Run No. 

1 is characterized by a peaked wave, and Run No. 2 by the repetition of waves - 
with fairly high amplitudes, and with a period of about 2 sec. Run No. 3 is 


characterized by a high wind pressure with small factuations 


oo 


Basic Design Wind Loads 


and Working Gravity Loads 


(a) Geometric Data and Wind Load Ratios to 
System Coordinates the Reference Value 


_—— of Model Frames.— —In view of the purpose of the present nonlinear 
‘dynamic analysis, fairly tall, slender frames are desirable as representative a 
particularly due | to the possibility of their involving greater P-A effect. For — 
the purpose of parametric study, however, it is not desirable to involve too = 
“many members, since the hysteretic behavior of all the elements would have — 
to be analyzed and stored in the computer memory. The geometrical data shown — * 
in Fig. 3(a) have, therefore, been chosen here. 
Several frames, with different overall strength, have been designed consecu- 
tively to find the boundary strength, and each frame has been tested numerically 
_ for survival or collapse under the model disturbance defined here. The nodal 
_ force introduced by the pressure, corresponding to the 10 min mean wind velocity 4 
s 30 m/sec at a height of 10 m above ground, has been adopted as the reference 7 
_ value of the wind loads. The wind-load ratios to the reference value at each | 


node | are shown in Fig. 3(b). These ratios have been estimated on the ensle 
are 


n 
fees the spatial distribution assumed here. The “basic design wind loads’’ 
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defined by the reference and the wind-load ratios. Fig. 3(c) shows the 
basic design wind loads, together with the working gravity loads. Each frame ; 
_ has been designed using simple plastic theory so as to support the two sets 
of loads: One, consisting of the working gravity loads, and a times the lateral 
loads shown in Fig. 3(c); and the other, 1.4 times the working gravity loads 
_ only. _ The number, a, is called the ‘‘design load factor’’ in the sense that it 
represents the gust level for design, as compared with the 10 min mean level. 
> TABLE 1.—Cross-Sectional Areas and Fundamental Natural Period of Frames on 


Cross-Sectional Areas, in square centimeters 


period, 
| Fourth | Fifth | Sixth 
| @ | | ® | ® 
$4.29 97 | 12.54 
Column | 116.33 . 29.56 
Girder | 119. 12.54 
Column (116. 29.56 a 
Girder f 12.54 
29.56 
12.54 
29.56 
12540 


29.56 = 
A frame ‘designed for a a : times s the lateral loads Fig. 3(c)) is called FRa, eg. , 
the frame for a = 2.70 is designated as FR2.70. 
_ The required fully plastic moments of the members of each frame have been ; 
found under the two loading conditions by assuming over complete collapse 
_ mechanisms (10,19). The in-plane cross-sectional properties of a member may = 
; - be represented by those of an equivalent idealized sandwich member (19). The i 
equivalent depths of all the beams and all the columns adopted here are 40 
om and 30 cm, respectively. For a prescribed value of the yield stress, _ 
% area of an equivalent flange may readily be found SO as to satisfy the interaction — : 


Table 1 shows the cross- sectional areas sof the ‘members of several frames 
 FRI.00-FR2. 80, and the fundamental periods of their frames. All the other _ 
values have been omitted since they turned out to be almost equal to the values — 


The method of numerical analysis developed in Refs. 3, 12, and 18 has been 
applied to the present wind response analysis. The following approximations 
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1. A cross-section is idealized as an equivalent sandwich section consisting 
of two flanges, each possessing a concentrated area of half the actual total “= 

area, and located at the radius of gyration. 

Ss 2. Each flange obeys the bilinear hysteretic stress-strain relation, with Young’ s 
modulus | of 206 kN/mm’? (2, 100 ton f/cm? ), the yield stress of 2.35 N/mm? 
(2.4 ton f/cm’ ), and the strain- hardening | modulus of 0.005 times Young’s 
3. Each member consists ; of several beam-column elements connected in series, 
and only the representative, or averaged flange stress and strain of an element, 
- the relation in assumption 2. The basic element length used here is 40 
‘a Only the buffeting response in the plane of a frame is analyzed by employing 
the system coordinates shown in Fig. 3(a). The following approximations have 
4 been introduced fordymamic analysis, 
_ 4, The generalized mass matrix consistent with the displacement functions 
y of a member in the initial elastic range has been generated and held constant. _ 
_ §. The viscous damping matrix is 0.01 times the system stiffness matrix, 
4 _ and the aerodynamic damping has been neglected in comparison with the former. 7 
a A modified version of the Wilson-Clough method (1) has been utilized 


for direct integration of the equation of motion. After a convergence — 
‘the tim time step of 0.01 sec has been here. 
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a The lateral loads proportional to the design wind loads ety been applied © 
to the several frames and increased under the constant gravity loads. Fig. 4 _ 
_ shows the base shear-top lateral displacement diagram. The dropping portio 


f each curve has been obtained by the -controlled 
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analysis (il) and demonstrates the negative aie reflecting the effect « of 4 
4 » sy The ratio, B, of the maximum base shear to that due to the basic design 
Se, j wind loads is called the static collapse load factor. The values of B have been = 
= indicated by the numbers in the round brackets in Fig. 4. 4. The value of B wil 

“—1 each frame is slightly greater than its a value. In general, the B value for — 

4 a frame may not be equal to its prescribed « a value. 


The problem here i is to clarify what strengths a are to be assigned toth the cue rs 
uctuating wind disturbance during a “J 
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disturbances. When this disturbance is terminated, the frame is brought under 
. _ the mean pressure level again. Numerical simulation of this process was conducted 
on all the frames designed in the preceding section, and the dynamic collapse _ 
behavior and survival behavior of the frames were investigated. Some frames a 
may or may not exhibit the critical or precarious behavior during t the he application — 


a 
50. Wh 150. 200. -25. 08 so. 150 200. 
| 
te 
= _TOP LATERAL 
oo _ windstorm. For this purpose, the standard state of a frame is defined as that 7 
a state which is subjected to the gravity loads and to the mean wind pressure. 
_A ame_unde he andard state is subiected to one of th hree model __ 
— 
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would collapse during the autonomous state after the fluctuation was terminated, 
the numerical analysis of the behavior was - continued for 10 sec under the | 
“mean wind pressure. The analysis was also terminated top lateral displace- 
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FIG. 6.—Phase Plane Diagrams (Run No. a 4 


4 While frames have been subjected here, tentatively, to only three ore wid 
_ disturbances, it is necessary to conduct a similar numerical analysis with respect u 
to more realistic disturbances. Therefore, this work is viewed as a first step “A 
is toward the development of a more sound foundation for determining the safety 
factor of structures subjected to wind loads. 


Phase Trajectories. —Figs. 5, 6, and 7 show the phase trajectories on the 
top lateral velocity-top lateral displacement in these the symbol 


J | 


oun the time when the wind pressure ‘fluctuation is ov over, and each 


along a trajectory indicates the time in seconds. | 

_ The dynamic behavior of simple elastic-plastic structures has been studied ; 

__ previously with respect to the phase trajectories (8,20). Sun et al. (8) studies — 
the gravity effect on the dynamic behavior of an undamped single degree-of- _ 
“freedom inelastic system in free vibration, subjected to given initial conditions, 

and in forced vibration to an earthquake input. Yokoo et al. (20) has also analyzed — 
a definition of the dynamic stability limit of a simple bilinear hardening column 

- — and the engineering significance of stability under a finite disturbance. sui 
ae: ‘In this paper, phase trajectories will also play an important role in studying ~ 

- the dynamic behavior of frames subjected to model disturbances. But it would a 

be extremely difficult to investigate the dynamic behavior given here with the 


of Figs. 5-7 will therefore have to be attempted first. These trajectories are 
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[ (1) A process represented by a spiral in the neighborhood of a point on the 
displacement axis; (2) a process represented by a segment from one spiral to 

4 the next; and (3) a process such that both displacement and velocity components 

4A of the trajectory exhibit a tendency of divergence within the region of the 


response given here. The three processes correspond to the following patterns 
(1) Elastic Vibration-like Behavior.—Process 1 represents an approximate 

elastic vibration of the frame virtually without the work dissipated for the plastic 
deformation. . The pattern of the behavior of the frame represented by this process — 

will, therefore, be called ‘‘elastic vibration-like behavior.”’ | 

a (2) One-way Bounded Plastic Deformation Behavior.—Process 2 means thet 

; some plastic regions have been developed in some parts of the frame, but 

the one-way still bounded. The pattern of the 


| 
q 
_ analytic approach used in the preceding studies. A phenomenological examination te 
> 
1 
all 
substantially comnosed o me or gi) t fol ino three tyvnical nrocesses:_ 
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of the by this process, followed by process will therefore 
_ be called ‘‘one-way bounded plastic deformation behavior.” 
(3) One-way Unbounded Plastic Deformation Behavior.—Process 3 represents 
. weeny in which the plastic deformation of the frame, and the corresponding — 
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-s 10. —Time Histories c of Top and 4th Floor Lateral Di lacements (Run No. 3) 
‘nodal displacements, will keep o: on increasing unboundedly during the time interval — 
considered here. The pattern of the behavior of the frame represented by this 7 
= will therefore be called ‘‘one- wary unbounded plastic deformation beha- 


7 


- Run No. 1.—The trajectories of FR2.80, FR2. 75, and FR2.72 (Figs. S(a){c)) 
: mie represent the elastic vibration-like behavior. On the other hand, the 
trajectories of FR2.71, FR2.70, and FR2.60 (Figs. 5(d)-(/)) eventually represent 
the one-way unbounded plastic deformation behavior. Even while the wind 
4 ‘pressure varies at a relatively low level after 10 sec in (Fig. 2(a)), the trajectory . 
of FR2.71 represents the one-way unbounded Plastic « deformation behavior. The 
trajectories of FR2.71 and FR2.70 exhibit a c cusp-like process w with a a small loop, 
- followed by Process 3, after the Process 2, from 3-7 sec. This process indicates 
_ that there exists such a point as a saddle point for an autonomous system. ee by 
4g Run No. 2.—The trajectories of FR2.80, FR2.70, and FR2.65 (Figs. 6(a)-(c)) 
eventually represent the elastic vibration-like behavior. On the other hand, 3s 
aie of FR2.64, ‘FR2. 63, FR2. 60, FR2. 50, and FR2. 40 (Figs. ene 
the trajectories of FR2.64 pow FR2. 63, ‘there also exist such ‘cusp-like processes 
: as those observed in the trajectories of FR2.71 and FR2.70 (Run No. 1). FR2.60 : 
exhibits the elastic vibration-like behavior twice, while FR2.40 does so only» 


4 


SEC SEC 16.8 SEC 
—Deflected Shapes and Plastic Regions FR2.64 Subjected | to ‘Run 
once. The one loop process of FR2.50 is a degenerated version of the second — 
elastic vibration-like behavior of FR2.60. Among those frames whose a values’ 
are between 2.40 and 2.50, there will probably exist a frame whose trajectory _ 
exhibits a clearer cusp-like process as a limiting case of the one loop process 
_ Run No. 3.—The trajectories of FR2.40 and FR2. 31 (Figs. 7(a) and (b)) 
eventually represent the elastic vibration-like behavior. On the other hand, the > 
trajectories of FR2.30, FR2.29, and FR2.20 (Figs. 7(c)-(e)) eventually represent — 
the one-way unbounded plastic deformation behavior. The trajectory of FR2.30_ 
is mostly similar to that of FR2.31 until 21 sec, but the displacement component — 
of FR2.30 thereafter appears to he on increasing. In the autonomous State 


response nodal velocities have been very small compared to the mean | 


= 
@ 
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Displacement Response. —Figs. 8, 9, and 10 time of 
- fourth floor lateral displacements. The broken line ia each diagram indicates _ 7 
the time when each wind pressure fluctuation is over, 
3 The pene: characteristics have been observed in these diagrams. oat 
Even if a frame exhibits the one- way unbounded plastic 
 Sieviee. the lateral displacement of the fourth floor hardly increases. diiliens 
_ 2. A set of displacement response curves of the top floor looks as if ‘bundled io 
threads were becoming loose. For instance, the curves for Run No. 2, shown 
in Fig. 9(a), have been bundled initially, but they are divided into two groups — 
: at about 6 sec. One is composed of the curves whose displacement components a 


appear to keep on increasing rapidly, and the other is , composed of the subsequent _ 
gs -shaped curves. The latter is further divided into two groups. 


_ 3. Once the top lateral displacements exhibit a tendency of divergence, withia 
several seconds they exceed the prescribed displacement limit of 5 nm. lien a 


Overall Deformation and Growth of Plastic Regions. shows pr 
, drawn at several 
second intervals i in the s same scale a as the frame geometry then any amplification 
of displacement scales. A line drawn along a deformed member indicates the 
extent of the plastic region developed in the flanges of the member. The overall 
_ deformation proceeded almost with the fundamental mode. The plastic regions — 
- developed initially the third story through the sixth, and Subsequently in the 
fifth story. Eventually, the overall deformation involved some tendency of a 
_ story-sway mechanism. The maximum averaged flange strain at 18.4 sec is 9.3% 
in tension at the windward flange of the lower end element of the windward © 
column in the fifth story. The one-way unbounded plastic deformation processes — 
of frames are roughly similar. A typical process is shown i in n Fig. 11 for FR2.64. — 
o Efficiency of Numerical Analysis.—The computation | was carried out on a 
_ FACOM- M190 computer at the Data Processing Center of Kyoto University. 
_ The time required for one case was about 300 sec, and the memory required © 


(tend 
Dynamic Coutapse Bewavior AND Co.tarse Gust Loao F Factor 
In view of the numerical results, frames designed for various design load 
factors, a, are divided into two classes according to the dynamic behavior under 
a fluctuating wind pressure disturbance. One class of frames is characterized 
by the elastic vibration-like behavior at the final stage of the analysis. The 
- other is characterized by the e one-way unbounded plastic deformation behavior 
within the duration of the analysis. The latter is unbounded in the sense that 
both lateral displacement and velocity components of at least one story exhibit 


a tendency of divergence within the region of the response given saeane : : 


all the displacement components of a frame will be eventually bounded after 


the frame has completely crushed or overturned, the trajectory must | also be : 
_ bounded. If the analysis were continued a 7 


turn out to be bounded. For practical purposes, however, ‘theo one- -way a tion ol an 
ee deformation behavior may be regarded as a collapse behavior, possibly — 
leading to a final state, and is called the dynamic collapse behavior. On the A 


| 
Other nana, © ciasuc ViOrauion-like DOenhavior at the final Stage OF analysis 
| 


FRAME COLLAPSE 


“fl appears to: represent a survival behavior against the wind disturbance. 
This classification of the nonautonomous systems, with respect to their 

trajectories, roughly corresponds to the well-known classification of the trajec- 

tories of an autonomous system, and is supported by the phase trajectories — 
_which almost indicate the existence of such a point as a saddle point for the 
autonomous system. It is natural to expect, therefore, that there also exists 

a watershed trajectory dividing a saddle-like cluster of trajectories, for the present — = 
nonautonomous systems. ‘This potentially existent saddle point may be called 

a dynamic collapse point. ‘It is assumed that there exists a frame whose phase 7 
trajectory passes through such a dynamic collapse point. The corresponding 
value of a, denoted by a., may be called the dynamic collapse gust load factor 
of the frame to be designed with respect * the fluctuating wind pressure © 

disturbance. The frame designed for a. Value is a critical frome 


with respect to the overall plastic strength of their frames. The design inated 
factor a due to the simple plastic theory does not represent the true strength - 
of the frame. The static collapse load factor B, evaluated by the present refined 
numerical method, will provide a more accurate estimate of the effect. The — 
B value of the critical frame designed for the a, value is denoted by B.. Then — 
the value of a. may be called a nominal value of the ‘dynamic « collapse “gust 
load factor. However, for the purpose of the initial design, the nominal value 
Even if such a dynamic collapse point is not found numerically, it is possible y 
_ and, for practical purposes, sufficient to determine the upper and | lower bound = 
values of a, from a values of the neighboring phase trajectories. . The ¢ highest _ 
" lower bound and the lowest upper bound are denoted bya’ anda” , Tespectively, 
and ae a. < For the present, three numerical values of 


a* and a” are given by: a‘ = 2. Tl, = 2. 72 for Run No. 1; at = 2. 64, 
a, =2. 65 for Run No. 2; and at = 2. 30, a’ = 2.31 for Run No. 3. ‘Similarly, 
for a more accurate value, B., , of the dynamic collapse gust load factor, values 
of and BY are given by: = 2.90, Bo = 2.91 for Run No. 1; Bf =2.8 80, 
¢ = 2.82 for Run No. 2; and B £ = 2.40, 8.” = 2.41 for Run No. 3. is seheath 
24 The solid straight line and the broken straight line in Fig. 2 indicate the 
pressure level corresponding to the value and the a, value, respectively. 
_ The peaked wave characterizing Run No. 1 goes beyond the solid line. The 
__ Pressure level corresponding to the 8” value for FR2.72, which has not exhibited — 
; Z dynamic collapse behavior, ri is slightly above the solid line. The external 
- forces represented by the peaked wave do not have sufficient intensity ) and 
duration for inducing the dynamic collapse behavior of FR2.72. The frame has — 
therefore exhibited the elastic vibration-like behavior against the subsequent 
wind pressure fluctuation of lower level. In this sense, the increment appears 
to represent the toughness of the f: frame. Run No. 2 represents the wind fluctuation 
- obviously below the solid line. The decrement appears to represent the unfavorably 
_ dynamic effect of the wind pressure fluctuation on the frame. In Run No. 
43 the peaks of the wind disturbance are almost on the same level as the solid - 
line. In this case, it has been difficult to observe the effect of the wind pressure — 


| 
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"introduced for purpose of estimating the effect the resonant amplification 
of structural displacement response to gusts in simple static design criteria. wW 
5 The calculation involves a spectral analysis of the dynamic response of a linear 

7 system. On the other hand, the dynamic collapse gust load factor of a frame, 

_ with respect to a wind pressure disturbance as defined here, reflects the effect 


eS wind pressure fluctuation on the dynamic 4-4 behavior of the frame. 


_ A parametric numerical analysis has been conducted on strain-hardening 
multistory frames subjected to wind disturbances generated by amplifying a 
_Tecorded fluctuating wind pressure. The effect of spreading and diminishing 
_of the strain-hardening regions has been taken into account in the one-dimensional 
‘finite element formulation employed here. The following facts should be noted: 7 
a ' The dynamic behavior, including, in particular, the dynamic collapse 
_ behavior, has been extensively examined and the details are presented, for 7 
the first time, of strain- -hardening multistory frames | ‘subjected to fluctuating 
wind pressure disturbances. A phase trajectory is composed of some or > 
of the following three typical processes: (1) A spiral trajectory represents the 
_ elastic vibration-like behavior of the frames; (2) a segment of a trajectory from 
one spiral to the next represents the one-way bounded plastic deformation — 
behavior; and (3) a divergent trajectory within a given duration represents — - 
7 one-way unbounded plastic deformation behavior of the frame. The one-way 
_ unbounded plastic deformation behavior has been regarded as a dynamic collapse 
~ 2. Under a fluctuating wind pressure disturbance, frames designed for various — 
design loads are divided into two Classes: one in which frames exhibit the elastic 


“other in frames exhibit the or one- 1e-way unbounded plastic: def ormation behavior 
within the given duration. This classification of the nonautonomous systems 
with respect to their trajectories roughly corresponds to the well-known classi- 
fication of the trajectories of an autonomous system. There exist a couple of 
frames whose trajectories almost indicate the existence of such a point as a_ 
g _ saddle point for the autonomous system. The trajectories of those frames ; appear 
to be quite close to the watershed trajectory dividing the elastic vibration-like 
_ behavior and the one-way unbounded plastic deformation behavior, 
ia 3. The concept of the dynamic collapse gust load factor proposed by the 
_ senior writers (5) has been more thoroughly examined with respect to the three 
typical patterns of a recorded wind disturbance. The effect of a fluctuating 
4 wind disturbance on a frame may be represented by the required overall plastic 
strength, and therefore measured by a multiplier (load factor) of the basic wind» 
loads. There exists a critical frame which just exhibits the dynamic oulleges 
__ behavior under a wind disturbance. The corresponding critical value of the _ 
multiplier has been called the dynamic collapse gust load factor. There are 
two methods for evaluating the dynamic collapse gust load factor. The factor 
| evaluated by a refined numerical analysis is more accurate than the nominal Pd 
value due to the simple plastic theory. However, for the purpose of the initial - 
design, the Bominal value will be more convenient in practice. The dynamic 
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is, therefore, a more rational method of wind resistant 


‘ Clough, R. H., aint of Structural Vibrations 2s and Dynamic Response,” Recent 


_ Advances in Matrix Method of Structural and Design, R. H. Gallagher, Y. Yamada, J 
and J. eds. The ‘University of Alabama ‘Press, Alabama, is7l, PP. 


A.G., “Gust Loading Factors,” Journal of the Structural Division, ASCE, 7 
"Vol. 93, tag ST3, Proc. Paper 5255, June, 1967, pp. 11-34. 
4 3. Ishida, S., “A Method of Static and Dynamic Large Deflection Analysis of Elastic-Plas 
a tic Frames,” * (in Japanese), thesis presented to Kyoto University, at Kyoto, Japan, * 
7 in 1975, in partial fulfillment of the requirements for the degree of Doctor of Engineering. © 
4 Ishizaki, H., Murota, T., and Yoshikawa, Y., “‘On the Behavior of the Glass Pane 
under Natural Wind,” > Gn Japanese), Proceedings of the 2nd Symposium on Wind © 
Effects on Structures in Japan, 1972, pp. 67-72. 
. Ishizaki, H., Nakamura, Tsuneyoshi, Ishida, S., and Okuhira, T. , “The Analysis 
of Dynamic Collapse Behavior of Strain-Hardening Multi-Story Frames Subjected 
to Strong Gusts,”’ (in Japanese), Proceedings of the 4th Symposium on Wind Effects” 
on Structures in Japan, 1976, pp. 109-116. 
. Simiu, E., “Revised Procedure for Estimating Along Wind Response,”’ Journal of 
the Structural Division, ASCE, Vol. 106, No. ST1, Proc. Paper 15111, Jan., 1980, 
. Soize, C., “Gust Loading Terms,’ of the 
a Structural Division, _ASCE, Vol. 104, No. ST6, Proc. - Paper 18M, , June, 1978, pp. 
8. Sun, C., Berg, G. V., and Hanson, R. D., ‘Gravity Effect on Single-Degree ~~~ ; 
_ System,” Journal of ‘the Engineering Mechanics Division, ASCE, Vol. 99, No. EM1, 
9. Tanabashi, R., and Nakamura, Tsuneyoshi, “‘An Approach to the Last Hinge Point — 
= of Tall Multi-Story Frames,”” Proceedings of 11th Symposium on Bridge and 
_ Structural Engineering (1964), External Forces and Structural Design of High-rise and — 
Long-span Structures, 1965, pp. 169-179. 
10. Tanabashi, R., and Nakamura Tsuneyoshi, ‘“‘The Minimum Weight Design of A Class — 
of Tall Multi- “Story Frames subjected to Large Lateral Forces I, Il ‘y Transactions, 
_ Architectural Institute of Japan, Vol. 118, 1965, pp. 10-18, and Vol. 119, 1966, we 
1 


. Tanabashi, R., Kaneta, K., Nakamura, Tsuneyoshi, and Ishida, S., “To the Final _ 

State of Rectangular Frames,” Proceedings of the 4th World Conference on ee 

___ Engineering (1969), Santiago, Chile, Vol. 2, A-4, 1969, pp. 179-194. Pe 
12. Tanabashi, R., Nakamura, Tsuneyoshi, and Ishida, S., “Gravity Effect on the 
; __ Catastrophic Dynamic Response of Strain-Hardening Multi- Story Frames,”’ Proceedings 
of the 5th World Conference on Earthquake Engineering (1973), Rome, Italy, Vol. 

13. Vaicaitis, R., Shinozuka, M., and Takeno, M., ‘‘Response Analysis of Tall Buildings — 

4) _ to Wind Loadings,”’ Journal of the Structural Division, ASCE, Vol. 101, No. ST3, 

Proc. Paper 11194, Mar., 1975, pp.585-600. 
14. Vellozzi, J., and Cohen, E., ‘‘Gust Response Factors,’’ Journal of the Structural 
Division, ASCE, Vol. 94, No. ST6, Proc. Paper 5980, June, 1965, pp. 1295-1313. ; 
15. Vickery, B. J., discussion of ‘‘Gust Response Factors,’ by Joseph Vellozzi, and 


«Edward Cohen, Journal of the Structural Division, ASCE, Vol. 95, No. ST3, Proc. 


“Paper 6430, Mar. 1969, pp. 494-501. 
16. Vickery, B.J., “Wind Action on Simple Yielding Structures,”’ Journal of the Engineering 
Mechanics Division, ASCE, Vol. 96, No. EM2, Proc. Paper 7204, Apr., 1970, pp. * 


17. Wyatt, T. A., and May, H. a “The Ultimate Load Behavior of Structures under 


ail 


-__eollapse gust load 9° vides a rational basis for a practical safety or ff 
7 
jon 
© 


Loading,’ 
Buildings and Structures, Tekye, gp. 
_ 18. Yokoo, Y., Nakamura, Tsuneyoshi, Ishida, S., and Nakamura Takashi, “Cyclic 
-Load- Deflection Curves of Multi-Story Strain- -Hardening Frames subjected to Dead 
and Repeated Alternating Lateral Loads,”’ Preliminary Report of the International 
Association for Bridge and Structural Engineering Symposium on Resistance and 


Ultimate Deformability of Structures Acted on by Well Defined Repeated Loads, Lisbon, | 
Portugal, 1973, pp. 81-87. 


19. Yokoo, Y., Nakamura, Tsuneyoshi, and Keii, M., ‘““‘The Minimum Weight Design 


a _ of Multi-Story Building Frames Based upon the Axial Force-Bending Moment Interac- 7 
tion Yield Condition,’ Proceedings of the International Union of Theoretical and 
4 Appiied Mechanics Symposium on Optimization in Structural Design (1973), Warsaw, 

20. Yokoo, Y. , Nakamura, Tsuneyoshi, and Denda, M., ‘‘Analysis of Dynamic Stability 

> Limits of Elastic-Plastic Structures by Incremental Perturbation Method,”’ (in Japanese), 


Summaries of Technical Papers of Annual Meetings oS the Architectural Institute of a 
Japan, 1975, pp. 693-694. 


ala: - 


Bet 


= 
| 
J My | 


 Astability analysis of shells is generally carried out by comparing the expected 
load magnitude to a theoretical critical load divided by a safety factor. This 
_ safety factor used to be much greater than that applied for other structures. 
_ The reason is that there are several secondary effects decreasing the critical 


load of shells much more than that of other structures, and these secondary 


effects have not been clarified so far. | 


the stiffness of reinforced concrete decreases due to cracks while it increases" 
. to the reinforcement. _ Thus, an additional stiffness of shell caused by a 
double layer reinforcement has been evident but it could not be taken into 


_ so much that a dimensioning procedure, for which the safety factor will conform _ 
to that of common structures, can be developed. Such a dimensioning | procedure | “f 


* account. We believe that the research of these secondary effects has advanced — 


__ will be more correct and reliable for the various destabilizing secondary effects — 

i ie can be considered at their correct places with their effective values. Puanyientnn 

Pe _ When considering a shell of perfect geometry, the differential equations of 

a  egdiikdan and compatibility restricted to small deflections lead to a linear 

7 eigenvalue problem. The lowest eigenvalue obtained is the critical force based 

on linear theory, or bifurcation force of the shell (N"",). The ‘surface load 

7 belonging to that shell is the linear critical load of the shell. However, experimental 

wr observations indicated that cylinders and spherical domes buckle at a load 
_ magnitude much lower than the bifurcation one, and the buckling deflection _ 

is no longer small. Therefore, nonlinear terms relating to nonlinear deflections 

; were introduced to the diff ferential equations. Thus, equations became ——— 

ve to describe the linear as well as the nonlinear (postcritical) behavior. ol la 

Investigations by the energy method presented the following features: 

_ 1. The load carrying capacity of some shells does not decrease, but increases 
after reaching the bifurcation state. If there is an initial imperfection of amplitude — 
= Wey the shell has no bifurcation load any more, , and its carrying capacity increases 

"Struct. Lead Engr., Planning Office, BUVATI of Budapest, Budapest, Hungary. 

s Note. —Discussion open until May 1, 1982. To extend the closing date one month, 
_ awritten request must be filed with the Manager of Technical and Professional Publications, 
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ISSN 0044-8001 /81 /0012-2381/$01.00. 


BUCKLING OF REINFORCED CONCRETE SHELLS 
| 
7 
| 


DECEMBER 1981 


= with increasing deflections. Shells of this type are not sensitive to initia 
ne imperfections. Generally those shells exhibit such a behavior at which a — 
g change in the stress distribution | CocusTing at p“". stabilizes the shell, ‘thus i 
In ‘other c cases, carrying capacity of the shell until 
fe the deflection due to buckling reaches a value, Wenes | then it increases a am 


are called the lower critical 
and the lower critical load, p'°*, respectively [this notation was used by Karman : 


_ (20)] . Shells of this type are referred to as shells of decreasing carrying capacity. 
Generally, cylinders, cones, spheres, and domes belong to this type. However, 
~ not even this value of p'°* was confirmed by experiments. F ailure loads appeared 
between the bifurcation load and the lower criticalone. 
_ Further investigations pointed out that the reason for which cylindrical and 
spherical shells buckle at a load smaller than the bifurcation was their initial 


imperfection. The maximum amplitude of imperfection will be denoted by w,. 
TABLE 1.—Ratio of Snap- Through Critical Load and Bifurcation Load 


Investigations of shells with initial imperfections were carried out as well. The 
essense of the results was that the initial deformation increases with the load 
Ep intensity. The load magnitude reaches a maximum value, after that, it decreases 


with increasing deformations. This maximum is denoted by N Lynd 7 
critical cylinders and ‘spheres buckle at the snap-through critical 

The Ratio of Snap-through and Bifurcation Load. —Ref erring to radially loaded — 
spheres and axially compressed cylinders, can be estimated by using the results — 
found in Refs. L 4, 5, 8, 19, 21, 22, 23, 28 and 34. For different values of ‘ 

w,/t, the ratio is given in Table 1. Here ¢ is the shell thickness. 
An approximate ‘formula fo for | the considered fairly good fe for w, 0, 


| 
> 
| 
i 
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cylinders and radially loaded spheres F= 6, for radially loaded long cylinders — 
F = 0.6, and for short ones, F = 1.4. For hypar shells on square ground plan, 
the value of F is zero. For hyperboloids of revolution, the value of F depends: 
2 on geometrical ratios, namely, how much the current geometry differs from 
‘that allowing an inextensional deformation. _ don 
, os Initial deformation (imperfection) consists of two parts: (1) Deformation which 
can be calculated, w, .,,.; and (2) random initial imperfection, w, .,,. The average 


in which R denotes the radius of curvature. By equating the latter two expressions, 
the random imperfection can be expressed in terms of R: w,,,,, ~ R/3,500. — 
- This value is in a good agreement with those measured on shell models (2). _ 
The value of w, .,;. can be determined on the basis of the bending theory of 


shells. A coincidence of the calculated and random maxima is rather regain Ph 


_ ‘The construction of high-speed c pore has made it possible to develop 
: programs by which complicated shell buckling phenomena can be followed very “ 
accurately. In fact, even the special characteristic of reinforced concrete shells * 
can also be described (26). However, not every design engineer has access — 

to these programs; moreover, there is very often not enough time and money 4 ; 
to perform such comprehensive computations. It is therefore desirable to provide : : 
the designers with more simple, easy-to-survey and easy-to-use methods. We > 
have tried to comply with this demand following. The stability analysis of shells 

of increasing carrying capacity can be carried out as a customary dimensioning 
procedure limited by the bifurcation load. For this reason | the interest of of this 
paper is restricted to shells of imperfection-sensitive ee ee 


_ The method of analysis presented herein is equally ‘applicable to to elliptic (positive 
Gaussian curvature), parabolic (cylinder), and hyperbolic (negative Gaussian 
shells of concrete or reinforced concrete materials. 


Material properties of the reinforced concrete needed for shell buckling analysis 

. cannot be given by single numbers (moduli), as deformations depending on 

_ cracks in n concrete, the reinforcement and the creep effect are related to the | 

- Joading i in a nonlinear manner. That is why some more detailed investigations — 9 

on mechanical properties of reinforced concrete are presented herein. 
i: The material of reinforcement is assumed to be of ideal elasto-plastic behavior. 


The stress-strain diagram of concrete [Fig. 1(a)] is curved from the start _ 
_ of loading because of the viscous deformations nonlinearly connected to stress. ‘i 


It is considered as the plasticity of the concrete. In addition a ineat creep 


having random initial imperfections can be approximated by the formula 
crit 
| 
| 
— 
ii 
| 
| | 


the ideal elasto-plastic stress- strain p ten serving as the approximation of + 
the real one, depends on whether loads have been acting for a long or short — 
while. In the case of long time loading, viscous deformations can fully develop __ 
while this is not possible when the load acts instantaneously, or for a short ' 


“time (proof load). Thus E. pand 2; for long and short time loading, respectively, — a 
are the characteristic moduli for the ideal elasto-plastic model of material. Long _ 
time loading means more than one year of action. During such a while viscous = 


deformations can almost fully develop. geten 
In the case of instantaneous loading (e. Vibration), there is no viscous 


q 


> be determined from various deformation factors calculated for the two 
— parts of the cross section (tensile reinforcement, compression zone of concrete). 


FIG. 1 Stress @)- Strain of Concrete; (0) Change of Def of Deforma- 
tion Modulus ot Concrete E. Plotted against Stress ; J 
_Pe both t the reinforcement and concrete are regarded as ideal elasto-plastic 


materials, ‘the cracked reinforced concrete cross section also has a similar 
behavior. ‘Actually, the start of yield in the cross ‘section is influenced (besides — 
the geometrical data) by the quantity of reinforcement as well. This is automatically 

_ taken into account, however, in the usual computation of the reinforced concrete a 


_ In summary, the reinforced concrete differs — the elastic homogeneous : 
material in the following ways: 


1. Viscous deformation (linear creep) of the concrete. 
2. Plastic deformation (nonlinear creep) of the concrete. 
3. Tensile cracks in the concrete (stiffness of the cross section drops” off, 
the reinforcement becomes of particular relevance), 


Modulus of Elasticity (E.. o)-—The modulus of ‘elasticity, E..,, 


_ Characteristics of the cracked cross section of reinforced concrete are modelled 7 
> 
] al 
= 


instantaneous behavior, e. Its mean value can 


h 


= 6150 / mam? (694N/m = 1p (4) 


Here, ft denotes the specified compressive strength of the concrete, that can — 
sovemed for 80% of the cube strength, thus: 1! 


In practice, the deviation of "the modulus of elasticity is taken into ioe 

by using a nominal strength value in Eq. 4. The nominal strength is a threshold : 


belonging to the desired quantile. 


q 


Viscous Deformation (Linear Creep). —The buckling deformation of a concrete 


shell depends on viscous deformations of the concrete as well. . Although, viscous a 

-_ Sr does not alter the instantaneous modulus of elasticity, its effect — 
(buckling deformation increases with time while load remains unchanged) can a 

be taken into account by reducing the value of E., (see Refs. 9,35). 


_ The phenomenon of the linear creep is commonly described by ‘the creep 


creep elastic’ 


The asymptotical maximum of the creep function is the creep factord(?,. = = aq J : | 
4 Its value depends on the thickness of the structure, the ‘quantity and quality 
of the cement applied, the water- cement ratio, and the : surrounding of the concrete : 


; (6). ror normal circumstances $. depends only on the type of concrete of 
nab 


vided of for some cases are given in Table 2. 
~Modulus of Deformation E. -—The critical load influenced by linear creep 


‘This ‘modulus is generally computed as (6,9,35): 


If loads are applied on the structure a short while (i.e., a few 
q then only a part of the creep deformation develops (6). pe? aor) frond peat 


By considering this part of the creep deformation together with the straightening | 


of the =(e) diagram, the modulus of deformation for short time loading, 4 


poe = 4) 7 = 4,730 30 VF. gies: 


N/mm’ 


Later on in E. wil always denote the value of ‘the 1 modulus of 

a corresponding to the loading time. That is, E.= E.,, for long time 

loading, E. = E.,, for short time loading, and E. instantaneous looting 


Bffeet of Plasticity.—As is well known, the buckling of centrally compressed 


a 
| 
q 
| 
| 
| a 7 was determined exactly for some structures (9,15). [he results oDtaine ere, | 
; --_-—s confirm that the approximate consideration of the creep effect, by using a reduced — | 
— 
| 
| 
| 
| 
built OF nonlinearly elasic or plasic material Can be investigated Oy USINE | 
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theory I. of substituting the nonlinear stress-strain 


_ its tangent at the point in question. Shells built of materials of plastic properties — 
buckle in central compression similarly to straight bars (3). Thus, the theory 

4 of Engesser can also be applied for shells. In the stability analysis of centrally 
compressed structures, and for concretes of good quality used for shells, a 

ee for the elastic modulus E varying w with the stress is available in the } the 


The value of E is depicted by the curve 1 in Fig. 1b). For the ‘computation 7 
of -deformatio ons to eccentric instead of the elastic 


Newtons per square millimeters (pounds per square | inch) 
(0.69 x 10°) 


FIG. 2.—Stresses in Ecc 


4 the n modulus of deformation E. is s used. This value cuts off the p ane line & 
in the figure by a dotted line. _ 


approximating preceding condition. "line 2 in Fig. 1(b)}. In this formula 
§ means the axial membrane stress in elastic buckling. Rewriting the — 


7 
TABLE 2—Values of 6, Creep Factor, 
= 
— 


(pe 
=(p.) 


in which p®™ = critical load taking plastic features also into account; ge 
 snap-through critical load assuming linearly elastic material; and p,, = the strength 
_ of the shell if buckling didn’t occur. The same expressions were me obtained 
for shells of ideal elasto-plastic material. Thus, the critical 
load based on elastic theory 1 must first be obtained. 


STIFFNESS: FOR Concrete ano Reinrorcep Concrete Cross A 


Generally, the real dimensions of concrete and reinforced concrete structures 


- differ from those planned values. A long series of measurements carried out 
in the Hungarian Institute of Building Research (25 and 33) pointed out that 
the thickness of plates and shells are, in an average, 1.05 times bigger and 

are, at worst, by 10 mm (0.4 in.) less than their planned values. In buckling, 

the average value of the thickness is the relevant, while the influence of the 

weakest cross sections is to be felt only at rupture. Thus, in the stability ol 


of shells, similarly tc to » that of other structures s of reinforced concrete, ‘cfoss- sec- 


—-~xHere, B and j Dmean the flexural and extensional stiffness of the sel, respectively. 


7 = two stiffnesses preventing the buckling of shell, are combined in the 
‘shell K. For the homogene shell plate is K= V BD = Vv 
«Stil Stiffness data are computed on the basis of State Il. of the concrete 

ideal 
? a In cross sections of concrete ont reinforced concrete, the flexural and 
sagem stiffness depend on the reinforcement and on the eccentricity of 
the axial compression force, because the cross section cracks: more and more 
® Consider the cross sections, in Fig. 2. The upper one shows the stress diagram — 
of an eccentrically compressed, uncracked cross section of linearly elastic 
concrete. In this case, the flexural and extensional oe : 


neglecting the contraction factor (Poisson’s ratio), are: 


col 


| 
| 
| 
| 

g 

x at properties Mav oe taken account With thelr planned 
| i e given according to the explanation given in Eq. - 

q 
| 

| 

| 


which J = moment of inertia per unit width; and A 4 
Unit width. Thus, the shell stiffness of the uncracked cross section of concret - 

iss 


, also reinforcement is applied, the stiffness data tasulenn by taking into 


2G 


account the area of the steel with a factor n= E, /E,. .. Here, E, denotes the 
elasticity modulus of the steele 4 

TABLE 3.—Values of andy, 
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, ‘uncracked reintorced concrete 


suncracked concrete 


with double mesh 
ot 


~ Thus: 


~ conc ~ reinf.conc 


The shell stiffness of the uncracked cross section of reinforced concrete is: 


reinf. 


i 
) 
| 
‘ os =| 
| 
| 
1 - Her® lower index “‘reinf.conc’’ means that n times the area of the steel =f 


=—  [l+ - any 


: 


in which + + 3m (1 — 7 

a my which n = E,/E, (the ratio of deformation moduli for steel and concrete); — = 
= h’/t (the ratio of the « covering layer of reinforcement and the thickness ns 

of ‘the cross section); and » = the ratio of the area of steel and concrete. __ Se 
_ Single mesh or double mesh of reinforcement used to be applied to reinforced 7 

“concrete shells. In the case mesh, h, 5, double 


7 


coe 


teint conc. pe 
eft 


AIG. 4.—(a) Change Stiffnesses of Reinforced Concrete Cross Section 


this value is about = 0.2. With the use nen 22 we the 
Vo, of single mesh, and the values double mesh, for various steel 


| These values are shown in Table 3, and in Fig. 3 as well. a. 
When the cross section cracks, its stiffness decreases. The diminution on depends 
on the eccentricity of the compression force. 
_ Diminution of stiffness of a reinforced concrete cross section, for which — 
i = 0.20, is depicted in Fig. 4(a). When the eccentricity gets close to a “‘half 
cross section’’ magnitude, then stiffness become almost identical with those 


reinf.conc reinf.conc 


reinf .conc E. x A reinf. conc Et(é + r, mp). . (24) 


in which ,, = the factor of tension stiffening (see Appendix I). The stiffnesses, 
and differ from BY" and respectively, becau q 
the stiffness due to the uncracked concrete in tension gets lost as the concrete 


cracks. Thus, the shell stiff fness ben to iy eccentricities can be en 
as: 


By applying» the notations a/h =§ A,/bh and we restrict ourselves 


‘the commonly used symmetrical reinforcement, the values of and 


rein .conc 


: sl crack _. can be expressed with the use of the theory (13) of reinforced concrete 
Substituting these into Eq. 25, K (pure t bending) can be obtained. 


reinf conc 


"This value may be used in the cases _/t > 0,5. 
Introduce now the expression 


reinf.conc 
and let’ 's express the value from Eq. 26: 


+ §) (1 + 2&? + 2m? — 2€ — 


_ The relative height of ae compression zone in pure bending i is obtained from 
4 that condition, the statistical moment is t to ue am zero. From this condition: i 
computed values of , 2 Teferring ton = = 0.5 andy i 
shown in Table 3, and ‘also i in Fig. 3. oT et ial 
interpolation in the Table is allowed. Index 
(e, /t) = 0, while index o refers to the case (e, /t) >0. and 
however, to single mi mesh and double mesh of 
al there is no reinforcement in the cross section, tension force disappears 
when the tension zone of the cross section cracks. Thus, bending moment can | 
be balanced only by compression forces. For concrete cross sections, conse- 
- quently, the resultant compression force must be inside the cross section. Due | 
to shrinkage, cracks may also appear in concretes. Thus, for being on the safe 


side, the tension | strength of t the cross section is to be neglected even for the 


_ The compression area of a cross section of unit width is according to the __ 


D 


f 
~ 
| 
q 
| 


4 “centroid of the cross section. Sy~> shell buckling, we assume e = cw, c= 1.0 


For cross sections without 


‘second formula), we arrive at the value of 


ons 


1s, for the ca cases 0<e/t<0. the value of shell is as: 


conc conc ac 


Sites sites computed with the use of Eq. 34 are depicted in in Fig. 4b). 


For Sueus or Exastic, Concrete on REINFORCED 


Inside one buckle 1 wave vave the values of ‘the deflection and thus that of the -¥ 

"shell stiffness vary. However, it is mainly the shell stiffness at the maximum Pe 
deflection that influences the buckling deformation. The values of shell stiffness a | 
elsewhere are greater than that at the maximum deflection, thus, we are on 
. - the safe side when supposing that this minimum shell stiffness is valid throughout | 
the buckle wave area. This ensures that the differential equations are of constant 
‘This. procedure is demonstrated in Fig. 5 on the example of a cylindrical — 
shell panel for the “homogeneous, elastic” solution is to be found in 
The buckling atin for the perfect shell is depicted by a parabolic curve 
starting from the limit point of the bifurcation critical load, in Fig. 5. The 
minimum value of this curve is the lower critical load. When an initial deflection ; 
magnitude w, = 0.25 ¢ is present on the shell, the load-deflection curve is 
.. the lower part of Fig. 5 the shell shell stiffness J K of concrete and reinforced — 
concrete cross sections is shown for varying values of eccentricity e= w. The 


be reduced proportionally to the - diminution of the shell stiffness K (e). ast 
the ene concrete and concrete cross sections, respectively, the aa 


j 
a. 
| 
| 


cont continuous and dotted lines, in the upper part of Fig. 5, are 
_ maximum values of these curves give the critical loads based on elastic theory — 
belonging to the considered eccentricity, | for reinforced concrete and concrete 


When considering the por stiffness from start, we. arrive at the 


| 


—, 
eth | 


concrete 


| 


quner. 


concrete 


5. —Determination of for and Reinforced Concrete Shells 


G. —Comperteon of and nd Approximate | Values of p...,, for Concrete Shells. 


urve drawn by the thin dotted line, in the upper figure. As it is seen , the 

critical load so obtained is below that computed before. That is, we are again 


n the safe side when doing so. Goria 
_ By using this method, we computed the snap- -through critical loads of i isotropic pt 


concrete shells for different values of w,/t, and we itn 
very much on the values and 


7 wee 


— d 
q 
q 
— 
— 


e.g., the critical load of an an isotropic concrete shell (without tension 


low 


Strength) for which w,,,, = 2 and = 0.25 p.. The resulting curve computed 
the aforementioned method can be : in Big. 6 drawn by the continuous 


cracks in concrete appear at discrete distances, the: shell stiffness 


approximated with the following formula: 


a i in case of e, /t > 0.5 4 can be calculated by multiplying that of the homogeneous 
shell by the factor Thus 7 
oe 

7 For the cases in which e, /t <0. 5, the values of conc are between those 4 


Of Poonc ANd Prom: On the base of ‘more investigations, the 


Preint =, cone +. ron hom 
2 


“Supposing that 


-Preinf.conc ( =1 


the following manner: 
Values of B are With the use of of | B, 
Snap-through critical load of reinforced concrete shells can also be well 
THROUGH Crrmicat Load FoR REINFORCED ETE OF Puasnc ~ 


. Thus, when we finally consider all the effects (i.e., 


a 
“= 
‘Thus, Tinally, the shap-througnh critical load OF an elastic reinforced concrete 
shell can be computed as f 
a 
; _ Plastic behavior of reinforced concrete shells can be considered by using . 
the factor given in Eq. 11 
— 


2304 
"cracks, plasticity) described i in the foregoing, we arrive 


pe at the formula of the critical as is given in the following equation; 


In the formula, denotes the snap- -through critical load 


of a reinforced concrete shell, plastic behavior included. 


reint 


iret, 


0.2 


FIG. 7. —Velues of B Plotted against and e, /t 


Haas -v Koten 


Griegs 
Haas — v. Koten 


ich 


‘i In an earlier paper (11), we computed buckling loads for the copednaunal ia 


shells that had been published in Refs. 15, 17, 32, 31, 7, 16, 29, 14 and 27. 


_ We used Eq. 41 and made comparison between the results Its of « our ur computation | J : 


| 


| 
| 
| 
| 
| 
| 
y 


‘Siok ok 


and the experimental ones. These vahees plotted against ast R /t are seen, in Fig. 


while the standard deviation is 20%. This shows that our procedure presented 


The procedure proposed in this paper gives a probable mean value for the | 
_ buckling (failure) load. Uncertainties must be taken into consideration by applying 
safety factor. According to ‘the probability theory ‘the safety factor bs 


d foll 


= lowest possible value of the carrying. capacity be higher than the greatest 


domain of mode! Jena, Germany 5 ~ 
2. Jena,Germany 
ty 

AA, Hilling, USA 
pt 6. Hamburg, German 
7. Windward, PB ‘J r 


Wales, GB 


Atbuquerque, USA” 


= 10. Belgrade, Yugoslavia 
Beigrade, Yugoslavia Vv 
12. Algeciras, Spain 
Novosibirsk, su 


= 
| 15. Hungary 


T shells 17. Puerto Rico, USA 
18. Cleadon,GB 
FIG. 9. Comparison of Effective | Loads and Critical Loads of Several Erected Domes 
possible value of the load intensity. In those cases, where the structural model — 
used for computations is in fairly good agreement with the real structure, _k 
. i: is to be set to 1.25. Otherwise k, must be somewhat higher. ie 
—. _ The value of a depends on what safety levels are required. If “considering 
. a Gaussian distribution and the usual 2.5% value for the quantile, we obtain — 
£ denotes the resultant relative (variance) in in Eq. 42, 42, by 


pi + + 


— and Se the dared load, respectively, are under the square root. The standard 
deviation of the load values used to be about 10%, thus, Sioa = 0.1 0. 10. — 


| 
| 
4 
q H 
|_| 
| 
= 
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_ For the cylinder and dome experiments investigated i in Eq. 43, we found a 


deviation of 20%. As the number of cases investigated was rather low, we con- 
sider here a deviation of 1.1 x 20 = 22%, instead. With common reinforced con- _ 
crete constructions the deviations of the no-buckling eer. isabout 10%, thus — -~ 


we set = 0.1 and then S., can be expressed as (1.1 x 0. 0.2)’ - 0. 


a ‘Finally, the resultant deviation is: S,= V 0. 17(3) = 0.25. Thus, from Eq. 42 — 


ke 
— 2(0.2 


Consequently, the value of the safety factor must be 5 at 


‘it should be greater if there are any more uncertainties. 


_ We computed critical loads for several erected domes of great dimensions, 
and depicted these in Fig. 9. These data show that the use of safcty factors 


a For ‘Structures with a constant post- buckling load- bearing capacity are S., 
60.1, 8S, = 0. 17, and k = = 1.90. In the case of an increasing post- buckling 


 load- bearing capacity are S., 0, S, = 0.14, andk=1.75. 


 Let’s Griggs ( (15) spherical dome models, made of concrete an 


“rely built-in around their support. Griggs used the formula 


(840 ote. 
2, 134 


17 x 
28.8 x 10° KN/m 


to compute the critical load of his shell. : 
On plastic models just so formed as the aforementioned ones, experimental - 
- “results c = 0.57 and c = 0.63 were obtained. It conforms to the theoretical results 


+ = Ve: 0597 + 0.0677 = a 


As there was no reinforcement in the concrete coin 


_ and and fr from from Eq. 35, Pcon 


V3 > (7 16 x 10° x 0.25 ait 

= = 32.2 psi(222kN/m’) . 

= Pconc X Pern = 0.55 X 32.2 = 17.9 psi (123 KN/m? 
(1-28 17 x 10° 
Po 


= 0.55 . 


l 


and = 0.5 x 0.73 =0 


= results of the nonlinear sh shell buc buckling has been e to 
a to the reinforced concrete shells. A simple relation i is derived for the snap-through - 
¢ritical load of reinforced concrete shell. The solution contains the special 


“properties | of reinforced concrete, i.e., cracks, plasticity, creep, and quantity 


The main results of the research are are the 
1. The shell | buckling stiffness was seen to be is oe Le of two a the 
‘ie and the axial stiffness of the reinforced concrete cross section. inal 
_ 2. The critical load is reduced by the plasticity and the imperfections. __ 
wit The critical load of the reinforced concrete shells decreases by the cracks 
of concrete and and i increases ses by th the 1 reinforcement. “UT x BELT) 


Fromt 
| 
| 
| 
| 
| 
| | 
 &§ 
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The pe parameter, \,,, of tension stiffening. 
‘The parameter is as the ratio the } steel elongation, 


; 
‘steel elongation is €,= x /E, (neglecting of the effect of concrete 
in tension). The mean elongation of the reinforcement €. Situated in the ef oranda 


in which x, = the stress in in ‘the failure 
= =f, ye 2m = the stress in the steel reinforcement calculated in the cracked 


= 


0. 

= 20 N/mm’ (3 38 x Psi; f, N/mm’ 2. 76 10° psi); and 

immediately after "cracking; K, the coefficient which characterizes _ 

the bond properties of the bars (K, = 0.4 for high bond bars, K, = 0.8 for 


ahs bars); and f, = the yield strength of seentoccemnent. 


By applying the dates given in this this paper, , we “obtain: for both single mesh 
and double mesh r reinforcement 


On the basis of this relation the following formula can be written: oolu le oh 


Some of characterizing tension stiffening are given in Table 4 
= 20 N/mm’ (1.38 x 10° 400 N/mm’ (2.76 x 10° psi). 


| 
4 
| 
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modulus for concrete, at start of loading; 


moduli of deformation belonging to short time - loading and 
longtime loading, resp; 
compressive strength of concrete; .T wilo# 
compressive stre strength of concrete cube; godin 


wale 


= Jength of buckling waves in x and y directions; ions; 
= ratio of the moduli of steel and concrete, 


E 


a 
| 

| 

Kay BD 
| 


ent = 
we = lower critical force | and load of shell 
computed by nonlinear theory (with large deflections); 
| Snap-through critical force and load of the imperfect shell, 
v.e. maximal carrying capacity computed by nonlinear theory; 
snap-through critical load of reinforced concrete or concrete 


shell of linearly elastic behavior; 


reinf.cone = SNap-through critical load of reinforced concrete or concrete 


shell strength 


shell radius; ‘tn 
= maximum amplitude of buckling « detection; — 


calc 


factors used in the computation of shell (index 1 and s2 re wane 


factor; 
Pp 
factor of tension- stiffening; 
= ‘ratio > of ¢ cross sections of steel and concrete; and ‘ie 
factors used in computation of cies critical oad of 


“eg 
= _ random and calculated s of deformation: 
| 
| 
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SHEAR REINFORCEMENT FOR SLABS* 


By Walter H. Dilger’ and Amin Ghali,’ Members, ASCE 


Concrete flat plates are subjected to lntge bending moment and shearing force 


2 at their connections with columns. Their combined effect can cause failure 

In flexure, reinforced concrete slabs exhibit a great deal of ductility; extensive 
a deformations occur before their ultimate strength is reached. Design codes place 
ae reliance on this ductile behavior which enables slab systems to 


redistribute moments. Complete redistribution of f bending moments ¢ can an generally ‘é ; 


attention be given to both strength and ductility. 
Perhaps the greatest source of problems for slabs in practice, apart ton 
excessive deflection, is shear. . The so-called punching shear failure ‘of slabs — 
without shear reinforcement occurs suddenly without warning so that, normally, 
no remedial measures can be applied prior to failure. Several such failures 
have been reported in the literature (9,6). 
_ The development of the concrete design codes in North America in the last 
_ two decades and the extensive use of prestressed concrete in slab construction 
have resulted in the design of of thinner and thinner slabs. Particularly if we use 
prestressed concrete, we can eliminate problems related to deflections which 
may govern the design of reinforced slabs and, therefore, the main criterion 
‘for the design is ‘punching shear.’’ In addition, the architects may require 
openings for ducts in the vicinity of the columns which may substantially weaken | 
the punching shear resistance of the slab around the column. insiofiva joa ai 
Tf the design of a slab is found adequate for flexure, but | not for punching — 
_ shear, and an increase of the depth or use of higher concrete strength in the | 
whole slab are not desired, three possibilities exist to increase the shear strength: A 
f 1) Increase the slab thickness in the vicinity of the column (provide a drop | , 
“Presented at the April 14-18, 1980, ASCE Convention and Exposition, held at Portland, 
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achieved in a SiaD System prior tO lallure provide at shear ' 
prevented It ic therefore _imnortant that in desion hen considering nunching aah, 
 &§ 
— 


~ panei); Q if this is not ani flare the top of the ———- (i. e., use a column 
capital); and (3) provide shear reinforcement. 

_ In many cases drop panels and column capitals a are For example, 
‘in apartment bulidings where the bottom surface of the concrete slab is hol * 


ceiling of the rooms, a drop panel would be unfitting for esthetic reasons. 
a capitals are also not suitable for many types of buildings. In addition, — 
they “complicate the formwork and the construction process. In recent years, a 
the so-called ‘‘flying forms’? have been introduced which employ formwork 
panels of approximately the same size as the floor panels and are moved in 
their entirety from one floor to the next. If these large panel forms cannot 


be used for errr the slab or if they have to be modified to accommodate 


Dirrerent Types of SHEAR REINFORCEMENT - ve 


- for sabs have been many attempts to develop effective shear reinforcement 
a for slabs. Tt was, » however, recognized | that shear reinforcement for slabs creates 


types s of shear reinforcement have been investigated, “most of which are are sum- 
marized in an excellent paper by Hawkins (8), who gives a comprehensive 
list of references. The different types may be grouped into four different | ; 
categories: (1) Bent-up bars; (2) stirrups” (closed and U- |-shaped); (3) sheathead 
_ The first type of shear reinforcement was first ieventenied by Graf, and 
later by Elstner and Hognestad, Franz, Anderson, Rosenthal and Yitzhaki. Closed — 
_ and open stirrups were tested by Franz, Wantur, Carpenter, Keefe and Narasin- 
ham. Shearhead cages were used by Moe, Anderson and Leonhardt. Shearheads 
consisting of structural steel plates or structural steel sections, or both, — 
_ investigated by Tasker and Wyatt, Corley, and Hawkins. A few typical examples 
‘ are shown in Figs. | to 4 (exstensive lists of references are given in 8 and 9). 
_ In many of these tests it is found that the shear reinforcement consisting 
of bars is not fully effective, which means that it does not reach the yield 
strength at failure. The reason is that the anchorage of the bent bars or stirrup = 
[2 is not sufficient to develop the full capacity of the bars. In case of bent-up 
bars, the tensile force in the inclined leg is associated at the bends with compressive a 
: stresses higher than the concrete strength and with a pull-out effect of the 
horizontal portion of the bar. This leads to an increase in width of the shear 4 
_ crack, which when extending to the compression zone results in the destruction 
of the concrete shear compression resistance before the shear reinforcement 
reaches its yield capacity. Similar observations have been made on ‘stirrups. 
A 90° bend at the end of a vertical stirrup generates high concrete compressive 
- stresses under the bend and also a certain pull-out effect of the horizontal P| 


leg under high stress in the vertical leg. A vertical slip of more than 0.04 in. 


(1.0 mm) has been measured in beam shear tests. In beams, this is normally 7 
Rot critical, but i in thin slabs the slip explains why it it is difficult or, Perhaps, 


care for anchorage details for shear reinforcement i in n slabs has been emphasized ns 


many researchers. 
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—— Anchorage slip in shear reinforcement made up of welded 5 wire febtic was 
measured by Seible, et al. (12). To monitor the vertical movement of the shear 
~ teinforcement, small studs were tack welded to the bars and their vertical 
_ movement relative to the surrounding concrete was measured by means of dial 
on gages. A typical set of results is shown in Fig. 5. In order to reduce this pull-out 
| _ effect, Leonhardt (11) pooposed that the vertical bars be welded to ) rings of 
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er 


round bars as shown in Fig. 3. But to thread ad all the flexural steel through — 
"one or more such rings is a difficult task at the construction site. beisanaie adnan 

7 s In order to make normal stirrups reasonably effective, they have to enclose 
the top and bottom flexural reinforcement. For this reason, the installation 
= of these stirrups in situ is a time-consuming job and in thin slabs with 6-8 


(150-200 mm) thickness it is virtually i impossible. bod 


i 
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In view of the difficulties, type of shear 
was developed in the late sixties: The so-called shear head which resulted from 
a research of Corley and Hawkins (4) and is not part of the design codes 
in North America. In Europe, a closely related development is the Geilinger 
collar, shown in Fig. 4. It consists of steel channel sections forming a rectangle 
somewhat larger than the column ‘size and the load is the 
FLEXURAL 
_NOT 


SECTION oat 


3.—Vertical Bars Welded to Horizontal F Rings 


‘The. 


slab to the collar and from there to the column by iemeens steel 1 plates welded b 
to the channels. Because the bearing stress under these plates would be excessive 7 
an additional horizontal steel plate i is required 


‘expensive. Second, "they may interfere with ‘the: reinforcement of the column = 


_ and may be difficult to place. Third, the depth of the structural steel sections _ 


| 
| 


must be substantially smaller than the slab thickness as there are soutien two 4 


layers of reinforcing bars to be placed above and below the structural steel — 

_ sections. When the concrete cover is 3/4 in. (19 mm) and the reinforcing bar 7 
diameter is 5/8 in. (16 mm), the height available to accommodate the structural » 
steel sections is equal to t — 4 in. (= t — 102 mm), in which t = the slab thickness. 

In relatively thin slabs, the available height is not sufficient for the required 

4 structural steel section and shear heads cannot be used. This difficulty may 7 
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_*FIG. 5.—Anchorage Slip Versus Applied Load at at Top Ben Bent of Shear Reinforcement — 
equal to the slab thickness ai and butt-welding to it the flexural reinforcement 
interrupted by the collar. fi quia ye 
_ Another technique to increase . the shear strength of reinforced concrete flat 
"plates is by prestressing the zone around the column by means of high-tensile 
bolts passing through vertical holes in the slab tightened manually and anchored 
to steel plates at the top and bottom faces of the slab. This method was aees 
in practice to strengthen existing slabs which showed some sign of shear distress. _ 


Ghali, (7) demonstrated experimentally that substantial i increase in 


| 
| 
in ny 4 
| 
q | 
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i 
shear strength can be achieved t by ‘this | technique. However, in a normal case -_ 
_ it may not be architecturally acceptable to have the anchor plates and nuts 
"protruding outside the slab thickness. Also, application of the prestress adds oF 
_ The difficulties considered in context with the different types of shear 
im reinforcements lead to the following requirements for shear ae Gone in 
. Good anchorage at top and bottom of shear reinforcement. 
. No interference with placing the flexural reinforcement. 


. Installation in thin slab should be possible. 


. Anchor plates should not t project a above slab surface. a 
3 It must be economical. 
_ The type of shear reinforcement w which ‘satisfies all tl these requirements consists. 
5 of individual shear elements with small anchor plates at top and bottom. The — 


placing i in 1 the | forms easy, are arranged around the column with 1 no interference 
the flexural reinforcement. 

¥ 
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Types of Shear ELements 


of Calgary (1, 5,10, 12): = 
1. oe 2 
7 2. Nelson studs welded into steel strips (Fig. 6(b)). 
Shear elements with round and square plates at and Fie. 


4. Bent welded wire = 


Of these, the wire fabric did not behave satisfactorily 
ina specific way (12) which would complicate construction and will not, therefore, ; 
‘befurtherexamined here, 
_ Nelson studs with regular size heads did not turn out to be fully satisfactor 
because of the relatively small size heads (head diameter equal twice the bar 
diameter). Anchorage slip was observed and only a fraction of the yield strength © 
of the studs would develop. This result was not surprising considering the 
extremely high bearing stresses that develop underneath the anchor heads under 
-. Excellent results were obtained with shear elements (experiments are presently 
underway at the University of Calgary to determine the minimum required — 
q of the anchor head) for which the anchor head had a bearing area 9 or 10 
times the stem area, for yield strength of the stem = 43 ksi (450 MPa) —_ 
concrete strength = 5,600 ksi (39 MPa). 
_ The segments cut from standard I- beams are also well suited for the purpose a 
because the thickness of the flange increases towards the stem. However, lll 
a production point of view there seems to be no inexpensive way to cut reve 


| 
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FIG. 6. _— oa of Shear Reinforcement Elements : and Top Views Showing their 


small segments from a rolled I-beam. Practical solutions seem to be Nelson- “type q 

studs with a heads or vertical steel bars with square or ro round anchor pm. | 


+ 
To date, a cae of 40 slab column connections with new shear edntecvemste 
shave been Of these tests, six were Subjected to unbalanced bending 
hes the results are 
‘not fully etachesive. All other test s specimens \ were » subjected to axial force 
only. The flexural reinforcement ratio was approx 1.1% and the design concrete 


Most of the test specimens had a 6-ft (1.83-m) square slab, 6 in. (153 mm) 


thick, and a 12-in. (305- mm) square column. Four slabs had rectangular columns 
8 X 24 in.’ * (203 x 610 mm”) and 11 slabs had a 9.8-in. (250-mm) square column; 
the latter size was chosen in order to reach higher shear stresses. 
Perhaps the most important conclusions of these tests are that, with the good 
anchorage provided by the anchor plates, the full yield strength of the new 
_ types of shear reinforcement can be reached, and that only the shear in excess 
of 4 V St! (psi) (= 0.33 Vf. (MPa)) needs to be resisted by shear reinforcement. 
To back up these conclusions, a few characteristic results are listed in Table 


CONCRETE SLABS 

| 
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Ee 
TABLE —Gummery of 


in inches square 
s (in 
milli- 


meters) meters) 
(3) 


65.96 


= 65.96 
(1,676) | 


67.00 
= 702) 


(315) 


“b, = perimeter of a critical section at d/2 from column face. 


d0.33 Vf. (MPa)= 


_*A, = cross section area of shear reinforcement elements on the first row at a. ; 


1. In this table the: yield ca capacity of the : shear elements located within a a distance 


d from the face of the column_is added to the shear capacity of concrete, — 


ultimate load V,,.,, (d = the effective depth, b, = the perimeter of a critical ca 


section at a distance d/2 from column faces). ‘Only those tests in which the | 
‘shear elements within a distance d have reached yield are listed. The American — : 
Concrete Institute (ACI) Code 318-1977 (2) and the Canadian Code CAN | 
(A23. 3- idl (3), require that the shear reinforcement ‘carries the shear stress 


in excess. of 2 Vf (psi) (0.17 Vf! (MPa)) and that t the ‘total ultimate nominal | 
shear_stress by the concrete and steel does not exceed 6V f’.(psi) (0.5 
Pa In the graphs Fig. 7(a) and (6), the applied load V is partitioned into V, ..., 
Vern resisted by shear reinforcement and concrete, respectively. The part 
V ies, (Ordinate below the: curve) is the tensile force measured in the shear _ 
reinforcement located within distance d from the face of the column. The part 
of shear resisted by concrete is equal to V..... = V — V,,.,,. The fact that | 


the concrete its capacity to resist a shear stress of v.= Vi. (psi) 


=4 (psi) (psi) b,d = 0.33 Vf’. (MPa) b, dand with the experimental 


inch 
| (in mega- 
4 | @ _ 
MV6 (14) 445 66.1 | 0.46 
29.0) | (1,452) (113) | (298) | (325) 
4,680 S75 | 70.8 0.92 50.8 
490 | | | om | 
33.8) | (1,400) | | G22) | G94 | GOO 


CONCRETE SLABS 


in in _ in kips Description 
_| newtons reinforcement 
| | -segments cut 
(632) } from standard I- 


I-segments cut pn) 
standard 
beams [Fig. 6(6)] 
7 62 | 0.94  |I-segments cut 
(623) from standard I 
(1.28 Bars with square 
anchor plates 
[Fig. 6(c)] 


Bars with 


column face (spacing between rows = = 


= 0.33 33 VF (MPa) in combination with effectively anchored shear reinforced 
s here evident by noting that at ultimate, the ordinate Vorest is smaller than 
the value (V-—4V (psi) 6, d), represented by line AB in Fig. 7(a) and (b). 

_ The c current codes (2, 3) a assume for calculation of cross section area 


than 4 V f’ (psi), in the sections to follow a more conservative value is suggested 
the design of shear reinforcement (Eq. 3 or Eq. 4). 


_TOWS s of shear elements have yielded before the ultimate load i is ; reached. a 
_ Stress in the shear elements is indicated in Fig. 8 at three stages of loading 

for the slab reinforced by shear elements with square anchor plates. It should é 

be emphasized, however, that this test is perhaps somewhat extreme with regard 


$7.9 148.9 j 
400.0 | 1124 127.0 70 | 1.13 |Bars with square 
q Ri; 24 | | | 1.13 Bars with square 
A 
Vf, (MPa)}, d). Although the test results represented in Fig. 7(a) and 
— 
father fine distnoulion OF shear ciements, nonetheiess emonstrate 
a that well distributed and well anchored shear elements generate a ductile type 
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FIG. 7.—Partitioning of External Applied Force V into Inte Internal Forces Vs 0s and dV., 
Resisted by Shear Reintorcement and 


Uf 
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FIG. 8.—Quarter Top Views of Test Slab Indicating Stresses ir in Shear 
Elements as Fractions of Yield Stress at Load Levels: 90%, 95% and 100% of Ultimate | 
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“CONCRETE SLABS 
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boganr SLAB DEFLECTION, in(mm) 
anion FIG. 9.—Load Versus Deflection 


(MPa) 
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‘SHEER STREES, v, 


a 


be ob PROPORED FOR 


BETWEEN COLUMN FACE AND CRITICAL SECTION 


BY EFFECTIVE DEPTH 


Prevesti 


«AIG. 10. Proposed Design Values for Nominal Shear Stress that 


_ can be Resisted by Concrete at Critical Section outside Shear Reinforced Zone 
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hike stress, Yay "reached in the experiments ia Table 1 ‘is 
% 2 VS. (psi) (0. 59 VSe (MPa)}. In another test (to be reported ina separate 
paper), shear stress, as high as 8° (psi) [= 0.67 Vf". (MPa)] was reached. 


These results confirm the statement by Howkins (6) that the —— ACI code 


that V, cé can high 10 0. 83 (MPa) (MPa) ‘Tt be 

_ noted that in all cases the ultimate load, V,,, cannot exceed the value corresponding ~ 
— to the flexural strength of the slab which can be predicted by Johansen’s yield-line 
theory. In the authors’ tests previously mentioned the value of v, = 8 Vv f’. (psi) — 
[0.67 Vf! (MPa) } closely corresponds to the flexural ‘capacity | of "the specimens. _ 

_ A number of test specimens failed by punching outside the zone reinforced © 
for shear. Plotting the concrete shear stress v’ at failure in the critical section 7 
(at distance d/2 away from the last row of shear elements) as a function of 
the distance of the critical section from the face of the column, it is found 


that decreases with the distance from the face of column (Fig. 10). 


DistrisuTion of SHEAR REINFORCEMENT 


In earlier tests (10) the shear elements w were in a 
around the column. It is found, however, that this is not the ideal arrangement = 
because of some interference with the flexural steel. Placing the shear element 
parallel to the faces of the column is more oe _ and such arrangements <_< 


05! 


FIG. 11.- —Quarter Top View of Test Slab Indicating Stresses in Shear a ; 


Column. Load Level: 96% of Ultimate 
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i of shear failure. Fig. 9 indicates that even after ex 
_ the connection is still capable of transf 
that there is a redistribution of interna 


The question, whether a concentration of the in bends 
_ facing the column sides (Fig. 6) or uniform distribution in rows covering the 


entire area around the column (Fig. 8) is more effective, is dealt with in a 
: ; recent test series. Although the experiments are not totally conclusive the failure 


is outside the zone reinforced for shear, it is found that the reinforcement 
within a distance d from the face of the column is ‘Stressed to the same extent 
7 in both cases. There is, however , one ‘significant difference: The transverse 
expansion along the ‘diagonals of the slab with banded shear reinforcement 
3 ‘substantially higher ‘than in the slabs with uniformly distributed shear rein- , 
forcement. Accordingly, there is a tendency for the shear elements at the edges . 
— of each band to be stressed higher than those inside the band (Fig. 11). ges 
the overall behavior is ‘Rot s significantly | affected | by concentrating the ‘Shear 
reinforcement i in bands. Hain wodiw ‘vd od 
- T he distance of the first row of shear elements from the face of the column 
must not be excessive. In one experiment, the first row of shear elements a 
was placed at 3.2 in. (~ 0.70d) from the column and this resulted in a punching» 
shear failure e along | the face of the column before yie yielding of the shear elements. PS 


It is interesting to note that the relatively flat failure crack (20°-30°) rome 

observed in tests without ‘Shear ‘reinforcement, or outsice the shear- reinforced 


"FIG. 12.—Flat Crack (25°) outside Shear Reinforced Zone and Steep Crack (50°) in 
F zone, does not occur in the zones with shear reinforcement as long as the | 

_ shear reinforcement does not yield. Only relatively steep cracks (40°-50°), mostly 
bent-over flexural cracks can be observed prior to yielding of the shear — 

reinforcement (Fig. 12). A shear crack thus intercepts the shear reinforcement 

_ Situated within a distance approximately equal to d. Thus, in a design equation, 
= 
the shear force. V, resisted by steel in any one ‘critical section should only | 
account for an area of shear reinforcement located within a width = +d/ 2 


_ Shear reinforcement for flat plates must be securely anchored to the concrete | 
as close as possible to the top and bottom slab surfaces. The following are 

_ Suggested rules for the design of such reinforcement. The symbols used and 


the equations adopted comply as much as possible with current ACI and Canadian | 
Codes (2, 3). — 
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—should be based 


in which vy = the nominal shear stress due to factored load; ¢ = 0.85 = 5 : 


the strength reduction factor (2,3); and the nominal shear strength v,=8V Fo (psi) 
2. When shear reinforcement i is required, it shall be at arranged i in rows parallel 
to the column face and shall cover a distance away from the column face 
such that nominal shear stress, v,, in the critical section at d/2 from the outermost — 
Tow does not exceed > v/; in which v’ = the nominal shear stress that can 


be resisted by concrete at ultimate without shear reinforcement. The value of a 7 


_y! can be taken from the graph in Fig. 10 or calculated somal met otateiere 


or wi=0.17)1+— 
in which a = ratio of the distance between column face and critical section 
to the slab effective depth d. The values of a and B. in Eqs. 2 and 3 must 
be equal to Or greater than | and 2, respectively. Thus, at sections closer than 
d from the column face, use a = 1; similarly, when the column aspect ratio 


is smaller than 2, use B. = 2. 2 re 
(3 The shear reinforcement shall be sufficient to resist the nominal shear 


stress inexcessof 


(4) 


4, Thecross section ares area shear per ‘unit tlength of the perimeter 
of the pied section 


l 


an in which s = the spacing — rows meseunil perpendicular to the column 


i face; and f, = = ‘Specified yield ‘Strength o of the : shear reinforcement. . If the anchor 


_ shear reinforcement cannot develop and /, shall be replaced by a proportionally — 


_ The first row of shear reinforcement closest to the column shall have cross 
Section area } according to Eq. 6, using v, ata critical section at d/2 from column 


r reinforcement 


ae a a 1. With shear reinforcement, the design of the critical section perpendicular > 

y 

worse 
q 

1 

\ 
| 

or v= — + 0.08 Vf". (MPa)... a 
iS tian iV the Siem Cross Section. the Tull Sirenein OF Lie 


elements as the first row; otherwise, the area may be reduced at a section 


farther away from the column face in which v, v, is smaller, but the requirement 


of Eq. 6 shall be satisfied at all sections. 


_ 5. The first row of shear reinforcement shall not be closer than 2 4 from 
column face when 4 V f” (psi) <y (psi) [0.33 St’. (MPa) 
= 0.50 Vf. (MPa)] in v, is calculated at a section d/2 
. _ from the column face, the first row shall not be farther than d/2 from the 
face and s shall not exceed 3d/4. When 6v (psi) <v,/o= 8 Vf". (psi) 

[0.50 Vf”, (MPa) < v,/¢ =< 0.67 V (MPa) ], the first row shall not be farther 

than 0.35d from the column fee and s shall not exceed d/2. With rectangular — “i 
columns (or with a square equivalent it in area to a circular column), the a 


4 and Codes (2,3) i in the 


Higher nominal shear stress, v,, is allowed a 

reinforcementis used: (8 V f ) instead of 6 V (psi) 0.67 Vf 


“= In determining the zone to te « covered by : shear reinforcement, the nominal 
shear stress, vw’, that can be carried by concrete in the absence of shear - 
reinforcement is reduced gradually as the critical section is moved away from 
the column | face > (Fig. 10, Eqs. 3 or 4). When the critical section is at a distance 
di or “closer from the column face, the value of v’ given by Fig. _ 10 is the 
Same as in the current codes. 
3. The cross section area of shear reinforcement is designed to resist shear 
in excess of a value v. higher than in the current codes Q, 3). For example, 
- for a square column at a critical section d/2 away from ‘the column, Eqs. 


PALE. 
4or = 3 Vi. (psi) (= (= 0. 25 ) instead of (psi) (0. 7 


of the po tests which i in fact indicate ‘that larger values of 


be assumed to be carried by concrete. > 


It is to be noted that Eq. 6 gives the s area of shear reinforcement per unit — 
length of the perimeter of the critical section rather than the total area required 


for the whole 2 section as in the ct current codes 2, 3). ‘The writers 3 prefer t to iw 


the case » when moment is transferred between the slab and the column. In 
_ this case, the nominal shear stress v, to be used in Eq. 6 is that caused by 
the combined effect of the factored vertical shearing force and the “* unbalanced” 
- moment transferred by eccentricity of shear in accordance with the current - 
- eodes (Section 11.12.2.2 of Ref. 2 or Section 9.14.2 of Ref. 3). 
When no moment is transferred between the slab and the column, the foregoing 


‘dain equations can be put in similar form to the equations of the codes (2,3) 
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| 
| | 
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with | by Eqs. 4 or 5. ‘The nominal strength Provided by 


in which = 


For vertios shear in concrete flat plates” must 


be well anchored at the two ends as close as possible to the top and — 
7 surfaces of the slab. For practical considerations, the reinforcement should 
be easy to install even in thin slabs and should not interfere with the placing | 
- of the flexural reinforcement. A new type of shear reinforcement | which satisfies 


“4 
these r requirements is presented. me 
_ The main results of several test series, in which the new seinforcoment is 7 
‘used, justify the suggested revisions to the current codes. The revisions allow — 
_ higher shear stresses in slabs with well anchored shear reinforcement; thus, _ 
thinner slabs can 1 be built. ‘The equations Suggested for design 1 result in a smaller ; 


portion of the shear stress can be carried by concrete. Rules are 
for placement and spacing of shear reinforcement. Setetits 
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in laboratory, The of theese tests indicate that 
symbols are used in this paper: finns. ip 


= 


= specified compressive strength of concrete, in megapascals 
cy the Use ribs, xad special ~ 
= specified yield strength of reinforcement; 
a § = spacing between rows of shear reinforcement elements measured — 
perpendicular to column face; Mewura) failures aed 
‘ (a V = applied shear force in test; ee 
= nominal shear strength provided by concrete; 
‘= nominal shear strength (force); 
nominal shear strength provided by shear reinforcement; 
V, = ultimate shear strength in test, or shear force due to factored load» 
a = assumed nominal shear stress resisted by concrete at ultimate to be 
in design for calculation of required shear reinforcement; 
ve = me nominal shear stress that can be resisted byc concrete at critical section 
outside shear reinforced zc zone; 
o =. = nominal shear stress at ultimate; i 
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a = = ratio of distance between column face and critical section to ob a > 
B. = ratio of long side to shorter side of rectangular colt column; and 
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Le Loaps ON Compostre Deck- REINFORCED 
foot spans ie om By W. Roeder,’ M. ASCE 


A Composite; po reinforced floors, of the type shown in Fig. 1, are € commonly 
used in modern construction. Moreover, this construction is economical. These 

_ decks are typically designed for distributed loads. _ Various studies (1,2,5,6,7) _ 

: 4 have investigated the behavior of floors under this loading. However, they are 
"frequently loaded with concentrated loads such as wheels of forklifts, or 


_ machinery. The behavior of these elements under this = is not well 


This paper describes ar 
‘floor system under concentrated loads. The study was motivated by observed E 
damage in an existing building on a U.S. Government facility. . The building 
was inspected and a series of nondestructive load tests were performed. Test 


- specimens were designed to simulate this floor yor system and they were tested 
to failure in a laboratory. The results of these tests indicate that the continuity » 
of the floor system restrains and limits deflections. This results in a greater — 
‘ capacity for resisting concentrated loads. This capacity may be ‘several —, 


larger than that predicted by existing design procedures. ME ar 


_baol bailqge Yo ads aiqoloveb oala yam 
_ The behavior of individual components of composite form-reinforced decks 
7 “4 are generally well understood (1,2). Composite beams are known to be influenced 
4 by the size and geometry of the ribs, and special design provisions have been ; 
propos posed for these beams. Other “studies (5,6,7) have considered the behavior 
of slabs when. they ar are isolated from the beam. Shear bond failure (7) is recognized — 
_as the most common mode of Fannie, but flexural failures and punching shear oi, 
These existing results are valuable but they do not consider the full complexity — 
of the composite floor system. Complex states of stress are developed in the a 


"Assoc. Prof., Dept. of Civ. Engrg., Univ. of Washington, Seattle, Wash. 98195. 
Note. —Discussion open until May 1, . 1982. To extend the closing date one ais, 
7 - awritten request must be filed with the Manager of Technical and Professional Publications, — 
ASCE. Manuscript was submitted for review for possible publication on October _ 
1980. This paper is part of the Journal of the Structural Division, Proceedings of the 
American Society of Civil Engineers, @ASCE, Vol. 107, No. ST12, December, 1981. 
ISSN 0044-8001 /81 /0012-2421 / $01.00. Goal “ite 
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we Le, aie cadionab ab 
slab. Compressive, bending, and shear stresses are developed in the slab due — 
to bending of the slab and the composite beams. Local stresses and punching 
shear stress may also develop in the vicinity of the applied load. wed by the 
and felatively complex boundary conditions, which are introduced by the © 


7 to form a more > complex: structural system than has ‘been studied in agen dl 
_tesearch. This study considers this more complex behavior, 


com were spaced 5 ft (1.525 m) on center, and they spanned between — 
frames which were 28-32 ft (8.54-9.76 m) apart. The = beams were composite 
beams with 3/4 x 3 in. (19 x 76 mm) metal studs and simple Type 2 (3) 
end connections. The slab was 4 in. (102 mm) thick (maximum thickness), with | 
6 x 6 in. welded wire mesh with #10 wire. The metal deck was 20 gage steel 
in 36 in. 1. 0. 915 m) widths, and it was continuous over 4 span lengths. The 
a ribs were 1.5 in. (38 mm) deep with inclined ‘embossments, as indicated pl 


— Spee. ih 
> 
) 


groups of test specimens were so designed to simulate the above floor 
S&S (Fig. 2). The same metal deck, wire mesh, and headed metal studs 
Ss used. The studs were welded through the metal deck into the flange of — 
a W 8 x 48 beam. This beam size was chosen to satisfy similitude requirements. © 
The normal weight concrete was mixed in a single batch and was designed 
have a 2,500-3,000 psi (17,200-20,700 kPa) seven day strength. > daw 
Groups A and B used single panel widths, four span sections of metal deck. 
4 They were identical except that the metal deck and wire mesh were continuous’ 
for four spans in Group B, and spliced over each beam in Group A. Each 
"span was tested independently by a procedure analogous to a shear bond test 
(5). However, the specimens were different from typical shear bond specimens 


METAL DECK IS 
x4 


: 


nese of the metal studs and the continuity of te slab. The fourth : specimen 
_ of both groups (A-4 and B-4) had an 8 in. (203 mm) diam core cut through 
the | slab and metal deck after the concrete had set. - Specimen C Cc was as designed 


_—Sketch of Test Specimens 
= 


of this continuity on the behavior of the structure. The slab was 1s designed to. - 
- be a full scale model of the prototype structure, but the beams were 0.5-0. 33 n 


scale because of space limitations, WE Tid 
xg All loads were applied with a a 2, 500 kip (11, 125 kN) Baldwin Hydraulic Testing © 


aa at midspan through 2 a steel wide flange. Six dial gages were attached 
to these specimens to measure slab deflection. The loads were applied monotoni- 


¥ 
— 
— 
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we 


cally, except that a cyclic uaiien was used for a A-3 3 and B B-3. The 
_ cyclic loading was used to reproduce the effect of repeated loadings on a slab 
. _which had sustained previous damage. It was applied by first loading the specimen 
to a load just above the load of visible cracking in the slab. The specimen 
was then unloaded, and the load was repeated until failure. Specimen C was 
loaded with concentrated loads which were applied to the slab at various locations 
through steel plates of various sizes (3 x 3, 5 x 5, 10 x 10, and 15 x 15 
in.) (76 x 76, 127 x 127, 254 x 254, and 381 x 381 mm, respectively). jae 
thickness of the steel plates were chosen to minimize deformation of the sl 
- and to provide a uniform distribution of load over the applied area. = ; 
gages and dial gages were attached to the flanges of the beams, and at | ft 
305 m) intervals on the midspan of the slab Cc. 
@ Test Resut 


7 ‘shown in Fig. 3. All eight specimens exhibited similar behavior. “The o center 
: line —T of the slab increased linearly up to the cracking shear, V,, and 
deflection, . At this load, - visible: cracking progressed through the depth of — 
the slab, as . iacieiad in Fig. 4. In most shear bond tests (7), this cracking 
would be accompanied by immediate slippage between the concrete and the 


moet deck, and a reduction in load capacity. However, slippage and p< 


thousand 
pounds per pounds | per ‘pounds per | 
Square foot width | foot width inh inches in inches 
“pascals) per per meter) meters) meters) 


(22,500) | Qs) 


(27,000) | Q5) 
4,050 16 
4.200 
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TABLE 1.—General Test Results for Groups A and B 
 Speci- 
q men 


slab. Thus, a trussing action in Fig. 4. This action produced 
: 4§ significant increase in shear capacity at reduced stiffness until the ultimate y 


shear and deflection, V, and A, , were attained. At the u ultimate load, the concrete 
began to crack around the metal studs and slippage and uplift occurred. A 
Teduction i in load occurred after the slip. 
| A comparison of the results for Groups A and B provides a measure of | 
a the influence of the continuity of wire mesh and metal deck on the shear bond — 
iz . The comparison shows that the continuous metal deck pres an 
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> 3. —Typical For 


HIGH 

"STRESS _IN METAL 
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FIG. 4. 4.—Trussing Ac Action of 
-—aproxinat average 40-50% increase in shear bond capacity over the spliced do 
deck specimens of Group A. However, there was considerable scatter in the — 
“data. This increase is somewhat surprising, because continuity does not reduce 
maximum | shear forces. However, a visual inspection of the specimens at the a 
ultimate load provides a hypothesis for r explaining this phenomenon. Considerable __ 
_ ‘megative curvature was noted in the metal deck over the beams. The metal - 
__ studs connected the metal deck to the beam flange and the flanges added 
_considerabie rotational stiffness to the metal deck. The increased stiffness caused 
_ the metal deck to carry a larger portion of the shear force, as is predicted 


by beam (4). Thus, it is believed that this stiffening effect 


q 
4 


q 


broken at one corner prior to testing. The initial break extended down to the 


: metal studs and prevented the slip restraint and trussing action which was noted | 

_ in the other specimens. This led to a significant reduction in V, and V,, and 
provides further evidence of the importance of the slip restraint provided by © 
the studs and the blocking action of adjacent spans. Specimen B-3 was also © 
tested under cyclic loading. The cyclic loading produced some deterioration | 

the strength and stiffness of both A- 3 and B-3. However, the 


Progressed — in A-3. The more deterioration in A-3 was 


ome 


FIG. 5 of Grouse A A and B with Previous Test Results (5) 
to > the | initial damage noted for Suecinen A-3. With specimens A-4 and B- 4, 
a hole was cut through the floor in the zone of high shear stress. The holes 


= cut because it was suggested that this operation could cause premature 


separation and loss of shear bond capacity. Both specimens had lower strength 
_ relative to other members of their group (Fig. 5). However, further testing 


to properly assess the degree of this reduction. 


& 


‘The results of these tests on specimen groups A and B were also’ compared 
3 to previously reported shear bond tests without stud shear connectors (7). The 


comparison is shown in Fig. 5. The previous tests were on the same type metal 
deck as used in this study, but the deck was inverted. It is recognized that 
the results for the inverted deck may not be directly comparable to the the results ; 


a a was largely a function of the metal studs and the beam flange width and thickness. - 
| 
+ 
| 


‘than #0 a normal slab of the same thickness. ‘The shear embossments and bond 
characteristics are similar. Thus, the inverted deck is expected to have oa 
bending and shear stresses, and higher strength, if the end restraint conditions — 
are identical. This comparison shows that the continuity and restraint, which 
i owas provided by the studs and the continuous floor, significantly increased 
the capacity over those in prev previous tests (7). Increases of 50-100% were obtained 
_ with simple span metal deck and stud restraint, and adjacent span blockage. — 
_ Continuous decks with stud restraint provided increases of 100-150%. aattogte, 
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FIG. —Load- -Deflection Curves for Specimen C 

ss Specimen C was loaded with a series of point loads at different locations. 
Most of these tests were nondestructive. General conclusions which can be 
drawn from these tests are as follows. Maximum slab deflections occur near ao 
the applied load, but they are smoothly distributed over a considerable distance. - 
A study of these deformations indicated that a significant portion of the applied — ‘. 


load was transmitted across the weak direction of the slab. Only approx 50% " 


_ of the applied load was directly supported by the loaded deck panel. The remainder 
was distributed to adjacent panels. It should be noted that the deflections of | 
ae C were compared to those obtained in a load test of the prototype : 


 gmactane, and the maximum deflections and distribution of deflection were similar. 
ts Two locations of Specimen C were loaded to failure and these results are 


va of this study. However, the comparison of Fig. 5 is believed to be a conservative _ ’ 
ike ‘measure of the beneficial effect of stud restraint and blockage of adjacent spans. __ 
| 
| 
| 
| 
| 
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examined in detail. Both tests were the at the 
_ center of the beam span, and at the center of a deck panel. The deformation 
of the slab (deflection of the beam is removed) under the concentrated load ~ 
is shown in Fig. 6. No concrete cracks were visible prior to these tests. The 
- first load was applied through a 5 x 5 in. (127 x 127 mm) steel plate, and 
"the second load used a 3 xX 3 in. (76 X 76 mm) plate. Both tests —_ 
similar behavior. The force-deformation curve increased without incident > ; 
to the 20-24 kips (89-107 kN) load range. In this load range, the — 
emitted numerous cracking sounds, and the strains in the metal deck — i 
_ below the applied load showed a sharp increase. However, there was no ' —_ i 
_ cracking on the slab surface. This was believed to mark the initiation of the ; 
4 shear bond cracking in the panel immediately under the load, but complete 4 
shear bond failure could not occur because of the restraint provided by the 
metal studs and the continuous metal deck and slab. This observation provides 
further evidence of the ability of the slab to distribute the load in its weak 
_ direction since this load is approximately twice the crack initiation load for 
a single panel (6 x V,)asdetermined in GroupB. 
_ The load capacity continued to increase beyond the crack initiation range - 
because of the restraint provided by the continuous metal deck and the metal — 
‘studs. Punching shear failure occurred at 37.4 and 31.5 kips (166.4 and 140.2 
KN), respectively, for the 5 and 3 in. (127 and 76 mm) plates. _ This mode of 
failure was very abrupt. The load dropped off sharply to about 20% of its 
‘maximum. This small residual load was carried entirely by the metal deck. 
There was no visible cracking in the concrete prior to failure. The 37.4 kip 
load was nearly twice the ultimate shear bond ‘capacity of a single deck panel a 


is distributed across the slab to adjacent deck panels. This indicates much larger r 
_ concentrated load capacity than aagane by typical one way slab design methods. 


Summary “AND Conciusions 
Metal deck composite are re practical economicel. 

_ This study shows that these floors may have a much larger capacity for resisting 

- concentrated loads than suggested by existing design procedures. There are -, 
two sources of this additional strength. First, the continuity of the metal deck 

_ and the slab, and the slip restraint of the ‘metal studs, increases the shear bond + 
capacity. Second, the continuity of the slab restricts the deformation of the 

- floor system. This ‘tends to distribute the deflections, bending stress, and shear — 
stress across the weak direction of the floor to adjacent deck panels. 
_ While these tests clearly indicate that slabs can resist large concentrated loads, a 
the results are not readily usable by the design profession. These tests were ; 
performed with a single metal deck type, stud size, span length, and concrete 


depth. Further study is needed to evaluate the effect of these parameters = 


& develop strength predictions which can be readily used by the engineering 

Most of the experimental work, which is described in this pa paper, was funded — a 


| 
> 


by the U.S. Navy through Centrac Engineering Associates in Lynnwood, Wash. a 
N. Bass and A. K. Batra served as contacts with these agencies and — 
numerous suggestions. The writer also thanks N. M. Hawkins for his continuous 
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Note 
To a place within ASCE for publication of rechaical 


as yet, to the point where they warrant publication as 
publication of Technical Notes was ‘authorized by beads Board of —— on October 16-18, j 


1. An original manuscript and two copies are to be submitted to the a of Technical 
, 10017, 


4. The technical publications staff will prepare he ee for use in the earliest pothte eel 
issue of the Journal, after proper coordination with the author. 7 pee ge 
5. Each Technical Note is not to exceed 4 pages in the Journal. As an approximation, — 
each full manuscript page of text, tables, or figures is the equivalent of one-half a Journal 

Journal. 


8. The final date on which a Discussion should reach the Society is given as a footnote 
_ 9. Technical Notes will not be included in Saimin. 
10. Technical Notes will be included in ASCE’s annual eat cumulative subject and ll 


The manuscripts for Technical Notes mu: must n meet the 
4 Titles must have a length not exceeding 50 characters and spaces. 


2. The author’s full name, Society membership grade, and a footnote teees Stating present o q 
oer. must appear on the first page of the manuscript. Authors need not be ie 


3. The manuscript is to be submitted as an original copy (with two duplicates) that is typed ay 

double-spaced on one side of 8-1 /2-in. (220-mm) by 11-in. (280-mm) white bond paper. sia a » 
_ 4, All mathematics must be typewritten and special symbols must be properly identified. 


The letter symbols used must be defined where they first appear, in figures or text, and arranged 


alphabetically i in an Appendix.—Notation. 
_ 5. Standard definitions and symbols must be used. Reference must be made to the lists 
_ published by the American National Standards Institute and to the Authors’ Guide to the 

6. Tables must be typed double-spaced (an original ribbon copy and two duplicate copies) 
a one side of 8-1/2-in. (220-mm) by I1-in. (280-mm) paper. An explanation of each table 
must appear in the text. in back ink on one side of #-1/24 

7. must be drawn in black ink on one side of 1/2- -in. by 11-in. (280-mm)_ 


: a 9. Dual units, i.e., U.S. Customary followed by SI (International System) units in parentheses, 


should be used throughout the paper. 
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Orrimum OF I- BEAMS 
Abul K. Azad," M. 
> 


In designing three- -span 1 I- beams or om designers may encounter si situatio 
_ where the locations of the interior supports are not prefixed and can be varied, 
_if advantageous. A common example is design of three-span highway overpasses 
_ where the location of piers can sometimes be adjusted to produce an economical — : 
a design. In order to find the most economical design from the viewpoint of 
minimum the optimum locations of the | interior supports must be 


6, 7,9) are covered in Ref. 7 | which also provides a comprehensive | list. of references ; 
_ This paper presents a a plastic design ‘method to determine the minimum-weight _ 
design that is of a symmetrical three-span I-beam of given length and is subjected 
to uniformly distributed load. The beam is symmetrical in cross section, and a 
has an unchanged web of constant height. Flange areas are varied considering 
only the : maximum positive and negative moments, thus, introducing = one a 


restricted number of cutoff points for relatively short span beams is not uncommon c 
since too many cutoff would the cost of 


Le 


F 
plastic section can be written 


2A, + ht 


"Assoc. Prof. of Civ. Engrg., Univ. of Petroleum and Minerals, Dhahran, Saudi Arabia. — : 
Note. —Discussion open until May 1, 1982. To extend the closing date one month, 
= written request must be filed with the Manager of Technical and Professional Publications, 
¥ ASCE. Manuscript was submitted for review for possible publication on January 21, 4 
7 a 1981. This paper is part of the Journal of the Structural Division, Proceedings o of the J 
American Society of Civil Engineers, ©ASCE, Vol. 107, No. ST12, December, , 1981. 
ISSN 0044-8001 /81 /0012-2433 / $01.00. ! 


i 
the flange plates are varied to save material. Considerable research has been [im 
= jj. done in the area of optimization of I-shapes using various assumptions and — 
mnroaches Refs If) are a renrecsentative of such wark hased on 
_| 


tok 
in which F, = the yield steel. 13, ‘1-3, the ‘total 


_ in which p p= icon unit weight of steel; M >, 
length, 4; and 


Referring to the symmetrical I-beam of leng 
Yo edt o: sab 
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7 | 7 a The value of Z, required of a beam section to resist a design moment, M > if 
jo 
7 : = the design moment for the segment __ 
. s, n = 3. Eq. 4 shows that, fora | 
n 1 valu Sas. 7 
= 


q 


y, as defined 6, is a value to be 


0. 
FIG. 3.—Values of F.. and y Versus Values of B_ 

. 5-11, the: objective function, F = =M,, Z for ¢ all se: eral 

V 0.58? — 4y } 


For a given of B, setting aF/ay = 0, ‘the ‘condition for 


42 25p? — 
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For an assumed value of B, the value to produce Frain be ‘obtained 
from Eq. 14 by trial. The values of i’ and F.,,, for various values of B are— 


shown in Fig. 3. While the optimum value sche = and the corresponding 
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the support locations and the corresponding. support moment, and each | span , 
4 weight is dependent upon its span length. The optimum segment lengths for 
the end span agree with the results presented in Ref. 9. — ; - 
Baw, Optimum Web.—Once the optimum design moments and the corresponding — 
segment lengths are known, the next step is to find the optimum web for the - 
minimum weight design. The flange areas required for various segments to provide ¢ 
the bending strength can be determined from Eqs. 2-3 because A and ¢ are 
a. Since the web thickness must correspond to commercially available plate 

thickness, and it must provide the shear resistance of the beam, the oo 

- acceptable web thickness is selected which would be capable of carrying the 

shear force within the guidelines of the code practice. The procedure for such © 
a selection of web thickness would be similar to the method in Ref. 1. For 
a selected web thickness t, the optimum web height is obtained “it setting a ae 
This height is if it the design restrictions on h from the 
practical view point and on the ratio h/t. In the case in which there is no 
= height limitation, the optimum height to satisfy the ratio h/t = a 
a ‘is obtained by replacing t with h/k and setting 5 = 0, and 


The v: value of h heen! ‘Eq. 16 may have to be modified in design since the 
value of t from t = h/k may not yield a plate thickness commercially available. 
Knowing A, t, and optimum segment lengths, the of the beam 


can be finalized. 


A plastic design method to find the minimum-weight design of a symmetrical — 
thr 


ree-span continuous I-beam of given length is presented. The beam is ~~ og 


= on one cutoff point in each end span and two cutoff points in the center 
Span, the optimum ratio of the Conter span to the end apen is determined as 
Zz 25. The method provides the | optimum segment lengths, the a 
moment distribution, and the web height for the beam. 
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sti2 
Fin = 0.007124n t an economical design 
j q can be obtained within the range of 6 = I.2-1.3. While the exact relationship _ a 
Po between B and y is described by Eq. 14, Fig. 3 shows that the value of y 4 
pct linearly with the weisht of the heam denends upon 
| 
| 
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_AEROELASTIC STABILITY OF SUSPENSION BRIDGES 


te By ‘Augustin Petre! and Ramiro 0 Sofronie® 


_ The aim of this note is to develop a theory on the sinemaieaiie stability of 
7 7 suspension bridges that takes into account the drag effects of aerodynamically 
_ shaped decks (1,2 ,3). The theory is applied both to classical suspension bridge _ 
systems, and to recently developed prestressed suspension bridge systems (4) ot 
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Fig. 1). Practical methods to solve reliability the stability problems of suspension - 
Basic Assumptions. —The only two deformations considered significant in 
problems of the type treated herein are : bending and torsion a) ee. 
_ The elastic forces corresponding to the previously defined deformations wil 
be expressed by differential operators. Since the elastic lines are ae 
, the expressions of elastic forces are independent of one another, i.e. 
slastic coupling terms do not appear. 
_ ‘The structural damping forces are considered ‘according to Voigt’ s model (3). 


included in the expressions of elastic forces. The imaginary component of these se 

The inertia forces are generated by the mass of decks, cables, and hangers 

mada to be > uniformly distributed spanwise. Due to the vertical symmetry 


7 _of bridge cross-sections, the mass centers coincide with the elastic centers. — 

_ The aerodynamic forces are considered according to the unsteady i incompressi- 
be flow the serogynemic operators of type. 


uspension System 


7 assumption that : tol operators are applicable is justified by the increasing — 

use of streamlined suspension bridge decks, the aerodynamic — 
of which are similar to those of airfoils (6). flea om ; 
Ri To calculate the equilibrium stability between the previously described ai 

_ itis assumed that the bridge deck is subjected to an arbitrary but small perturbation 
consisting of a vertical displacement, w or a transverse rotation, , or both. 

; Consequently, under the wind action, aerodynamic forces of lift and drag develop. ; 


f From the interaction between these aerodynamic forces and the elastic restoring a 


forces there results a a motion ora deformation, according to whether the influence 


perturbation decrease and vanish, ‘the resulting aeroelastic phenomenon i is stable. 

_ On the contrary, if the initial perturbation increases, the phenomenon is unstable. 
; _ Between the two equilibrium states there is a limit state in which the aerodynamic _ 
- forces are perfectly balanced by the elastic forces. In this state, any additional 

perturbation, however small, does not encounter any opposition from the 


‘defines the onset of the seroelastic | phenomena of flutter or div ergence (1). ol 


« Governing Equations.—In the most general case, to the initial perturbation — 


| 
| 
if 
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T1200 
of bending and torsion previously assumed, there corresponds 1 the motion > al. } 
- = (Mo) = 7% . (la) 
a + ig.) Gl, = 


in which in each equation the first three terms represent elastic forces with 


_ structural damping, inertia forces, and lifting airloads. The last terms introduce 
the drag effects and, for a simply supported deck, ‘the bending moment M, 


ass sociated with the unifor uniform running drag has the expression den wad | 


in which V, ub 
For or simple harmonic” motion of circular 


If, in Eq. 1, the inertia forces are reel initial perturbation reduces to — 
abending-torsional deformation (7), 
Note that the two components of the motion ( 1) are always aerodynamically 
_ coupled. But as Eqs. 3a and 3b, associated with the boundary conditions coe | 
= = o and (0) = « =0 form Sturm systems, they can 


q 


w? mw(t) + — — de 


In these equations G and G are s influence functions 


= bending and torsion, respectively, of the considered suspension bridge ~ 


iw 


¢ 


_ The canonical expressions of aerodynamic force and its moment, with respect 7" 
to the elastic axis in the adopted strip-theory, are given in Ref. 7. 
“technique does not change w when other expre ions al are 
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order to the st stability equations are to 
be first nondimensionalized. The unknown parameters are ical 


and Z = (I + where 
‘for dynamic stability, and = 
2 


— — 


Classical Suspension Bridge System. “The special features wis this 

_ as treated herein, are the following: (1) The cables hang freely in vertical planes 
following the funicular lines of the loads applied on the deck; and (2) the 
cross-sections of the bridge system (which includes the cables) are open Fig. 
la). For this type of bridge with a deck of constant cross-section, Green’ Ss Pa 
influence functions, G,,(x, €) and G, (x, €), and integral Eqs. 4a and 4b of — 

stability, as well as those of static ‘stability, are givenin Ref.7. end 
Prestressed Suspension Bridge System.—The special features of this structure, , 

_ as treated herein, are the following: (1) The cables follow the geodetic a . 
of the funicular surface of vertical loads applied on the deck; and (2) the 


cross-sections of the bridge system (which includes the cables) are closed ( 10) 

: (Fig. 15). For this type of bridge, with a deck of constant cross-section, Green’s _ 
influence functions and the integral Eqs. 4 of dynamic Stability, 4a and 4b % 

as well as those of static stability, are given in Ref. Bein footer et) o} ig 
Solution Methods.— —The two static and dynamic problems o of aeroelastic stability 
previously described ar are, in fact, eigenvalue problems. Since by | using Green’ s 
influence functions the governing equations were put into integral form, it is 
; more practical to transform them into algebraic equations. The operation can > 


be done, for example, by punctual collocation, with the aid of the tangent 
‘method. Thus, integral equations 4a and 4b, for -aeroelastic flutter, take the 


0, 
(A, IG, + 116] 


i! 


| 
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in which the drag were ‘jatroduced wy the operators D,, and 

the ratio W,/0, contains the phase difference between the two components ; 

of motion; and, obviously, the normalized vectors, , [w] and [4], replace the =) 

_ complex mode shapes, w(x) and (x). Such algebraic. systems also could be 

directly obtained by a finite element technique. However, for aeroelastic stability | 

_ problems, the previously described method is much simpler and reliable. ae et 

_ In order to solve this linear system, two methods can be used, the direct | 
method and the iterative one. The well-known method of assumed modes is 

7 not so adequate in the case of suspension bridges because it is difficult to _ 

predict the common deformations. of the cables and decks. In the direct method, fl 


— — 


the matrix of a linear system is transformed into the Hessenberg matrix and © 
the eigenvectors and eigenvalues of this matrix are then determined (3). The 
method is useful when the explored field is completely unknown, although it 
eS more computing time, and, for plotting the V-g curves, each eigenvalue - 
has to be pursued. The iterative method is, in fact, the inverse power method. 
tis. is a much faster method but not always convergent enough, especially when 


the critical velocities of divergence and flutter, affected by drag, assume a 


= aa Both methods were used by the writers to investigate the aeroelastic ely 


4 of suspension bridges. For example, by asouming for dimensionless quantities 
(7) the values bee 0. 001, r? = 0. 25, m = 1, C = 0 (hence without diag), | 

p = |, and A = 10, according to Eq. 5, the dimensionless critical ‘aibeley 


a 


a takes for the classical funicular bridge, the value V, = 4.3447, and for o 
prestressed geodetic bridge, V, = = 13.673, ie. 3. 147 times larger. If, for the 

same parameters, ji increases from | to 10, the critical velocity of the funicular 
bridge becomes V, = 13.937, i.e. 3.208 times higher. But as a consequence — 

of aerodynamic drag, the critical weneny of the same funicular bridge decreases 

to V, = 7.244 for C = 0.02, and to V, = 1.7552 for C = 0.04, Le. ‘it eniaell 
14 and 7. 7.940 times lower, respectively. 


WP 1. The aeroelastic stability of suspension bridges is strongly influenced by | 
the aerodynamic drag of span decks. . By taking into account the effects of 
> 
drag in stability analyses, the two components of deformation or motion are J 
always coupled. A consequence of drag is the reduction of the absolute values 
_ of divergence and flutter critical velocities with respect to the values calculated : 
_ only for lifting forces. In addition, the ratio between the two critical velocities — 
also can be modified. Thus, if without drag the critical velocity of divergence | 
_ is higher than the corresponding critical velocity of flutter, as a consequence _ 
of the drag effect, the inequality sign can be reversed. 
_ 2. The governing equations, determined as presented herein, allow a compari- 
son between the two structural systems of suspension bridges dealt with in | 
this note. Thus, while their bending components are formatly identical, the 
are essentially different. Classical suspension bcidiges 


| 
| 
| 
il 
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more sensitive to torsion atu prestressed suspension bridges. For this reason, 2 
_ the effect of drag is greater on the former than on the latter. Also, while funicular — 
bridges behave quite differently in bending and torsion, geodetic bridges have 
a similar behavior in both degrees of freedom. 
3. The analysis technique | based on Green’s influence functions use¢ used herein 
is an efficient tool for the study of aeroelastic stability of the two types of _ 
suspension bridges considered herein. In addition, it allows one to discover — 7 
the aerodynamic and structural parameters responsible for any forms of static © 
or dynamic instability and to modify them conveniently. The same analysis — 
ee can be used for the aeroelastic optimization of suspension bridge — 
1. Petre, A., Teoria aeroelasticitatii, Statica, Ed. Academiei R.S. R., Bucuresti, Romania, 7 
2. Petre, A., Teoria aeroelasticitatii, Fenomene dinamice periodice, Editura Academiei_ 
By Petre, A., Ashley, H., ‘‘Drag Effects on Wing Flutter,” Journal of Aircraft, 
Sofronie, R., Suspension Bridge,” United States Patent No. 3,857,130. 
5. Petre, A., ‘‘On the deformation of one-dimensional structures,’ ” Rev. Roum. Sci. ~~ 
Techn. -Méc. Appl. Vol. 25, No. 1, pp. 85-98, Bucharest, Romania, 1980. 
6. Scanlan, R. H., Tomko, J. J., Airfoil and Bridge Deck Flutter Derivatives, _ Journal 
of the Engineering Mechanics Division, ASCE, Vol. 97, No. EM6, Dec., 1971, pp. 


7. Petre, A., Sofronie, R., “Aeroelastic | Stability of rll Bridges,’’ Wind Engi 
ing, Proceedings of the ‘Fifth International Conference, Fort Collins, Colorado, U. .S.A., 
July, 1979, J. E. Cermak, Ed., Pergamon Press, Oxford, U.K., 1980. 


_ 8. Sofronie, R., “‘Vertical Deflection of Suspension Bridges,”’ Rev. Roum. Sci. Techn 
_Méc. Appl., Vol. 24, No. 3, pp. 485-505, Bucarest, Romania, 1979. 
Sofronie, R., ‘Torsional Deflection of Suspension Bridges,’’ Rev. Sci. Techn.- 
_ Méc. Appl., Vol. 24, No. 4, pp. 647-664, Bucarest, Romania, 1979. Note ap 
10. Sofronie, R.. of Prestressed Bridges,” ” Rev. Roum. 
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REINFORCED CONCRETE BEAMS IN PURE’ Torsion 
By Himat T. Solanki," M. ASCE = 
_ Design provisions for torsion are now available in various codes (2,4,5), but 


the —_ differ considerably. Out of these approaches, space truss theory 


hater 


*=0.66+035 sus 
wed, 


FIG. 1. —Ettect of y/x on Ultimate Torque 


and skew bending theory pi predict | the true Nehanten of the reinforced concrete. 

beams. Past research (16) shows the comparison of torsional strength of reinforced . 

| a concrete beams with different theories (Fig. 1). By considering only the under- — 


| 
> 
3 


ULTIMAT 


reinforced beams with stirrup spacing, s < (x, + y,)/4, the American Concrete 


Institute (ACI) method has been studied in this paper and then modified. The 4 


ee. Engr., United Engrs. and Contractors, 30 S. 17th St., _ 15M3, ‘Philadelphia, Pa. 19101; < 
_ Formerly Chf. Design Engr., Transportation Dept., Engrg. Div., Sarasota, Fla. 2 
Note.—Discussion open until May 1, 1982. To extend the closing date one ents 
a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on October 18, 
1979. This paper is part of the Journal of the Structural Division, Proceedings of the 
American Society of Civil Engineers, ©ASCE, Vol. 107, No. ST12, December, 1981. _ 
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under- reinforced beams were determined by the Mitchell and 4 Collins approach 
ultimate torque of a rectangular beam can 
) 


in which T, = ultimate ceneies T.= = torque resisted by the Ahad at — 

and T, = torque resisted by the reinforcement at ultimate. "79 
using a ‘‘full’’ elastic stress distribution as proposed by Kuyt (1D), the 
torque resisted by concrete can be expressed as 


= cylinder compressive strength of concrete; smaller overall: 


sturreart 
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2.—Experimental a, as Function of y, 


éeeesicn ao of rectangular section; and y = larger overall dimension of rectangular 


According to the philosophy of ACI 318- -77, th the torque resisted by the 


which A, of one leg of a = yield s strength of 
oe 4s, enalles center to center dimension of a closed rectangular stirrup; 7 


ne larger center to center - dimension of a closed rectangular stirrup; s= 
“spacing of stirrups; anda, = acoefficient. 
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; 
For beams having an ‘volume of longitudinal and and ‘Stirrup steels, 
3 may y be mi multiplied by ry Vm, | in which m=o0f,/ef5 = ratio of yield force 
of longitudinal reinforcement to yield force of closed stirrups per unit volume. 
Substituting the value of T. and T, into Eq. 1, the ultimate torque of a onmeie 
reinforced concrete section unequal volume and stirrup 


to the ACI 318-77, the a, is a function 


fore, all the computed values of a, were plotted 2 as a function of y,/x, in 


Fig. 2 by using test and | 
Mr 


It can be seen that the is a of y,/x, and statistically 


1 the correlation between the strength predicted by Eq. 
7 4 and available test data is good. Only under- reinforced square or rectangular 


beams with s < (x, + y,)/4 are included in this analysis. When both longitudinal q 
and stirrup steels yield at the approach of failure of beam, it was classified 
as an under-reinforced beam. Table | includes hollow box beams (7,10,13,14,15, 
| 19), beams with lightweight concrete (21,26) and asymmetrical reinforced : 
4 beams (13,14,17). Based on Ref. 6, the value ine was determined for asymmetri- 
ss cal reinforced beams. For European beams, the cylinder compressive strength 
was chosen equal to 0.85 Feuve © based on actual comparison (9,24,25) instead 
(0.76 + 0.2 log 2840) and was chosen equal tof, 
Minimum Reinforcement.— —The criteria for minimum reinforcement is 


- Substituting for m in Eq. 8; assuming fy = f =f, and raranging 
|= — 0.25x) VF" fed ai 


- Maximum Reinforcement.—The criteria for maximum reinforcement are Per 
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1.—Correlation of Results Under- Reinforced Beams 


Source | 


= 


a 


wt 


Virginia 


West Virginia 


University of Washington 


University of 
University of 


q 
inkip- | 


(4946 | 0. 
06 
02. 


: 
Beam inches | 
| | 197 | 209] 
Gas 80.34 | 5.4 106 
D¢ =| 2046 | 0B 
| 
3 102.6 | — 
- 


3 University of Alberta 


=" 
Cement & Concrete Association ¢ rere 


Delft University of Technology 


Eidgenéssische Technische 


Hochschule pA olin 


Bending, and Shear” 


7B. 


University o of New w South Wales | 7&7A | 105 


ot of Ras Plata Age 


| 109.6 | 0.75 
ri? 1,270 1,173.8 | 
2° 1,266 1,185.2 T 
3 907 935.1} 0.97 
var | | 232] Om 
VQ2 im | 23.4 0.83 
| «(373 | 426.0 | 0.86 
10 
130 
= | | 


DECEMBER 
1.—Continued 


853 781.7 
058 


kip-in. = 0.113 KN-m. 


Asymmetrical concrete beams for beams -1, T-3 and T- -8 are 0.99, 


Taking 1/3 x (y 0.25 25x) 2 Vf‘. and using Eq. 


—x"(y—' x) +a, mA,x at 


for m i in ll, , assuming fy =f, and rearranging 


q 


reinforced concrete beam decreases with an increase in aspect ratio, y/x. 
ie 2. The coefficient a, varies with y,/x, and can be calculated by using Eq. 6. 
“4 Ea. satisfactorily predicts torsional moments of under-reinforced square 
Eq. 9 isa a simple “expression for minimum reinforcement while Eq. 
_ §. Because the beams with stirrup spacing s = (x, + y,)/4 are considered bs 
in Table | ‘it is necessary to increase the value of a, , compared to the ACI 
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Discussions may be submitted on any a paper or technical note published in any 


: ~ Journal or on any paper presented at any Specialty Conference or other meeting, the Proceedings — 
4 ‘ of which have been published by ASCE. Discussion of a paper/technical note is open to 


anyone who has significant comments or questions regarding the content of the paper / technical 
note. Discussions are accepted for a period of 4 months following the date of publication 
of a paper/technical note and they should be sent to the Manager of Technical and Professional 
Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. The discussion period may 
be extended by a written request from a discusser. _ eS ee nS 
The onginal and three copies of the Discussion should be submitted on 8-1 /2-in. (220-mm) _ 
by Il-in. (280-mm) white bond paper, typed double-spaced with wide margins. The length of | 
a Discussion is restricted to two Journal pages (about four typewritten double-spaced pages 
of manuscript including figures and tables); the editors will delete matter extraneous to the a 
_ subject under discussion. If a Discussion is over two pages long it will be returned for shortening. 


7 All Discussions will be reviewed by the editors and the Division's or Council's Publications 


” ienieals for Discussions are the same as those for Proceedings Papers. A .—— is 
_ subject to rejection if it contains matter readily found elsewhere, advocates special interests, 
_ is carelessly prepared, controverts established fact, is purely speculative, introduces personalities, 
_ or is foreign to the purposes of the Society. All Discussions should be written in the third 
7 person, and the discusser should use the term ‘“‘the writer’’ when referring to himself. The 
of the original paper/technical note is referred to as ‘‘the author.” 

_ Discussions have a specific format. The title of the original paper /technical not - appears 7 > 
5 at the top of the first page with a superscript that corresponds to a footnote indicating dl 
month, | year, author(s), and number of the original paper/technical note. The discusser’s full » 
2 should be indicated below the title (see Discussions herein as an example) together with — 


his ASCE membership grade (if applicable). 
‘The discusser’s title, company affiliation, and business address should appear on the first ; 
page of the manuscript, along with the Proceedings paper number of the original paper /technical 
note, the date and name of the Journal in which it appeared, and the original author’s name. 


_ Note that the discusser’s identification footnote should follow consecutively from the a 


-paper/technical note. If the paper/technical note under discussion contained footnote aumbers = 
. and 2, the first Discussion would begin with footnote 3, and “Subsequent Discussions would — 


_ Figures supplied by ‘the discusser should be cide’ by letters, starting with A. This also 


appeared in the original paper/technical note use the same number used in the original. _ 
It is suggested that potential discussers request a copy of the ASCE Authors’ Guide to 
the Publications of ASCE for more detailed information on preparation and submission of 
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Seismic Desicn beut FOR SeconpaRY Systems" 
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- ‘The writer would like to thank ‘ew and O’ Connor for their interest in 
the paper, and for their thoughtful and interesting comments. Indeed, if design — 

: ground motions were defined as random processes and characterized in some 
appropriate form, the need for direct generation of floor spectra approaches ; 
would not arise. However, the definition of design earthquakes by some 
characteristics of a random process is not a straightforward and easy task. 
The definition in terms of a (design) power spectral density function (PSDP) 
is rather very tempting, a and has in fact been used byt many investigators including — 
the writer to study dynamic behavior of structures under earthquake type of ‘ 
excitation. ‘This form of characterization is preferred because analytical 

7 3 manipulations are relatively simple and manageable. Mathematically speaking, | 
_ however, the earthquake motions cannot be characterized by a PSDF as they — 
are not stochastically stationary i in 1 character. _ Then, does ik mean that many 


‘assumption that a PSDF exists are e all ‘wrong? Probably ‘not ‘entirely, and at 


least not qualitatively. Quantitatively, however, there is room for Saee 


a A realistic earthquake motion model must, of course, be nonstationary. But we 


difficult to answer. Stationary models modulated by deterministic time functions © 
- (also called envelope functions) have often been used as nonstationary models 
a for earthquakes. These modulation functions are generally assumed to start 
with a motion buildup phase followed by a strong motion (stationary) phase 
_and then by a decaying phase. Other forms of modulation functions have also 
been used. However, to prescribe models for designs, a more precise definition 
of the modulation function parameters an and PSDF will required. To obtain 
_ these parameters directly by a statistical analysis of available ground records _ 
is at best a delicate task. Whereas to obtain these parameters from a given 
i spectrum and then use them as input in the analysis is quite unnecessary especially 
when a spectrum can be directly used in a so-called single-step direct approach. — 
_ The discussers refer to the topic of composite damping and modelling uncertain- — 
ties. In a direct approach, these problems also don’t seem insurmountable (15,16). s 
, ‘The writer, however, fully supports the views of the discussers that the research 
"must continue to define a realistic stochastic model for earthquake ground motions. j : 
‘February, 1980, by Mahendra P. Singh (Proc. Paper 15207) 
* Assoc. Prof., Dept. of Engrg. Sci. and Mechanics, Virginia Polytechnic “kre 7 and 


State Univ., Blacksburg, va. 24061. Te: 
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“15. Singh, M. P., “Seismic Response by SRSS for Nonproportional Damping,” Journal J 
of the Engineering Mechanics Division, ASCE, Vol. 106, No. EM6, Proc. Paper 15948, 
16. ‘Singh, M. P., “‘Seismic Response of Structures with Random Parameters,”’ 
_ at the Sept. 8-13, 1980, 7th World Conference on Earthquake Engineering, held at 
Istanbul, Turkey, 


The should be made to ‘the origi ‘nal 


‘Page 2259, line 5: Should read “‘relief’’ instead of ‘‘refief’ 


~~ paper by mR provides a comprehensive summary of the next — 
g of critical problems that must be closely defined and solved in order to 
' improve the present state of understanding of the stability of metal members 
and frames. . Task Group 4, ‘‘Frame Stability and Columns as Frame Members,” 
of the Structural Stability Research Council has attempted, over the past several © 
- years, to collect, discuss, and critically examine the research related to the 
stability of steel building frames and the writer has been requested to communicate, — 


on behalf of Task _ Group + 4, to discuss the progress and future directions sin 


_ As outlined by Bjorhovde, the general thrust of the Task Group’s effort - 


has been to encourage studies and assess results that will be useful in moving — 
(November, 1980, by Dennis Clark (Proc. Paper 15794). 
ba December, 1980, by Reidar Bjorhovde (Proc. Paper 15926). 
"Prof. of Engrg., Univ. of Alberta, Edmonton, Canada. 
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DISCUSSION 
the eventually from the “K factor” approach to design to the pind 
Delta’”’ appro 
used to more readily reflect those factors and influences ‘that have ; a significant : 
effect on the behavior of the building frame. In fact, the design standard adopted 
Z by the Canadian Standards Association has been written to encourage the use 
_ of the P-Delta approach, with encouraging results. (Canadian Standasds J Associa- 
tion, ‘‘Steel Structures for Buildings—Limit States Design,’””» CAN 3— ~$16.1— 


_ The Task roid recognizes, however, that while the supporting g evidence > 


‘p- Delta me method. "Specifically, since an elastic de is used to predict the 

4 story sway deflections at the factored load level (implied values: of these 

deflections are used even in an allowable stress design) it is important to know 

a - the extent of the underestimation of these deflections and the pce oder ol 

Stability effects. Thus, the parametric studies referred to dsl Bjorhovde are 


Secondly, it is important to be able to the “‘P-Delta’”’ approach to 

- different structual systems. This implies the ability to predict the story sway 72) 
deflections of such systems. Bjorhovde has mentioned studies of semirigid 
systems and infilled frames. In addition, the stiffening influence of standard 
. connections, the member width effect, etc., provide other areas where research ~ 


information can be used directly in the design method. 
a. And finally, it is essential to come to an understanding of how "secondary 
5 ef fects other than P-Delta influence the behavior « of the . structure. How important — : 
is the influence of axial load on member stiffnesses? Are there practical systems 
where these member effects are significant? At the present time it appears 
_ that these cases will be rare indeed but as a part of the continuing study of 
o behavior of steel frames this important question should be addressed. epee 
_ Efforts to improve the understanding of the behavior of steel Structures will 
never be completed. In the present period the Task Group has also attempted : 
- to encourage studies that will assist in the application of the P-Delta approach ‘ 
to particular classes of structures. One such study is sponsored by the American ; 
Institute of Steel Construction and has consisted of an extensive literature review — 
and computer study. The study will recommend the adoption of the P-Delta 
for inclusion in design “specifications and will focus on the factors 
, required for a conservative application of the approach to both the combined 
loading case and the gravity load only case. GA, 


_ The writers would like to commend the author for bringing out an extensive | 

; report on the research needs in stability of metal structures. The lateral torsional : 
*Full Prof., Fachgebiet Stahlbau, Technical Univ. of Berlin, Berlin, West Germany. 


° Alexander von Humboldt Fellow, Fachgebiet Stahlbau, Technical Univ. of Berlin, 


_ the ‘‘K-Factor’’ approach, additional information is required relating to the jf 
| 
— 
ig | 
i 
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buckling of beams are important in construction but 
during the life of the structure (69), 
- Tes author has listed: (1) Residual stresses on beams; (2) effect of ee 
stress on inelastic lateral torsional buckling strength of beams; (3) parametric 
_ Studies on effect of residual stress and initial out-of-straightness; (4) laterally — 


_ supported continuous beams; (5) effect of cross-sectional distorsions; and (6) 


effect of coped ends as research needs for laterally unsupported beams. The a 
writers are of the opinion that the following should also be included as research ~ 


needs for the case of on tent 


we? Effect of shape beams: Most : of ‘the been done is 
& on symmetrical or unsymmetrical I- Beams. Therefore research is warranted 


in channel and Z Sections, which are quite often used as purlins,§ = 
Effect of sheeting on lateral torsional buckling: It has been shown that 
the diaphragm bracing, if properly utilized, can eliminate the need of other 
7 types of bracing (59). The lateral al torsional | restraint by the roof or f floor decking a 
depends on a number of factors. Much research has to be done to include 


the restraint provided by the sheeting, which will provide economical solutions 


3. Intermediate restraints: The intermediate restraints substantially | increase 


“the critical moments on beams. The effect of these intermediate restraints are 
equally important as that of the effect of end restraints i 
i... The effect of bottom flange restraint : Very little is known about the effect 
of bottom flange restraints (57), though certain amount of literature is available 


on the effect of mixed end restraints (67), 


5. Effect of openings on lateral buckling of beams: For aesthetic appearance 


and to save e material, castellated beams : are being increasingly used. Practically 


lateral stability of: a slender -castellated beam (58,61). 


‘The design formula recommended i European Convention for Constructional _ 


steel work (ECCS) for unsupported beams is given (64) | 


M, Mg) “M, = = elastic critical moments; and n = = system factor. The 
recommended a single value of 2,5 to ensure that the resulting design 
— to a meanvalue for rolled beams. The choice of one curve charac- 
terized by n = 2,5 is mainly controverted by the Japanese researchers, who 
_ propose to use several design curves with respect to the type of beams: rolled, 
i welded or annealed. They have made tests on a lot of welded beams, - which 
show that the experimental results may be below the ECCS design curve (60,65). 
Thus, many more tests results or computer simulations (or both) are — 
to fix the system factor n (Ref. 33 of the original text, and =. 


pe 
a 
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J | the load should be studied for the preceding cases. J 
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Prorosep SPECIFICATIONS FoR STEEL Box GIRDER BripGes" 
Discussion by Charles Seim,* F. ASCE and Steve Thoman® 


The Author’s proposed specification advances the art of bridge engineering — 
_ and fills in a gap in the design of long : span steel box > a >= In the 


a “December, 1980, by Roman Wolchuck (Proc. Paper 15942). 
“Senior Vice-Pres., T. Y. Lin International, San Francisco, Calif. cai 


a Sal Engr., T. Y. Lin International, San Francisco, Calif. oot 
= 


q 
| 
| 
| 
| 
= 
] 


writers’ opinion, at the time, steel are not competitive for 
spans shorter than about 500 ft with a concrete deck and possibly 750 ft Q29 
"7 m) with a steel orthotropic deck. Whereas a prestressed concrete bridge of 
: these span lengths requires a box girder configuration. Many times engineers 
or owners compare the cost of steel box girders with concrete box girders 
because of their similarity in appearance. However, if the comparison is —— 
on the type of bridge considering only its economic span range, the steel plate F 
gitder o1 or steel truss type may be more competitive with the concrete box. 


The writers consulted on the designers of a composite steel box girder bei bridge 


having five continuous spans of 310 ft, 400 ft, 450 ft, 400 ft, and 310 ft. The 
trapezoidal box section consisted of three webs having a maximum depth of © 
21.4 ft and a minimum depth of 10.4 ft. The report “Design Specifications 
_ for Steel Box Girders’’ submitted to the FHWA by the firm of Wolchuck and 
_ Mayrbaurl was utilized as the basis of the design criteria for the webs and 
- flanges. It became apparent that in the hands of a designer not familiar with 
buckling design or behavior that the proposed specification required effort 
considerable beyond the techniques employed for the design of other elements - 
_of the bridge. A shortened version of the proposed ae was abstracted 
to simplify and ensure a more manageable design, 
_ Perhaps the final version of the specification could be divided into a shorter — 
simplified section and a section in the more complex form as the proposed _ 


‘Specification now W exists. Thus, the designer cc could choose the degree ys refinement - 


4 


types rather than by sefinemem of stiffener arrangement for a steel box. 


- pr ertenersane the total strength of a web panel as sum of the beam Puce strength | 
and the tension field strength. The rigidity required of the girder flanges to 
_ adequately anchor the tension field forces has not been extensively test proven. 
Therefore the proposed specification makes the assumption of neglecting the 
flange” Tigidity, corresponding to a lower limit of the tension field strength. : 
A simplification can be made by ignoring g the tension field strength which greatly 
_ reduces the design complexity without greatly increasing costs. 
_ Since the proposed specification has not had the test of application by skilled | 
designers to a variety of span lengths and box sizes, the specifications should — 
be used with discretion until these parametor studies are completed and the _ 
Bay has accepted the specification as a new design tool. The fl 
has made an excellent start and the writers hope that some form of the proposed 
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EpGE on FOR Memsers* 
_ Discussion by Srinivasan Sridharan* 


so _ The authors are’ to be congratulated for a very interesting and valuable paper 
on the subject. There appear, however, some questions on their 
"analysis. ” The computer approach (8) they have adopted does not appear to ff 
have taken into account the nonlinear contribution of in-plane displacement — &g 
5 of plate elements in the strain displacement relations, very vital in modelling | 1 
: the type of buckling and postbuckling behavior considered. ‘This comment is 
8 of particular significance regarding the lip of the channel which does suffer » 
large scale in-plane displacements. This difficulty is reflected in Fig. 6, which " 
> indicates very nearly the same buckling stresses for a spectrum of d/t ratios 
 (d/t = 18.2, 15, 9.9, and 8.0), a trend clearly contradicted by Fig.4. = 
The authors do not mention the boundary conditions and the mode of application — 
of load employed in the analysis. Was the load applied by prescribing uniform y 
end shortening to simulate the experimental set up? a ae 
; ‘a Finally, what was the nature of failure? For the proportions of lipped channels ~ 
— tested (w/t = 45.6, Fig. 11), one would expect failure by an interaction of 
initial imperfections and plastic yielding. If so, would they comment on the 


value of a purely elastic postbuckling analysis? = ne 
a ~ “February, 1981, by Thomas P. Desmond, Teoman Pekoz, and George Winter (Proc. ; 


“Asst. Prof., Dept. of Civ. Engrg., ‘Washington Univ. in St. Louis, Louis, Mo. 
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Original papers should. be submitted in triplicate to o the Manager of Technical and Sueded 


Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. Authors must indicate the _ 
Technical Division or Council, Technical Committee, Subcommittee, and Task Committee (if any) 
to which the paper should be referred. Those who are planning to submit material will expedite — 
the review and publication procedures by complying with the following basic —— 
Titles must have a length not exceeding 50 characters and spaces. 


; 2 The manuscript | (an original ribbon copy and two duplicate copies) ‘should be. double-spac -spaced 
on one side of 8-1/2-in. (220-mm) by 11-in. (280-mm) Paper. Three copies | of all couse and 


tables must be included. 


- Generally, the maximum length of a paper is 10,000 word-equivalents. As an approximation, 
each full manuscript page of text, tables or figures is the equivalent of 300 words. If a particular 
subject cannot be adequately presented within the 10,000-word limit, the paper should be accompanied 
by a rationale for the overlength. This will permit rapid review and approval by the Division 
or Council Publications and Executive Committees and the Society's Committee on Publications. 
= contributions to the Society’s publications are not intended to be discouraged by this 


- aa The author’s full name, Society membership grade, and a footnote stating present sdidinaaaiae 
must appear on the first page of the paper. Authors need not be Society members. | 7 Wee 


aq 
 §, All mathematics must be typewritten and special symbols must be identified properly. The — 
letter symbols used should be defined where they first appear, in figures, tables, or text, and - 
arranged alphabetically in an appendix at the end al the paper titled Appendix. --Notation. = ~ 


oe: Standard definitions and s symbols should be used. ‘Reference should be made to the lists = 
published by the American National Standards Institute and to the Authors’ Guide to the Publications 


me Figuees should be a. in black ink, at a size that, with a 50% reduction, would have 

a published width in the Journals of from 3 in. (76 mm) to 4-1/2 in. (110 mm). The lettering 
must be legible at the reduced size. Photographs should be submitted as glossy prints. Explanations _ 
and descriptions must be placed in text rather than within the figure. Se ie ake 
8. Tables should be typed (an original ribbon copy and two duplicates) on one side of 8-1 /2-i rin. 
20am) by | 1-in. (280-mm) paper. An explanation of each table must appear in the text. os 


9. References cited in text should be arranged in alphabetical order in an sgquedin at the © 
ow of the paper, or preceding the Appendix. —Notation, as an Appendix.—References. i ¥ 


a 10. A list of key words and an information | retrieval abstract of 175 words should be provided 


12. A set of conclusions must end the paper. 


13. Dual units, i.e., U.S. “Customary followed by ( aternational System) units in parentheses, 
should be used the Paper. 


14. A applications section included also, if appropriate. 
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